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EVALUATION  OF  EXPLOSIVES  STORAGE  CRITERIA 


COMPUTER  PROGRAM  USER'S  MANUAL 

I.  INTRODUCTION 


The  FORTRAN  computer  program  described  herein  was  developed 
for  the  Armed  Services  Explosive  Safety  Board  in  partial  ful¬ 
fillment  of  the  requirements  specified  under  contract  DAHC04- 
69-C-0095.  This  study  [1]  was  directed  toward  developing 
additional  data  for  use  in  evaluating  quantity-distance 
criteria  for  blast  damage.  Updated  analytical  and  experi¬ 
mental  evidence  has  provided  a  basis  for  storage  distance 
requirements  applicable  to  a  number  of  specific  civilian 
targets*  not  previously  considered.  It  will  be  convenient  to 
have  a  copy  of  Reference  1  available  to  use  in  connection  with 
this  manual. 

Computer  programmed  analytical  methodology  has  been  formulated 
which  will  develop  a  minimum  separation  distance,  for  a 
selected  target  and  quantity  of  explosive,  which  minimizes 
the  probable  risk  that  the  blast  damage  resulting  from  an 
accidental  detonation  will  exceed  a  predetermined  acceptable 
level.  These  procedures  were  derived,  where  possible,  through 
application  of  results  from  more  recent  analytical  and  experi¬ 
mental  blast  damage  programs  in  which  the  structures  or  struc¬ 
tural  elements  subjected  to  blast  could  be  related  through 
similar  structural  properties  and  construction  techniques  to 
the  civilian  targets  of  interest.  This  computer  program  was 
used  to  derive  the  results  presented  in  Reference  1.  In  that 
reference,  plots  showing  the  computed  separation  distances 
are  presented  for  a  number  of  target  parameter  variations , 
and  the  defined  damage  level  in  each  target  when  exposed  to 
blast  from  five  charge  sizes.  These  data  may  be  compared  with 
present  quantity  distance  standards,  which  are  shown  on  each 


*The  term  "target"  is  used  throughout  this  report  to  refer  to 
personnel  and  to  public  structures  which  normally  require  pro¬ 
tection  from  accidental  detonation  of  stored  explosives.  Such 
targets  include  civilian  personnel,  various  types  of  private 
dwellings  and  buildings,  mobile  houses,  commercial  aircraft, 
automotive  vehicles,  other  storage  magazines,  etc. 


1 


plot  for  the  convenience  of  the  reader.  Plots  of  the  incident 
pressure  and  impulse  values  associated  with  these  distances 
have  also  been  prepared  for  each  target  type.  These  latter 
plots  are  clearly  of  the  hyperbolic  form  which  has  been  pre¬ 
dicted  by  earlier  workers  in  the  field  [2]. 
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II.  SELECTED  TARGETS  AND  DAMAGE  LEVEL  SPECIFICATIONS 


Ten  representative  or  generic  types  of  civilian  targets  and 
variations  thereof  may  be  analyzed  with  the  existing  computer 
model.  These  are: 


1.  Modern  split-level  house  constructed  of  a  combina¬ 
tion  of  masonry  and  frame. 

2.  A  modern  church  with  laminated  beams  and  an  A-frame 
roof. 

3.  A  modern  elementary  school. 

4.  A  multi-story  office  building. 

5.  A  passenger  bus  moving  around  a  curve  on  a  highway. 

6.  A  camper-pickup  unit. 

7.  A  stationary  mobile  home. 

8.  Personnel. 

9.  The  front  wall  of  a  standard  storage  igloo. 

10.  A  commercial  jet  aircraft. 


Specific  examples  of  these  targets  have  been  described  in 
precise  detail  [1]  so  that  the  functional  relationship  between 
their  actual  component  weight,  elasticity,  strength,  dimen¬ 
sions,  etc.,  and  the  dynamic  response  to  blast  might  be  deter¬ 
mined.  Modifications  of  these  targets  may  be  treated  by 
changing  the  suggested  input  data  parameter  variations  (such 
as  Young's  modulus,  modulus  of  rupture,  etc.,  for  specific 
target  components)  or  the  fixed  parameter  specifications 
within  the  program  itself  (such  as  the  angle  of  inclination 
of  the  house  roof,  etc.) . 


Because  separation  distance  is  usually  an  expensive  require¬ 
ment,  the  acceptance  of  some  minor  or  incipient  structural 
damage  to  the  targets  is  necessary  in  order  to  minimize  these 
costs.  Thus,  specific  degrees  of  acceptable  damage  have  been 
sought  for  each  target  which  minimize  the  risk  of  serious 
injury  to  occupants  without  requiring  separation  distances 
that  are  excessive. 


In  establishing  "acceptable  damage"  criteria  for  each  of  the 
above  targets ,  a  number  of  unacceptable  degrees  of  damage 
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were  first  considered  for  each  structure.  The  explosive 
storage  separation  distance  for  a  given  explosive  quantity 
was  then  established  with  relationship  to  a  "worst  case"  or 
first  damage  level.  For  the  stationary  targets  (home,  church, 
school,  office  building,  and  the  mobile  home) ,  two  levels  of 
unacceptable  damage  or  risk  levels  were  considered.  The  first 
was  damage  which  constituted  a  critical  hazard  to  occupants 
from  high  velocity  glass  fragments  from  broken  windows.  The 
second  was  blast  damage  sufficient  to  cause  an  immediate 
hazard  to  occupants  of  the  building  as  a  result  of  collapse 
of  the  roof  or  walls  of  the  structure.  In  the  analysis  of 
these  structures,  levels  of  pressure-impulse — for  each  selec¬ 
ted  explosive  quantity — which  cause  each  of  these  types  of 
damage  were  established.  The  greatest  separation  distance  at 
which  unacceptable  damage  occurs  denotes  the  critical  damage, 
and,  thus,  is  designated  as  the  "worst  case". 

Acceptable  damage  to  each  structure  is  identified  as  all  struc¬ 
tural  damage  of  a  degree  up  to  but  not  including  the  specified 
critical  damage.  For  example,  if  a  house  roof  is  predicted  to 
fail  before  the  brick  or  frame  walls--and  also  before  flying 
glass  becomes  a  significant  hazard--then  the  acceptable  level 
would  include  all  damage  up  to,  but  not  including,  roof  collapse. 
In  this  case,  a  distance  corresponding  to  the  blast  forces 
required  to  crack  the  roof  rafters,  but  not  sufficient  to  struc¬ 
turally  collapse  the  roof,  is  identified  as  the  recommended 
separation  distance. 

The  defined  acceptable  damage  levels  for  all  targets  are  dis¬ 
cussed  in  the  individual  treatment  of  each  target  in  Section 
IV. 


A  separate  model  of  blast  response  has  been  developed  for  each 
structural  component  of  each  target  which  can  reasonably  be 
expected  to  suffer  damage  within  the  range  of  charge  sizes  and 
scaled  distances,  AD,  considered  ( 4<A D <1 0 0 ) .  Thus,  for 
example,  the  house  is  analyzed  in  terms  of  the  behavior  of 
its  frame  roof,  frame  walls,  and  masonry  walls.  The  first 
failure  of  any  of  the  three  components  (to  the  established 
maximum  acceptable  damage  level)  defines  a  minimum  distance 
requirement  for  the  entire  target. 

The  dynamic  interaction  models  consider  both  the  elastic  and 
plastic  deformation  involved  in  target  response  and  several 
of  these  models  rely  heavily  on  the  concept  of  requiring  a 


4 


critical  impulse  to  be  imparted  to  the  target  member  in  no 
more  than  one  fourth  of  its  fundamental  resonant  period  [5]. 
Where  the  blast  forces  act  against  both  front  and  back  sur¬ 
faces  of  a  target  or  component,  both  forces  have  been  con¬ 
sidered.  The  models  rely  upon  the  work  of  the  Ballistic 
Research  Laboratories  for  descriptions  of  the  blast  wave 
phenomena  in  terms  of  peak  incident  pressure,  incident  impulse, 
and  positive  duration  [3].  These  values,  together  with  appro¬ 
priate  reflection  factors  and  target  dimensions,  have  been 
used  to  define  the  reflected  pressure  and  impulse  to  which  a 
particular  target  element  must  respond  as  a  function  of  time. 

This  computer  program  provides  output  data  relative  to  the 
critical  periods  of  each  component,  associated  critical 
impulse,  required  separation  distance  to  meet  the  maximum 
acceptable  damage  criterion,  and  the  blast  exposure  of  the 
other  target  components.  The  program  will  accept  a  series  of 
structural  properties  as  inputs  to  the  program  so  that  the 
effects  of  variations  in  tensile  strength,  elasticity,  mass, 
etc. ,  can  be  shown  for  each  target  element. 

While  the  models  and  computer  program  elements  are  extensively 
concerned  with  target  structural  properties,  it  will  be  valu¬ 
able  to  remember  that  reasonable  protection  of  occupants, 
rather  than  the  protection  of  structure,  forms  the  basis  for 
this  program. 
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III.  COMPUTER  PROGRAM  ORGANIZATION 


A.  SUBROUTINES 

The  computer  program  has  been  formulated  as  a  series  of  linked 
subroutines  in  order  that  it  might  also  be  used  on  smaller 
storage  capacity  computers.  A  card  input  is  utilized.  Flow 
charts  and  a  listing  are  presented  in  Section  VI.  All  input 
data  pertaining  to  blast  parameters  and  target  component 
parameter  variations  are  read  in  as  an  initial  step  by  the 
HEADR  subroutine  (see  flow  charts  -  blocks  1  through  17) . 
Thereafter,  a  single  control  card  and  the  CALLR  subroutine 
direct  the  program  to  the  desired  target  analysis  subroutines. 
Each  target  is  evaluated  separately  with  the  exception  of  the 
office  building.  In  its  treatment,  separate  subroutines  are 
used  to  analyze  the  half  wall  with  windows  and  the  solid  arch¬ 
ing  type  block  wall  with  no  openings. 

The  appropriate  subroutines  and  applicable  blocks  of  the  flow 
charts  for  each  target  analysis  are  shown  in  Table  I. 

As  linked  subroutines  are  used,  separate  sets  of  monitor  con¬ 
trol  cards,  Common  and  Dimension  Statements,  and  output  format 
cards  are  provided  in  each  subroutine  (see  source  deck  listing. 
Section  VI).  Only  one  set  is  necessary,  of  course,  if  the 
program  is  modified  to  one  continuous  program.  Some  duplica¬ 
tion  in  the  numbering  of  output  format  statements  has  also 
been  used;  this  must  be  avoided  if  modification  to  one  con¬ 
tinuous  program  is  made. 


B.  INPUT  DATA 

In  order  to  facilitate  the  discussions  of  the  analysis  of  each 
target,  the  required  computer  input  data  cards  are  outlined 
in  Table  II. 
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TABLE  I 


TABLE  I. 

Subroutine  Identification  and 
Flow  Chart  Sections 

Applicable 

Subroutine 

Number 

Subroutine  Name 

Flow  Chart 
Block  Numbers 

1 

Input  Data  -  HEADR 

1-17 

2 

Control  Program  -  CALLR 

18-24 

3 

HOUSE  ANALYSIS  -  HOUSE 

25-120 

4 

CHURCH  ANALYSIS  -  CHRCH 

121-182 

5 

SCHOOL  ANALYSIS  -  SCHOL 

183-253 

6 

OFFICE  HALF  WALL  ANALYSIS  - 

OHAFW 

254-296 

7 

OFFICE  FULL  WALL  ANALYSIS  - 

OFULL 

297-337 

8 

BUS  ANALYSIS  -  BUS 

338-381 

9 

CAMPER-PICKUP  ANALYSIS  - 

CAMPR 

382-425 

10 

MOBILE  HOME  ANALYSIS  - 

TRLER 

426-466 

11 

STORAGE  IGLOO  ANALYSIS  - 

DOOR 

467-498 

12 

PERSONNEL  ANALYSIS  -  MAN 

499-529 

13 

AIRCRAFT  ANALYSIS  -  AIRC 

530-558 

7 


TABLE  II.  Computed  Input  Data 

Card  No. Description  and  Format 

1  Charge  Sizes  -  (FORMAT  5F12.2)  :  Five  charge 

sizes  must  be  specified.  No  limitation  on  charge 
size  is  made;  however,  only  scaled  distances 
4<Ad<100  are  considered.  Therefore,  only  dis¬ 
tances  applicable  to  this  range  are  considered. 
In  Reference  1  charge  sizes  of  1,000,  10,000, 
100,000,  1  million,  and  9  million  pounds  were 
used.  (Array  Name  -  CS(I),I=1,5) 

2-11  Peak  Incident  Overpressures  -  (FORMAT  5F12.2) : 

Ten  cards  are  used  to  input  the  incident  over¬ 
pressures  corresponding  to  scaled  distances 
4<Ad<. 40;  (1)  and  45<.Xj}_<100  ;  (5)  *  .  These  values 
may  be  determined  from  Reference  3  or  10. 

(Array  Name  -  PO(I) ,1=1,50) 

12-21  Scaled  Positive  Durations  -  (FORMAT  5F12.2): 

Ten  cards  are  used  to  input  scaled  positive 
durations  corresponding  to  the  scaled  distances 
4<Ad<.40;(1)  and  45<Xj-><.100 ;  (5)  .  These  scaled 
durations  are  of  course  multiplied  by 
3, - 

/  charge  size  to  obtain  the  actual  positive 
duration  for  use  in  the  program.  Scaled  dura¬ 
tions  in  the  above  range  may  be  obtained  from 
Reference  1  or  3.  (Array  Name  -  TAU (I) , 1=1 ,50) 

22-31  Reflected  Peak  Overpressures  for  Non-Normal 

Surfaces  -  (FORMAT  512.2)  :  Ten  cards  are  used 
to  indicate  the  reflected  peak  overpressures  for 
surfaces  with  approximately  a  70-degree  slope 
(as  measured  from  the  vertical) .  These  values 
are  used  in  the  house  roof  and  aircraft  fuselage 
analyses.  They  again  are  given  for  scaled  dis¬ 
tances  in  the  ranges  4<Ad£40  ;  (1)  and  45<Xd£ 

100;  (5).  These  values  may  be  obtained  from 
Reference  1.  (Array  Name  -  CPR(I) ,1=1 , 50) 


*This  notation  should  be  interpreted  as  the  consideration  of 
all  integer  values  of  \D  in  the  range  of  4£Xq£40  and  all  inte¬ 
ger  values  in  increments  of  5  in  the  range  45<AD<.100. 
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TABLE  II  (Continued) 


Card  No. 


32 


33 


34 


35 


36 


37 


Description  and  Format 

Young's  Modulus  of  the  House  Rafter  -  (FORMAT 
412.2):  Up  to  4  values  of  E  may  be  specified. 

For  Douglas  fir  lumber,  values  in  the  range  of 
1.1  to  1.5  million  psi  are  generally  applicable 
[4],  The  number  of  variables  specified  here 
must  correspond  to  the  number  used  on  cards  33 
and  34.  (Array  Name  -  HRE (I) ,1=1,4) 

Weight  of  16-Inch  Roof  Section  Per  Lineal  Foot 

of  Rafter  -  (FORMAT  4F12.2)  :  Four  values  of  the 
estimated  weight  (per  foot  of  rafter)  of  each 
16-inch  wide  section  of  roof  may  be  specified. 

For  the  composition  roof  considered  in  Reference 
1,  values  from  10.4  to  11.3  lbs/ft  were  used. 
(Array  Name  -  HRW (I)  ,1=1 , 4) 

Modulus  of  Rupture  of  House  Rafter  -  (FORMAT 
4F12.2):  Up  to  4  parameter  variations  of  the 

modulus  of  rupture  may  be  specified.  Reference 
4  indicates  that  for  Douglas  fir,  the  range  of 
8,500  to  12,000  psi  is  generally  applicable. 
(Array  Name  -  HRS (I ) , 1=1 , 4) 

Young 1 s  Modulus  of  House  Frame  Member  -  (FORMAT 
4F21.2) :  Four  estimates  of  the  Young's  modulus 

for  the  2-  by  4-inch  stud  may  be  indicated. 

These  values  correspond  to  those  of  the  roof 
rafter.  The  number  indicated  must  agree  with  the 
number  used  on  cards  36  and  37.  (Array  Name  - 
HFE (I) ,1=1,4) . 

Weight  of  16-inch  Frame  Section  per  Lineal  Foot 

of  Wall  -  (FORMAT  4F12.2) :  Up  to  4  estimates 
may  be  made.  Values  in  the  range  of  18.3  to 
20.0  lbs/ft  are  applicable.  (Array  Name  - 
HFW(I) ,1=1,4) 

Modulus  of  Rupture  of  House  Frame  Member  (FORMAT 
4F12.2) :  Up  to  4  values  are  used.  Values  used 

for  the  house  rafter  are  appropriate.  (Array 
Name  -  HFS (I) ,1=1,4) 
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TABLE  II  (Continued) 


Card  No. 


38 


39 


40 


41 


42 


43 


Description  and  Format 

Young's  Modulus  for  House  Brick  Wall  -  (FORMAT 
4F12.2) :  Up  to  4  estimates  may  be  specified.  A 

range  of  1  to  1.3  million  psi  is  considered  gen¬ 
erally  applicable  for  8-inch  thick  brick  walls. 

The  number  of  values  used  here  must  correspond 
to  the  number  used  on  card  39.  (Array  Name  - 
HBE(I) ,1=1,4) 

Tensile  Strength  of  the  Brick  Mortar  -  (FORMAT 
4F12.2):  Up  to  4  values  are  specified.  Values 

from  30  to  70  psi  are  appropriate.  (Array  Name  - 
HBT(I)  ,1=1,4) 

Young's  Modulus  of  the  Church  Roof  Beam  -  (FORMAT 
4F12.2) :  Up  to  4  values  may  be  specified  for  the 

laminated  Douglas  fir  beams.  The  same  range  of 
values  as  indicated  for  the  house  rafter  is  con¬ 
sidered  to  be  applicable.  The  number  of  values 
used  here  must  correspond  to  the  number  used  on 
cards  41  and  42.  (Array  Name  -  CRE(I) ,1=1,4) 

Weight  of  the  15-Foot  Roof  Section  Per  Lineal 

Foot  of  Church  Beam  -  (FORMAT  4F12.2) :  Up  to  4 
parameter  variations  of  the  weight  of  each  15- 
foot  wide  section  of  church  roof  (including 
beam)  are  specified.  Weights  of  200  to  230 
lbs/ft  are  appropriate  for  the  church  considered. 
(Array  Name  -  CRW(I) ,1=1,4) 

Modulus  of  Rupture  of  the  Church  Roof  Beam  - 
(FORMAT  4F12.2) :  Four  estimates  may  be  specified. 
The  same  range  of  values  used  for  card  34  is 
appropriate.  (Array  Name  -  CRS (I) ,1=1 , 4) 

Young's  Modulus  of  the  Church  Roof  Decking  - 

(FORMAT  4F12.2):  Up  to  4  values  are  specified; 
the  total  number  used  must  correspond  with  the 
number  used  on  cards  44  and  45.  For  the  church 
considered,  values  in  the  range  of  1.3  to  1.55 
million  psi  are  applicable.  (Array  Name  - 
CDE(I)  ,1=1,4) 
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TABLE  II  (Continued) 


Card  No. 


44 


45 


46 


47 


48 


49 


50 


Description  and  Format 

Weight  of  the  Decking  Member  and  Roofing  Per 

Lineal  Foot  of  Decking  -  (FORMAT  4F12.2) :  Up  to 
4  values  may  be  specified.  Values  in  the  range 
of  13  to  14.5  lbs/ft  are  appropriate.  (Array 
Name  -  CDW(I) ,1=1,4) 

Modulus  of  Rupture  of  the  Church  Decking  Member  - 

(FORMAT  4F12.2):  Up  to  4  estimates  are  speci¬ 
fied.  The  range  of  1.1  to  1.4  million  psi  is 
applicable.  (Array  Name  -  CDS ( I ) , 1=1 , 4 ) 

Young's  Modulus  of  the  School  Stone  Wall  -  (FOR¬ 
MAT  4F12.2):  Up  to  4  values  are  specified;  the 
number  treated,  however,  must  agree  to  the  number 
of  variables  used  on  card  47.  Similar  values  as 
on  card  38  are  appropriate.  (Array  Name  - 
SBWE(I) ,1=1,4) 

Tensile  Strength  of  the  School  Wall  Mortar  - 
(FORMAT  4F12.2):  Four  parameter  variations  are 
used.  Similar  values  to  card  39  are  applicable. 
(Array  Name  -  SBT ( I ) , 1=1 , 4) 

Young's  Modulus  of  the  School  Roof  Beams  -  (FOR¬ 
MAT  4F12.2) :  Four  estimates  may  be  denoted-- the 
member  used  must  correspond  to  the  number  used 
on  cards  49  and  50.  Values  from  1.3  to  1.55 
million  psi  are  appropriate  for  this  school  con¬ 
struction.  (Array  Name  -  SRE ( I ) , 1=1 , 4 ) 

Weight  of  2-Foot  Wide  Section  of  School  Roof  Per 

Lineal  Foot  of  Beam  -  (FORMAT  4F12.2) :  Up  to  4 
values  are  specified.  For  this  asphalt  construc¬ 
tion,  the  weight  ranges  between  30  and  40 
lbs/ft.  (Array  Name  -  SRW (I) , 1=1 , 4) 

Modulus  of  Rupture  of  School  Roof  Beam  -  (FORMAT 
4F12.2):  Up  to  4  estimates  are  specified. 

Values  in  the  range  of  1.1  to  1.3  million  psi 
are  thought  applicable.  (Array  Name  -  SRS(I), 
1=1,4) 
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TABLE  II  (Continued) 


Card  No. 


51 


52 


53 


54 


55 


56 


57 


Description  and  Format 

Young's  Modulus  of  the  Office  Hall  Wall  -  (FOR¬ 
MAT  4F12.2):  Up  to  4  values  are  indicated,  the 
number  of  which  must  agree  with  the  number  used 
in  card  52.  Values  used  on  card  38  are  appro¬ 
priate.  (Array  Name  -  OHE  (I) ,1=1 ,4) 

Tensile  Strength  of  the  Half  Wall  Mortar  -  (FOR¬ 
MAT  4F12.2):  Up  to  4  values  are  specified. 
Similar  values  to  card  39  may  be  used.  (Array 
Name  -  OHT(I) ,1=1,4) 

Young's  Modulus  of  the  Office  Full  Wall  -  (FOR¬ 
MAT  4F12.2) :  Up  to  4  parameter  variations  are 
denoted.  The  number  indicated  must  agree  to 
that  used  on  card  54.  Values  used  for  the  half 
wall  may  be  used.  (Array  Name  -  OFE (I) ,1=1 , 4) 

Tensile  Strength  of  the  Office  Full  Wall  Mortar  - 

(FORMAT  4F12.2) :  Up  to  4  values  are  used. 
Appropriate  values  lie  in  the  range  of  30  to  70 
psi.  (Array  Name  -  OFT  (I)  ,  1=1 , 4) 

Bus  Center  of  Gravity  Location  -  (FORMAT  4F12.2) : 
Up  to  4  estimates  are  used;  the  number  used  on 
this  and  card  number  56  must  agree.  These  values 
indicate  the  height  above  the  ground  plane  of  the 
c.g.  and  probably  lie  in  the  range  of  38  to  42 
inches.  (Array  Name  -  BCG (I) ,1=1 ,4) 

Bus  Weight  -  (FORMAT  4F12.2):  Up  to  4  weights 
of  the  bus  may  be  indicated.  For  this  bus, 

28,000  to  32,000  lbs  is  a  reasonable  approxima¬ 
tion.  (Array  Name  -  BW(I),I=1,4) 

Camper  Center  of  Gravity  Location  -  (FORMAT 
4F12.2):  Four  estimates  of  the  height  of  the 

c.g.  above  the  ground  plane  may  be  indicated. 

The  number  used  must  agree  with  the  number  of 
parameter  values  used  on  cards  58,  59,  and  60. 
These  values  lie  in  a  range  of  52  to  58  inches. 
(Array  Name  -  CCG(I) ,1=1,4) 
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TABLE  II  (Continued ) 


Card  No . 

Description  and  Format 

58 

Pickup  Truck  Center  of  Gravity  Location  -  (FOR- 
MAT  4F12.2):  Four  c.g.  locations  are  indicated. 

Appropriate  values  are  contained  in  the  range  of 
30  to  34  inches  above  the  ground  plane.  (Array 

name  -  CTG(I) ,1=1,4) 

59 

Truck  Weight  -  (FORMAT  4F12.2) :  Up  to  4  truck 
weights  are  specified.  A  range  of  3,600  to 

4,800  pounds  is  considered  appropriate.  (Array 

Name  -  TRW  (I)  ,1=1 ,4) 

60 

Camper  Weight  -  (FORMAT  4F12.2):  Four  camper 
weights  may  be  used.  A  range  of  2,000  to  2,200 
pounds  is  suggested.  (Array  Name  -  CPW ( I )  , 1=1 , 4 ) 

61 

Mobile  Home  Center  of  Gravity  Location  -  (FORMAT 
4F12.2):  Up  to  4  c.g.  locations  are  considered. 

Appropriate  heights  above  the  ground  plane  lie 
in  the  range  of  60  to  72  inches.  The  number 
used  must  agree  with  that  used  in  card  62. 

(Array  Name  -  TRGG(I) ,1=1,4) 

62 

Mobile  Home  Weight  -  (FORMAT  4F12.2) :  Four  esti- 
mates  may  be  specified.  A  range  of  12,000  to 
14,000  pounds  was  considered  in  Reference  1. 
(Array  Name  -  TRWT(I) ,1=1,4) 

63 

Modulus  of  Rupture  of  the  Storage  Igloo  Steel 

Door  -  (FORMAT  4F12.2) :  Up  to  4  values  may  be 
used.  Appropriate  values  are  considered  to  lie 
in  the  range  of  50  to  55  thousand  psi.  (Array 

Name  -  DORS (I) ,1=1,4) 

64 

Weight  of  Man  -  (FORMAT  4F12.2) :  Four  weights 
of  the  personnel  target  may  be  entered.  The 
total  number  of  values  should  correspond  to  that 
used  for  card  65.  (Array  Name  -  WMAN(I) ,1=1,4) 

65 

Presented  Area  of  Man  -  (FORMAT  4F12.2) :  Up  to 

4  presented  areas  may  be  used.  The  standing 
"average"  man  presents  an  area  of  5.5  to  8.0  sq . 
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TABLE  II  (Continued) 


Card  No. Description  and  Format 

ft.  depending  upon  his  orientation.  (Array 
Name  -  PMAN(I) ,1=1,4) 

66  Young's  Modulus  of  Aircraft  Frame  Member  -  (FORMAT 

4F12.2):  Four  values  may  be  used.  For  the  frame 

member  considered  values  from  10.6  to  11.6  mil¬ 
lion  psi  are  used.  (Array  Name  -  AIRE(I) ,1=1,4) 

67  Dynamic  Overpressure  Factors  for  Aircraft  -  (FOR¬ 
MAT  5F12.2) :  Two  cards  are  used  to  input  these 
values.  The  ten  values  used  must  be  .25,  1.25, 
1.42,  1.54,  1.68,  1.72,  1.75,  1.78,  1.81.  These 
numbers  are  used  to  indicate  the  relationship  of 
static  loading  to  dynamic  loading  in  the  aircraft 
analysis.  (Array  Name  -  AIRF ( I ) , 1=1 , 10 ) 

68  Program  Control  Card  -  (FORMAT  1212) :  The  first 
two-digit,  fixed-point  entry  on  the  card  is 
denoted  as  NOT,  the  total  number  of  target  sub¬ 
routines  to  process;  the  following  11  entries 
(IA(I) 1=1, 11)  indicate  the  order  of  targets  to 
be  processed. 


The  control  card  is  used  to  direct  the  flow  of  the  program. 

The  first  entry,  NOT,  on  this  card  indicates  the  total  number 
of  target  subroutines  to  process.  The  order  in  which  they  are 
to  be  analyzed  is  denoted  in  the  preceeding  two-digit,  fixed- 
point  fields  (denoted  as  IA(I) 1=1,11) .  The  order  of  targets 
is  shown  in  block  24  of  the  flow  charts.  To  illustrate, 
assume  the  three  targets:  the  house,  the  school,  and  the  man 
translation  are  desired  to  be  analyzed.  To  achieve  this  pro¬ 
gramming,  NOT  is  designated  as  3,  IA(1)  is  entered  as  1,  IA(2) 
is  set  equal  to  3  and  IA(3)  is  denoted  as  9.  No  entry  is  made 
for  IA(4)  through  IA(ll). 

It  also  must  be  noted  that  the  computer  program  multiplies  all 
Young's  modulus  (E)  values  by  10^  and  all  modulus  of  rupture 
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(a)  values  by  10^  whenever  used.  This  must  be  accounted  for 
in  the  entry  of  the  input  data.  To  reflect  the  desired  E 
values  of  1  and  1.3  million  psi  for  the  house  rafter,  for 
example,  it  is  necessary  to  enter  1.00  and  1.30  on  card  32. 
Analogous  adjustments  must  be  made  for  all  modulus  of  rupture 
values.  All  other  values  are  entered  with  no  modifications. 
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IV.  ANALYSIS 


A.  SPLIT-LEVEL  HOUSE 

The  selected  representative  type  of  residential  home  construc¬ 
tion  is  illustrated  in  Figure  1.  It  is  described  in  detail 
in  Reference  1.  In  brief  summary,  the  house  roof  is  composed 
of  2-  by  8-inch  Douglas  fir  rafters  on  16-inch  centers.  The 
rafters  are  17  feet  long  and  have  a  5/12  slope.  Each  16-inch 
section  of  composition  roof  (including  rafters)  supported  by 
a  rafter  weighs  approximately  10.4  to  13.3  pounds  per  lineal 
foot  of  rafter.  The  7.5-foot  high  brick  walls  are  of  8-inch 
solid  brick  construction;  there  is  no  reinforcing.  The  frame 
wall  is  7.5  feet  high  containing  a  door  and  2-foot  8-inch  by 
4-foot  6-inch  windows.  Two-  by  4-inch  studs  on  16-inch  cen¬ 
ters  support  the  7.5-foot  high  frame  wall.  Each  16-inch  sec¬ 
tion  of  f iberboard-plas terboard  wall  weighs  approximately 
18.3  to  20.0  pounds  per  lineal  foot  of  stud. 

Three  exterior  surfaces  are  considered  as  possible  vulnerable 
components  in  the  blast  damage  assessment  model  for  this 
house:  the  roof,  a  portion  of  the  brick  wall  containing  no 

openings ,  and  the  frame  section  containing  window  and  door 
openings.  For  all  of  the  components,  the  risk  level,  or 
unacceptable  damage,  is  defined  as  structural  collapse. 
Therefore,  damage  levels  which  structurally  degrade  each  com¬ 
ponent  but  do  not  cause  their  structural  collapse  have  been 
defined.  For  the  roof  and  frame  wall  sections,  this  damage 
is  regarded  as  the  cracking,  but  not  complete  severence,  of 
one  or  more  of  the  main  2-  by  8-inch  rafters  or  2-  by  4-inch 
wall  support  members.  For  the  brick  wall,  a  deflection  of 
approximately  5  inches  at  its  center  is  considered  sufficient 
to  seriously  damage  the  wall  through  extensive  cracking,  but 
not  to  cause  partial  or  complete  collapse. 

In  the  computer-programmed  methodology,  the  minimum  distance 
for  each  specific  charge  size  which  relates  to  the  first-to- 
occur  of  the  above  three  criteria,  is  determined.  Then  the 
required  impulse  and  calculated  impulse  achieved  at  this 
specified  distance  is  determined  for  all  parameter  variations 
considered  for  each  of  the  other  structural  components,  in 
order  to  insure  that  this  truly  is  the  minimum  distance. 
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The  analysis  of  the  wood  portions  of  the  house  (roof  and  frame 
wall)  and  other  targets  is  based  on  Sewell's  concept  of  achiev¬ 
ing  a  "critical  impulse  within  a  critical  time  period"  [5] . 

His  premise  is  that  if  a  blast  wave  impinging  upon  a  structural 
surface  creates  a  velocity  in  the  elements  of  the  structure, 
relative  to  fixed  points  on  the  surface,  greater  than  a  "criti¬ 
cal  working  velocity",  then  permanent  deformation  of  the  struc¬ 
ture  will  result.  A  further  assumption  is  that  this  critical 
velocity  must  occur  within  one-quarter  period  of  the  fundamental 
resonant  frequency  of  the  structure.  Since  critical  working 
velocity  is  not  known  for  many  materials ,  Sewell  advocates  use 
of  the  critical  impact  velocity;  i.e.,  that  velocity  causing  ten 
sile  failure  in  the  material.  Because  this  criterion  requires 
a  minimum  impulse  in  a  fixed  time  period,  it  accounts  for  both 
the  contribution  of  reflected  pressure  and  impulse  in  that  the 
delivery  time  limitation  requires  the  reflected  overpressure  to 
be  relatively  high  in  order  to  produce  the  necessary  impulse 
within  the  time  limit.  Additionally,  the  characteristics  of 
the  target  are  given  consideration  in  that  the  critical  time  is 
a  function  of  the  resonant  frequency  of  the  target  and  the  criti 
cal  impulse  is  expressed  as  a  function  of  target  characteristics 
These  two  considerations  are  considered  necessary  in  accurately 
modeling  target  response-blast  force  interactions. 

In  the  treatment  of  wood  structures,  the  acceptable  level  of 
damage  is  based  on  simply  cracking,  but  not  completely  break¬ 
ing,  a  particular  wood  member.  Therefore,  the  elastic  phase 
of  deflection,  up  to  and  including  the  yield  point,  is  of  pri¬ 
mary  interest.  In  the  application  of  Sewell's  criterion  to 
wood  structures,  it  is  necessary  that  the  critical  impulse  be 
imparted  to  the  wood  structure  so  that  it  will  deflect  beyond 
its  yield  point  within  one-quarter  of  its  natural  period. 

1.  Analysis  of  the  Roof  Rafter 

A  typical  17-foot  long  2-  by  8-inch  roof  rafter  and  the  16- 
inch  wide  section  of  roof  it  supports  are  analyzed  as  a  simply 
supported,  uniformly  loaded  beam.  Its  natural  frequency  is 


where : 


radians /second 


(1) 


E  =  Young's  modulus  (psi) 

I  =  the  area  moment  of  inertia  (in^) 
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m  =  the  mass  per  inch  of  beam  (lbs) 


l  =  the  length  of  the  rafter  (in) . 
Its  critical  period,  Cp ,  becomes 


C 


P 


1000  2n 

— - -  •  —  msec. 

4  p 


(2) 


The  critical  impulse,  Ic,  which  must  be  delivered  to  the  beam 
is  that  impulse  which  deflects  the  uniformly  loaded  beam 
beyond  its  point  of  maximum  elastic  deflection  within  Cp  msec. 

The  calculation  of  the  deflection  of  a  loaded  beam  is  based 
on  the  deflection  of  a  loaded  spring.  The  energy,  U,  stored 
in  a  deflected  linear  spring  is 


U  =  P  •  |  (3) 

where  P  is  the  force  required  to  cause  the  deflection  A.  In 
order  to  use  equation  (3) ,  an  expression  for  A  must  be  obtained. 
If  the  load  function,  f (w) ,  and  sufficient  boundry  conditions 
are  known,  this  can  be  obtained  from  a  successive  integration 
of 

,4 

El  =  f  ( w )  (4) 

dx 

where : 

A  =  deflection,  and 
x  =  distance  along  the  beam. 

In  the  analysis  of  the  simply  supported,  uniformly  loaded 
beam,  the  bending  moment  equation,  M,  is  used  in  (4)  and  A 
obtained  from 
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El 


-M 


(5) 


The  equation  thus  obtained  for  A  will  involve  a  load,  w.  To 
find  the  maximum  permissible  value  for  w,  wQ,  a  maximum  tensile 
fiber  stress,  omax ,  is  obtained  and  substituted  in  the  equation 


w  °maxl 

Mmax  -  n 


(6) 


The  Mmax  thus  obtained  is  substituted  for  M  in  the  bending 
moment  equation  along  with  an  x  that  corresponds  to  the  loca¬ 
tion  along  the  beam  of  maximum  moment.  In  the  case  of  the 
uniformly  loaded  beam 


M  = 


wx^ 

2 


The  maximum  moment  occurs  when  x  =  £/2 ,  so  that 

£  £  w  £2  w£2 

Mmax  ~  2  W  2  ~  2  4  ~  8 

_  8  Mmax  _  8_  cmaxI 


(7) 


The  internal  energy,  U,  is  determined  by  integration  of  equa¬ 
tion  (4)  over  the  length  of  the  beam,  where 


p  =  wQ  dx 


(8) 


In  this  case,  the  following  simple  expression  is  obtained 


U 


240EI 


(9) 


The  necessary  unit  impulse,  Ic,  to  provide  this  required 
energy  is  then  computed  from 


I  = 


1000- (2Um) 


1/2 


(area  of  the  beam  (in  ) ) 


,  (psi-msec) 


(10) 
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The  calculated  critical  period  and  critical  impulse  for  the 
roof  rafter  depend  upon  the  assumed  strength  and  elastic 
properties  of  Douglas  fir.  The  parameters  E,  w,  and  a  are 
known  only  within  a  range  of  values;  thus,  selected  parameter 
values  for  these  variables  are  entered  as  input  data  on  cards 
32,  33,  and  34.  A  minimum  separation  distance  for  each  selec¬ 
ted  charge  size  is  then  determined  for  each  resultant  set  of 
calculated  Cp  and  IQ. 

Blast  loading  models  presented  in  Reference  6  have  been  used 
in  this  computer  program.  Individual  loading  properties  and 
the  resultant  imparted  impulse  on  different  surfaces  of  a 
target  may  be  calculated,  as  well  as  net  loading  and  impulse 
values.  A  detailed  discussion  of  these  models  and  the  treat¬ 
ment  of  the  non-horizontal  roof  of  the  house  is  given  in 
Reference  1.  The  appropriate  blast  loading  model  is  used  in 
the  computer  program  to  determine  a  minimum  separation  distance 
for  the  house  roof  for  each  charge  size  and  combination  of 
selected  parameter  values.  At  each  of  these  specified  dis¬ 
tances,  similar  estimates  are  made  of  the  reflected  impulse 
imparted  to  the  house  brick  wall  and  frame  portion.  Compara¬ 
tive  values  are  indicated  following  a  brief  discussion  of  the 
frame  section  and  brick  wall. 

2.  Analysis  of  the  Frame  Wall  Section 

The  2-  by  4-inch  wall  support  member  is  analyzed  as  a  uniformly 
loaded,  simply  supported  beam.  Its  dead  weight  is  determined 
from  its  own  weight  and  the  16-inch  wide  section  of  wall  it 
supports.  Calculation  of  its  critical  periods  and  impulses 
are  made  using  equations  (1)  through  (10) . 

Selected  values  of  E,  w,  and  o  are  entered  on  input  data  cards 
35,  36,  and  37.  A  net  blast  loading  concept  is  employed  in 
which  the  net  impulse  is  determined  through  integration  of  the 
overpressure  function  from  time  t  =  0  to  a  time  at  which  the 
overpressure  function  decreases  to  0.2  psi. 

It  has  been  concluded  from  the  results  of  several  large-scale 
tests,  that  a  steady  pressure  of  approximately  0.5  psi  over  a 
time  period  equal  to  the  wall  calculated  critical  period  is 
required  for  a  simply  supported  8-inch  thick  brick  wall  in 
order  to  cause  tensile  bond  failure  in  the  mortar. 
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In  the  computer  model  of  the  brick  wall  behavior,  the  calcu¬ 
lation  of  the  imparted  reflected  impulse  is  based  on  these 
assumptions.  The  blast  loading  model  which  reflects  no  open¬ 
ing  in  a  perpendicular  wall  is  used;  it  is  thoroughly  discussed 
in  Reference  1.  The  required  unit  impulse  is  computed  as  the 
difference  in  the  reflected  impulse  imparted  to  the  outside  of 
the  frame  wall  within  its  one-quarter  period,  and  that  imparted 
to  the  rear  surface.  This  concept  is  also  detailed  in  Refer¬ 
ence  1. 


3.  Analysis  of  the  House  Brick  Wall 

The  derivation  of  the  resistance  functions,  the  work  required 
to  severely  crack  the  wall,  and  the  computed  critical  impulse 
necessary  to  do  that  work  are  discussed  in  Reference  1.  From 
this  study,  it  is  concluded  that  when  the  overpressure  decays 
to  a  low  level,  the  latter  part  of  the  decay  process  is  of 
little  consequence  to  the  effective  impulse  imparted  to  the 
wall.  This  point  of  minimum  effective  pressure  is  concluded 
to  be  in  the  range  of  0.1  to  0.3  psi.  In  this  computer  model, 
the  calculation  of  reflected  impulse  has  been  based  on  the 
integration  from  time  t  =  0  to  a  time  at  which  the  overpressure 
decays  to  0.2  psi. 

The  calculation  of  reflected  impulse  required  to  severely 
crack  the  wall  also  depends  upon  the  blast  loading  forces 
imparted  to  the  roof  which,  of  course,  change  with  variations 
in  separation  distance.  Thus,  in  the  program  output,  a 
required  unit  impulse,  based  on  the  blast  loading  on  the  roof 
for  each  given  charge  and  separation  distance,  is  indicated 
for  each  set  of  parameter  variations  considered  for  the  brick 
wall.  Up  to  four  values  each  of  Young's  modulus  and  the  ten¬ 
sile  strength  of  mortar  are  entered  on  input  cards  38  and  39, 
respectively . 

4.  Sample  Set  of  Computed  Values 

The  computer  output  for  the  house  analysis  is  arranged  so  that 
the  computed  critical  periods  and  impulses  are  printed  for  the 
house  roof,  the  frame  section,  and  the  brick  wall  section,  in 
sequential  order.  Then,  for  each  set  of  parameter  variations 
for  the  house  roof,  the  five  charge  sizes  are  considered  in 
ordered  to  determine  a  minimum  separation  distance  for  each, 
with  respect  to  the  roof.  As  each  of  these  values  is 
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determined,  the  incident  overpressure,  peak  reflected  pressure 
and  positive  duration  of  the  charge  is  noted,  and  a  probabi¬ 
listic  statement  concerning  the  likelihood  of  glass  fragment 
injury  to  occupants,  is  made.  The  other  two  house  components 
are  then  considered  at  this  distance.  The  required,  or 
critical,  impulse  and  the  reflected  impulse  achieved  at  this 
distance  are  printed  for  every  set  of  values  considered  for 
these  two  components.  Table  III  displays  a  typical  printout 
for  a  set  of  representative  parameters  of  the  house  roof  and 
a  charge  size  of  10,000  pounds. 

Table  III  indicates  that  at  the  minimum  separation  distance  of 
289.9  feet  for  a  10,000-pound  charge,  the  calculated  reflec¬ 
ted  impulses  on  the  frame  portion  and  the  brick  wall  of  the 
house  at  this  distance  are  far  less  than  the  necessary  criti¬ 
cal  impulse  values  computed  for  the  different  sets  of  parameter 
variations  considered. 


B.  CHURCH  BUILDING 

The  A-frame  church  building  is  shown  in  Figure  2.  The  lami¬ 
nated  Douglas  fir  wood  beams  are  7.6  by  16  inches  by  39  feet 
long  spaced  on  15-foot  centers.  The  upper  32  feet  support 
the  decking  and  roofing  while  the  lower  7  feet  support  no 
loading.  The  4-inch  thick  8-inch  wide  double  tongue  and  groove 
decking  is  considered  to  be  simply  supported  on  15-foot  centers. 
The  decking  is  covered  with  asphalt  roofing  and  has  a  combined 
weight  of  from  12  to  13.5  pounds  per  lineal  foot  of  decking. 
Excluding  the  windows,  the  large  (90  feet  long)  roof  is  con¬ 
sidered  the  only  vulnerable  component  of  the  church.  The 
acceptable  damage  levels  are  considered  to  be  the  cracking 
of  the  laminated  beams  or  cracking  of  the  decking  members. 

The  minimum  separation  distance,  therefore,  is  based  on  the 
first  of  these  to  occur. 


1.  Analysis  of  the  Laminated  Beams 

The  39-foot  laminated  beam  and  the  section  of  the  roof  it 
supports  is  analyzed  as  a  step-wise  uniformly  loaded,  simply 
supported  beam  with  a  dead  load  determined  from  its  own 
weight  and  the  weight  of  the  15-foot  section  of  roof  it  sup¬ 
ports  (approximate  weight,  200  to  230  lbs/ft  of  beam) .  Equa¬ 
tions  (1)  through  (10)  are  used  by  the  computer  program  to 
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Figure  2.  Typical  Modern  Church 


determine  the  critical  period  and  impulse.  Up  to  four  esti¬ 
mated  values  each  of  Young's  modulus,  the  weight  supported 
per  foot  of  beam,  and  the  modulus  of  rupture  are  entered  on 
cards  40,  41,  and  42. 

The  church  roof  has  an  angle  of  inclination  of  approximately 
23  degrees  from  the  vertical;  the  treatment  of  this  non¬ 
normal  surface  in  the  calculation  of  the  reflected  impulse  is 
given  in  Reference  1. 

2.  Analysis  of  the  Church  Decking 

The  decking  members  are  also  laminated  Douglas  fir  and  have 
similar  strength  properties  as  the  other  wood  members  con¬ 
sidered.  Up  to  four  values  of  E,  w,  and  a  may  be  entered  on 
input  cards  43,  44,  and  45.  The  critical  period  and  impulse 
calculations  are  made  using  equations  (1)  through  (10).  The 
blast  loading  model  on  a  perpendicular  surface  with  no  back 
pressure  on  the  interior  side  of  the  decking  is  used  to  calcu¬ 
late  the  impulse  imparted  to  the  decking  within  its  critical 
period  for  each  set  of  chosen  parameter  variations.  This 
model  is  discussed  in  detail  in  Reference  1. 

3.  Comparative  Results 

In  the  previous  analysis,  the  separation  distance  at  which 
the  church  roof  first  fails  for  a  10,000-pound  charge  is 
1,184.9  feet  (for  the  set  of  parameters  E  =  1.3  x  10^  psi, 
w  =  215  lbs/ft  of  beam  and  amax  =  11,000  psi).  The  calculated 
impulse  imparted  to  the  decking  at  this  distance  is  approxi¬ 
mately  33  psi-msec.  The  average  critical  impulse  required 
for  the  decking  is  about  110  psi-msec  when  comparative  values 
of  E,  w,  and  a  are  used.  Thus,  for  this  and  all  other  charge 
sizes  considered  in  Reference  1,  the  church  roof  beam  was 
found  to  crack  prior  to  the  decking. 


C.  ELEMENTARY  SCHOOL  BUILDING 

The  roof  of  the  elementary  school  (Figure  3)  is  supported  on 
3-  by  14-inch  by  28-foot  long  fir  beams  on  2-foot  centers. 

The  roofing  is  of  built-up  tar  and  gravel  construction  placed 
on  1-inch  insulation  board.  A  light  tile  ceiling  is  placed 
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Figure  3.  Typical  Elementary  School 


on  the  bottom  of  the  rafters.  The  total  roof  weight  of  a 
2-foot  section  is  approximately  30  to  45  lbs/ft  of  beam. 

The  end  wall  is  8  inches  thick  with  stone  mortared  together 
with  no  reinforcing.  It  is  10  feet  6  inches  high  and  contains 
no  windows.  It  has  the  same  simple  support  conditions  as  the 
house  brick  wall. 

The  windows  in  the  front  wall  vary  somewhat  in  size,  the  maxi¬ 
mum  being  approximately  3  feet  6  inches  by  5  feet.  All  win¬ 
dows  are  1/4-inch  plate  glass  supported  by  aluminum  frames. 

The  school  roof  and  stone  wall  are  the  most  likely  structural 
components  to  collapse  under  a  given  blast  loading.  The  win¬ 
dows,  of  course,  will  fracture  and  present  an  appreciable 
glass  fragment  hazard;  however,  the  loss  of  the  windows  and 
window  frame  portion  of  the  front  wall  will  not  cause  collapse 
of  the  roof  because  of  the  orientation  of  the  roof  support 
beams.  Blast  damage  levels  corresponding  to  cracking  the  roof 
rafters  and/or  severely  cracking  the  stone  wall  therefore  have 
been  defined  as  acceptable  criteria  on  which  to  establish  mini¬ 
mum  separation  distances.  To  reemphasize,  a  deflection  of  the 
center  of  the  8-inch  stone  wall  of  approximately  5  inches  is 
considered  sufficient  to  seriously  crack  the  wall  but  not  to 
cause  partial  or  complete  collapse. 

1.  Analysis  of  the  Stone  Wall 

The  resistance  functions  and  the  unit  impulse  required  to 
severely  crack  the  stone  wall  were  computed  using  the  model 
for  masonry  walls  discussed  in  Reference  1.  As  in  the  treat¬ 
ment  of  brick  and  concrete  walls,  a  steady  pressure  of 
approximately  0.5  psi  was  assumed  to  be  required  to  cause 
tensile  bond  failure  in  the  mortar.  In  the  programmed  ana¬ 
lytical  model,  then,  a  requirement  of  at  least  a  0.5  psi 
average  overpressure  throughout  a  blast  loading  interval 
equal  to  the  calculated  critical  period  was  stipulated.  As 
in  the  analysis  of  the  other  types  of  wall  section,  the 
effective  reflected  impulse  was  estimated  as  the  integral  of 
the  overpressure  function  from  time  t  =  0  to  a  time  at  which 
the  overpressure  level  decreases  to  0.2  psi. 

Equations  (1)  and  (2)  are  used  to  obtain  the  critical  period 
for  a  1-inch  wide  column  of  stone  wall,  analyzed  as  a  uniformly 
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loaded  beam.  From  one  to  four  values  each  of  Young's  modulus 
and  the  tensile  bond  failure  stress  for  mortar  must  be  entered 
on  cards  46  and  47. 

The  lack  of  openings  in  the  stone  wall  prevents  back  pressure 
from  building  up  on  the  inside  of  the  stone  wall;  thus,  the 
blast  loading  on  the  front  face  of  the  wall  is  the  only  con¬ 
sideration.  As  indicated  in  the  treatment  of  other  wall 
sections,  the  required  critical  impulse  for  each  set  of 
parameter  values  depends  to  a  large  degree  upon  the  blast 
loading  on  the  roof.  Under  the  assumption  that  the  stone 
wall  supports  one-half  of  each  roof  section,  the  computer  pro¬ 
gram  calculates  each  unit  impulse  requirement  based  on  one- 
half  of  the  contribution  (per  inch  of  stone  wall)  of  the  cal¬ 
culated  maximum  average  pressure  on  the  roof.  This  procedure 
is  clearly  defined  in  Reference  1. 

2.  Analysis  of  the  School  Roof 

The  school  roof  beams  are  analyzed  as  being  uniformly  loaded. 
Their  critical  periods  and  critical  impulse  calculations  are 
straightforward  using  equations  (1)  through  (10) .  Up  to  four 
parameter  estimations  of  E,  w,  and  o  can  be  entered  on  cards 
48,  49,  and  50.  As  is  the  case  in  all  target  analysis  sub¬ 
routines  in  this  model,  once  a  separation  distance  is  deter¬ 
mined  which  results  in  the  defined  damage  for  one  of  the  vul¬ 
nerable  components  of  that  target,  the  computed  impulse  on 
the  other  component (s) ,  at  that  distance,  and  its  respective 
critical  impulse,  for  all  parameter  variations  considered,  are 
compared.  For  the  parameter  variations  considered  for  the 
school  wall  and  roof  in  Reference  1,  it  was  found  that  the 
school  wall  suffered  the  defined  acceptable  damage  in  all 
cases  prior  to  achieving  the  defined  damage  to  the  school  roof. 


D.  OFFICE  BUILDING  -  HALF  WALL  WITH  WINDOWS 

The  analysis  of  the  concrete  block  walls  (Figure  4)  is  con¬ 
ducted  by  the  computer  program  under  the  same  general  set  of 
assumptions  as  used  in  the  treatment  of  the  other  masonry 
walls.  Selected  values  of  E  and  the  tensile  strength  of  the 
half  wall  are  entered  on  cards  51  and  52.  A  net  loading 
concept  is  used  in  which  the  pressure  buildup  on  the  rear 
surface  of  the  half  wall  is  subtracted  from  the  overpressure 
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function  on  the  front  surface.  Again,  the  reflected  impulse 
values  are  calculated  by  integrating  the  net  overpressure 
function  from  time  t  =  0  to  a  time  at  which  the  overpressure 
decays  to  0.2  psi. 


E.  OFFICE  BUILDING  -  FULL  WALL  WITH  NO  OPENINGS 

The  full  office  wall  is  an  arching  type  wall  under  blast  load¬ 
ing  in  that,  since  it  is  constructed  tightly  between  two 
immovable  concrete  columns,  some  mortar  must  be  crushed  in 
order  to  deflect  the  wall  to  a  point  at  which  serious  cracking 
will  occur.  The  models  for  the  analysis  of  an  arching  wall 
are  taken  from  Reference  7  and  are  described  in  detail  in 
Reference  1.  Based  on  mortar  sample  crushing  tests,  a  steady 
pressure  of  approximately  5  psi  must  be  maintained  for  a 
significant  period  so  that  sufficient  crushing  of  the  mortar 
at  the  top  and  bottom  of  the  wall  may  take  place.  In  the  com¬ 
puter  analysis,  it  is  assumed  that  this  average  pressure  must 
be  maintained  over  a  time  interval  equal  to  its  calculated 
critical  period.  Up  to  four  values  of  E  and  tensile  strength 
are  entered  on  input  cards  53  and  54.  The  reflected  impulse 
for  each  charge  separation  distance  and  set  of  parameter 
values  is  calculated  using  a  front  surface  blast  loading  model 
in  which  the  overpressure  function  is  integrated  to  a  time  at 
which  the  function  decreases  to  0.2  psi. 

Comparison  of  separation  distances  for  the  full  office  wall 
and  half  wall  can  only  be  obtained  by  requesting  both  sub¬ 
routines  with  the  control  card. 


F.  PASSENGER  BUS 

The  prevention  of  overturning  is  the  basis  on  which  acceptable 
damage  criteria  have  been  established  for  the  vehicular  tar¬ 
gets — the  passenger  bus,  camper-pickup  unit,  and  the  mobile 
home.  The  pressure-impulse  requirements  to  overturn  these 
targets  are  therefore  determined  in  this  analysis,  and  the 
acceptable  damage  level  arbitrarily  established  as  being  that 
damage  associated  with  80  percent  of  the  reflected  impulse 
necessary  to  overturn  each  vehicle.  With  such  forces  imparted 
to  the  vehicle,  some  deformation  of  the  exposed  sidewalls, 
window  breakage,  and  the  dislocation  of  some  internal  fix¬ 
tures  will  probably  occur,  especially  for  the  smaller  charges. 
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These  effects,  however,  should  not  present  an  unacceptable 
hazard  to  the  occupants. 

The  bus  analyzed  is  a  rather  standard  intercity  highway  bus 
(Figure  5) .  The  external  skin  sections  of  the  bus  are  attached 
to  a  rather  light  but  extensive  space  frame  composed  of  many 
small  members,  primarily  square  or  rectangular  steel  tubing. 
There  is  no  main  frame  such  as  is  commonly  found  in  trucks  and 
older  cars.  The  outer  skin,  either  steel  or  aluminum,  and  the 
plywood  passenger  and  baggage  floors  are  rigidly  attached  to 
the  framework  and  provide  part  of  the  structural  strength. 

The  empty  bus  weighs  between  18,000  and  22,000  pounds,  and 
when  fully  loaded  with  40  people  and  baggage,  an  estimated 
weight  of  about  30,000  pounds  is  attained.  The  center  of 
gravity  is  estimated  from  40  to  46  inches  above  the  ground 
plane . 

The  reflected  impulse  needed  to  overturn  a  stationary  target 
requires  a  determination  of  the  mass  distribution  of  the  tar¬ 
get  and  the  location  of  the  center  of  gravity.  The  height  of 
the  center  of  gravity  above  the  point  of  rotation.  A,  in  the 
ground  plane  is  designated  as  h  and  the  distance ,  in  the 
ground  plane,  from  the  center  of  gravity  to  the  point  of  rota¬ 
tion  as  d.  The  distance,  dQ  =  /d^+hg^-hg /  represents  the 
distance  the  center  of  gravity  must  rise  so  that  it  is  directly 
above  point  A,  at  which  point  gravitational  forces  will  over¬ 
turn  the  target.  For  most  targets,  point  A  is  a  point  in  the 
ground  plane  directly  below  the  outside  surface  plane  of  the 
vehicle . 

The  work,  W,  done  in  overturning  the  target  is 

W  =  dQ  •  (weight  of  the  vehicle,  w)  .  (11) 

When  a  sufficient  impulse  is  applied  to  the  target  rapidly, 
it  will  give  the  target  an  angular  velocity,  w,  great  enough 
to  permit  inertial  forces  to  complete  the  overturning  action. 
The  angular  velocity  will  be  sufficient  when  the  kinetic 
energy  is  greater  than  the  work  required  from  equation  (11) , 
where 

KE  =  1/2  IAo>2  .  (12) 


32 


GLASS 


33 


The  value  IA  is  the  moment  of  inertia  about  point  A  and  is 
given  as 

( b2+h2  2  ^ 

*A  =  m  +  c  j  (13) 

where : 

m  =  mass  of  the  target 
b  =  width  of  the  target 
h  =  height  of  the  target 

c  =  the  transfer  axis  distance;  c  =  dQ  +  hg  . 

By  equating  the  required  work,  W,  with  the  kinetic  energy  and 
substituting  equation  (13)  into  (12) ,  an  expression  for  the 
required  angular  velocity,  w,  to  overturn  the  target  is 
obtained: 


a)  =  [2W/Ia] 


1/2 


2  dQw 
JA 


„  V2 


(14) 


The  required  unit  impulse  (psi-msec) ,  H,  required  to  produce 
this  angular  velocity  is 


1000 • Ia*  w 

^  hc- (presented  area  of  the  target)  (15) 

where  hc  represents  the  height  above  the  ground  where  the  cen¬ 
ter  of  the  blast  pressure  is  applied. 

The  computed  reflected  impulses  required  to  overturn  the  sta¬ 
tionary  bus  were  all  well  within  the  present  quantity-distance 
criteria  for  the  five  charge  sizes  considered  in  Reference  1. 
The  bus  therefore  was  analyzed  under  more  vulnerable  condi¬ 
tions,  that  of  being  in  a  460-foot  radius  turn  at  50  mph. 

In  order  to  select  a  typical  turning  condition  for  the  bus, 
it  was  experimentally  determined,  that  for  normal  safe  driving, 
the  resultant  force  exerted  on  a  vehicle  by  the  road  (composed 
of  the  centripetal  force  applied  in  the  horizontal  plane  and 
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the  weight  reaction  force  applied  vertically)  would  be 
inclined  a  maximum  of  20  degrees  from  the  vertical.  This 
would  occur  at  different  speeds  for  different  radius  curves, 
but  at  50  mph  is  associated  with  a  460-foot  radius  curve.  The 
centripetal  force  corresponding  to  20  degrees  is  w*tan  20 
degrees  =  .364  w.  A  condition  of  "unstable  equilibrium"  can 
be  used  which  is  represented,  for  a  stationary  bus,  by  a 
tipping  to  the  point  where  the  center  of  gravity  is  directly 
above  point  A;  for  a  bus  in  a  turn,  unstable  equilibrium  will 
occur  when  a  blast  on  the  inner  side  of  the  curve  has  tipped 
the  bus  such  that  the  sum  of  the  moments  about  A,  of  the 
centrifugal  force  and  the  weight,  equal  zero.  Thus,  the  bus 
in  a  50-mph,  460-foot  radius  turn  is  at  unstable  equilibrium 
when : 


.364w  ( /(j2+j1^2  •  sin0 )  -  w  ( /d2+h^2cos0 )  =  0 
. 364w  =  cot 0  (16) 

0  =  70  degrees 

where  0  is  the  angle  between  the  ground  plane  and  a  line  con¬ 
necting  the  displaced  center  of  gravity  with  the  point  of  rota¬ 
tion  A,  and  d  and  hg  are  as  previously  defined.  Therefore, 
gravity  must  rise  a  distance 

dQ  =  (sin0/d2+hg2)  -  hg  . 

The  distance  the  c.g.  moves  in  the  direction  of  the  centrifu¬ 
gal  force  becomes 

dc  =  cos0  /d2+hg2  . 

The  work  done  by  the  centrifugal  force  is  therefore 

Wc  =  dc  (.364  w)  .  (17) 

The  above  work  is  subtracted  from  the  total  work  requirement 
calculated  with  equation  (11)  to  determine  the  work  required 
only  by  the  blast  force  to  overturn  the  bus.  Using  this  modi¬ 
fied  work  calculation  in  equations  (13) ,  (14) ,  and  (15)  results 

in  a  determination  of  the  reflected  impulse  calculation  required 
to  overturn  the  bus  moving  in  the  460-foot  radius  turn. 
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In  the  computer  analysis  the  location  of  the  c.g.  of  the  bus 
and  its  weight  are  treated  as  random  variables.  Up  to  four 
values  of  these  parameters  are  entered  on  cards  55  and  56 . 
The  net  loading  concept  is  used  in  determining  the  reflected 
impulse  imparted  to  the  bus  for  each  specified  quantity- 
distance  combination  and  set  of  parameter  values.  As  in  the 
treatment  of  walls,  the  "effective  impulse"  is  assumed  to  be 
that  impulse  obtained  by  integrating  the  net  overpressure 
function  to  a  point  at  which  it  is  £  .2  psi.  This  concept 
was  used  for  all  vehicular  targets. 


G.  CAMPER-PICKUP  UNIT 

The  structural  detail  of  the  camper-pickup  unit  is  illustrated 
in  Figure  6.  The  side  walls  are  composed  of  a  very  light, 
non-structural ,  aluminum  outer  skin;  a  very  irregular  glued 
and  stapled  framework  of  1-  by  2-inch,  2-  by  2-inch,  and  2- 
by  3-inch  wood  members  with  spacing  to  accommodate  the  windows, 
cupboards,  etc.;  and  an  inner  covering  of  1/8-  or  3/16-inch 
plywood  glued  and  nailed  to  the  framework.  In  the  interior 
there  are  numerous  cupboards,  seats,  counters,  etc.,  which 
contribute  considerable  strength  and  stiffness  to  the  side 
walls . 

The  minimum  separation  distance  for  this  unit  for  a  specific 
charge  was  determined,  as  previously  stated,  as  the  distance 
at  which  80  percent  of  the  reflected  impulse  required  to  over¬ 
turn  the  unit  is  attained.  Again,  the  permanent  deformation 
of  the  side  walls  and  glass  breakage  which  will  occur  at  this 
point  are  considered  acceptable. 

The  reflected  impulse  required  to  overturn  the  camper-pickup 
unit  can  be  determined  through  straightforward  use  of  equa¬ 
tions  (11)  through  (15)  presented  in  the  previous  section. 

Separate  estimates  of  the  camper  c.g.  location,  truck  c.g. 
location,  camper  weight,  and  truck  weight  are  used,  and  up  to 
four  values  each  are  specified  on  input  data  cards  56  through 
60.  The  net  blast  loading  concept  is  used  for  the  camper- 
pickup  unit  and  the  computed  reflected  impulse  achieved  for  a 
given  charge-distance  combination  computed  from  the  resultant 
net  loading  function. 
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H.  MOBILE  HOME 


The  mobile  home  (Figure  7)  has  a  main  frame  consisting  of  two 
light  but  deep  channels,  with  cross  channels  welded  between 
these  main  frame  members  in  the  same  plane.  On  top  of  the 
cross  channels,  2-  by  4-inch  wood  members  run  lengthwise  on 
16-inch  centers;  they  are  stabilized  by  1-  by  3-inch  wood  mem¬ 
bers  placed  at  right  angles,  on  24-inch  centers.  Covering 
these  are  1/2-inch  plywood  sheeting  and  non-structural  finish 
flooring  materials. 

The  side  walls  have  a  non-structural  outer  metal  skin,  2-  by 
3-inch  studs  on  16-inch  centers  with  several  1-  by  2-inch 
stringers ,  and  an  inner  covering  of  nailed  and  glued  3/16-inch 
plywood.  The  various  cupboards,  partitions,  etc.,  have  not 
been  illustrated.  The  roof  has  light  truss  cross  members  com¬ 
posed  of  2-  by  2-inch  chords  and  metal  web  members.  There  is 
a  2-  by  4-inch  wood  member  running  lengthwise  through  the 
center  of  the  trusses  to  stabilize  them.  The  roof  is  not 
particularly  vulnerable  to  blast  loading  because  the  shock 
wave  is  assumed  to  impinge  on  the  side  of  the  trailer  and  the 
reflected  pressure  on  the  roof  is  only  very  slightly  higher 
than  the  incident  pressure. 

Up  to  four  values  of  the  c.g.  location  and  mobile  home  weight 
must  be  entered  on  input  data  cards  61  and  62.  The  net  load¬ 
ing  concept  is  also  used  in  determining  the  reflected  impulse 
imparted  to  the  mobile  home  for  a  specified  quantity-distance 
combination. 


I.  EXPLOSIVE  STORAGE  IGLOO 

The  explosives  storage  igloo  is  supported  by  a  No.  1  gage  cor¬ 
rugated  steel  arch  with  a  2-foot  minimum  earth  cover.  The 
igloo  has  a  13-foot  radius  and  is  59  feet  long  (inside  dimen¬ 
sion)  .  The  front  concrete  wall  is  12  inches  thick  with  two 
5-  by  10-foot  steel  doors  opening  outward.  The  outside  of 
the  door  is  covered  with  3/8-inch  steel  plate.  Insulation 
and  16-gage  steel  plate  comprise  the  inside  portions  of  the 
door.  This  igloo  construction  is  shown  in  Figure  8. 

The  3-inch  thick  igloo  door  is  composed  of  a  3/8-inch  outer 
steel  plate  welded  to  a  3-inch  deep,  6-lb/ft  steel  channel 
forming  the  perimeter  of  the  door.  One  inch  of  insulation  is 
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Figure  7.  Mobile  Home 
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Figure  8.  Standard  Explosive  Storage  Igloo,  Front 


glued  to  the  inside  of  a  3/8-inch  plate,  and  a  16-gage  plate 
is  tack-welded  on  the  backside  of  the  channel.  The  door  is 
held  in  place  by  three  hinges  with  1-inch  diameter  pins.  Pins 
within  the  door  of  1-3/8-inch  diameter  extend  1-1/2  inches 
into  reinforced  holes  in  the  door  frame. 

The  acceptable  damage  levels  were  defined  as  severely  crack¬ 
ing  the  igloo  concrete  wall  or  dislodging  the  igloo  doors-- 
whichever  occurs  first.  In  Reference  1  three  types  of  failure 
were  evaluated;  these  are 

1.  Deformation  of  the  door  such  that  the  pins  pulled 
from  their  sockets  in  the  door  frame, 

2 .  Shearing  of  the  pins  between  the  door  and  the 
socket,  and 

3.  Failure  of  the  concrete  wall  around  the  door. 

Analyses  of  shear  pin  failure  and  of  fracturing  the  concrete 
wall,  by  Johansen's  yield  line  theory,  indicate  that  deforma¬ 
tion  of  the  door  so  that  the  pins  are  dislodged  will  occur 
first . 

An  8-inch  deflection  of  the  door  is  necessary  in  order  to 
dislodge  the  pins  from  the  holes  in  the  door  frame.  Analyzed 
as  two  rigid  halves  with  a  plastic  hinge  in  the  middle,  the 
work  W  required  to  achieve  the  stipulated  deflection  becomes 
W  =  M0 ;  where  M  is  the  plastic  moment  needed  to  deform  the 
central  section  of  the  door,  and  0  is  the  central  deflection 
angle.  The  required  unit  impulse  is  then  determined  by  using 
equation  (10) .  Values  of  Young's  modulus  of  the  steel  door 
are  entered  on  input  data  card  number  63. 

The  reflected  impulses  which  result  in  deformation  of  the 
steel  door,  so  that  the  pins  are  dislodged,  are  used  as  a 
basis  for  establishing  quantity-distance  specifications  in 
this  program.  Using  E  =  5  x  10^  psi  in  Reference  1,  a  criti¬ 
cal  impulse  of  139  psi-msec  was  obtained. 

J.  PERSONNEL 

The  "standard  man"  as  used  in  several  military  target  descrip¬ 
tions  and  vulnerability  studies  is  illustrated  in  Figure  9. 
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Figure  9.  Standard  Man 
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The  quantity-distance  specifications  for  the  personnel  target 
have  been  established  so  as  to  prevent  his  being  thrown  to 
the  ground,  down  stairs,  or  against  nearby  structures.  Love¬ 
lace  Foundation  studies  [8]  have  indicated  that  the  transla¬ 
tion  velocity  at  which  man  begins  to  be  subjected  to  these 
dangers  is  approximately  3.5  fps.  The  "acceptable  damage" 
criterion  thus  is  based  on  the  reflected  impulse  which  results 
in  the  translation  of  a  man  at  a  maximum  velocity  of  3.5  fps. 

Dynamic  overpressure  is  the  main  factor  in  the  translation  of 
objects  in  a  blast  wave.  In  order  to  translate  a  man  of  given 
weight  and  presented  area,  the  impulse  associated  with  the 
integration  of  the  dynamic  overpressure  function  must  be 
greater  than  the  product  of  the  mass  of  the  man  and  his  trans¬ 
lation  velocity  or 


(18) 


PA 


where : 


PA  =  presented  area  of  the  man  (in^) 


q(t)  =  dynamic  overpressure  function 
m  =  mass  of  the  man 
v  =  translation  velocity. 

The  above  equation  allows  velocity  to  be  expressed  as  a  func¬ 
tion  of  time;  the  computer  program  integrates  that  function 
to  obtain  a  translation  distance  associated  with  a  maximum 
velocity  of  3.5  fps.  The  maximum  acceleration  amx  over  the 
calculated  translation  distance  becomes 


(19) 


where : 


Fmx  =  force  given  by  Q0*Pa  ,  and 


Q0  =  peak  dynamic  overpressure. 
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Up  to  four  values  of  the  man's  weight  and  his  presented  area 
are  entered  on  cards  64  and  65.  In  the  analysis  a  peak  reflec¬ 
ted  overpressure  of  5.0  psi  is  used  as  the  threshold  of  eardrum 
rupture.  In  cases  where  only  eardrum  hazards  occur,  the  maxi¬ 
mum  distance  from  the  charge  where  5.0  psi  reflected  pressure 
occurs  is  designated  as  the  separation  distance  criterion. 


K.  COMMERCIAL  AIRCRAFT 

The  Boeing  707  aircraft  is  used  as  the  representative  large 
commercial  jet  aircraft.  A  detailed  drawing  of  each  section 
shown  in  Figure  10  is  presented  in  Reference  1. 

The  fuselage  shell  is  of  the  semi-monocoque  type  with  aluminum 
skin  and  clad  7075-T6  longitudinal  stringers.  The  shell  is 
stiffened  circumferentially  by  bulkheads  and  7075-T6  "Z"  sec¬ 
tion  frames.  The  wing  is  a  full-cantilever,  semi-monocoque, 
cellular  structure  tapering  in  both  platform  and  depth.  The 
skin  coverings  vary  in  thickness  from  0.06  to  0.45  inch.  The 
horizontal  stabilizer  and  vertical  fin  of  the  tail  are  of 
two-spar  construction.  There  is  no  other  spanwise  stiffening. 

The  University  of  Dayton  Research  Institute  has  conducted  an 
extensive  study  of  the  static  and  dynamic  loadings  required 
to  fail  the  most  vulnerable  skin  panels  and  other  components 
on  this  aircraft  [9] .  The  first  step  in  their  study  was  to 
perform  stress  analyses  on  several  of  the  weaker  fuselage  skin 
panels.  Since  the  skin  panels  transmit  the  blast  loading  to 
the  internal  structure,  the  location  of  the  critical  skin 
panels  gives  an  indication  of  areas  in  which  the  substructure 
should  be  investigated.  Complete  analyses  were  conducted  on 
these  components.  Based  on  their  conclusions,  a  large  fuselage 
skin  panel  located  between  fuselage  stations  1060  and  1080  on 
the  underside  of  section  C  fuselage  shown  in  Figure  10  was 
shown  to  be  the  most  critical  or  vulnerable  skin  panel.  Per¬ 
manent  deformation  of  this  panel  occurs  under  a  2.3  psi 
reflected  overpressure. 

It  is  highly  probable  that  the  aircraft  can  undergo  serious 
deformation  of  many  skin  panels  without  approaching  a  state 
causing  loss  of  control.  Therefore,  threshold  deformation 
of  the  substructure  supporting  the  weakest  skin  panel  has 
been  defined  as  the  acceptable  damage  level. 
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Figure  10.  Boeing  707  Aircraft 


The  static  loading  required  to  permanently  deform  the  frame 
section  supporting  the  vulnerable  skin  panels  described  above 
is  7.23  psi.  The  dynamic  response  of  this  panel  is  described 
as  a  function  of  the  ratio  of  the  positive  duration  of  the 
blast  wave  to  the  natural  period  of  the  frame  member.  The 
natural  period  of  the  frame  member  T  is  given  by: 

T  =  2JI/n  ,  (20) 

where : 


n  =  2.68  -  (21) 

V  ApR 

where,  as  in  (21) ,  and  as  previously  defined: 

E  =  Young's  modulus 

I  =  cross-sectional  moment  of  inertia 
g  =  gravitation  constant 
A  =  cross-sectional  area 
p  =  density  of  the  frame  material 
R  =  radius  of  the  frame. 

The  relationship  between  the  ratio  of  the  positive  duration 
and  the  natural  period  of  the  frame  and  the  dynamic  loading 
factor,  a,  is  shown  in  Reference  1,  page  94  and  Reference  9. 
Ten  values  are  entered  on  input  card  67.  These  correspond  to 
the  range  l/2<a<5  in  1/2  unit  increments.  The  ten  values 
listed  under  the  card  format  definition  must  be  used. 

Up  to  four  estimates  of  E  for  the  aluminum  frame  member  may 
be  entered  on  card  66.  These  values  lie  in  a  range  of  10.0 
to  12.0  million  psi.  For  E  =  10.6  x  10^  psi,  a  natural  period 
of  36  msec  is  calculated  for  the  frame  member. 

A  quantity-distance  specification  for  each  charge  size  is  com¬ 
puted  as  the  distance  at  which  the  incident  overpressure,  PQ , 
in  (21)  satisfies  the  inequality, 

DF  •  (RF*P0)  >  7.23  psi  (22) 
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where : 


Dp  =  dynamic  response  factor  and 
Rp  =  reflection  factor. 

For  a  normal  side-on  impingement  of  the  blast  wave  on  the  fuse¬ 
lage,  the  reflection  coefficient  for  the  panel  under  investiga¬ 
tion  is  determined  from  an  angle  of  incidence  of  approximately 
63  degrees. 
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VI.  FLOW  CHARTS  AND  SOURCE  PROGRAM 


Flow  charts  and  a  source  deck  listing  of  this  program  are  pro¬ 
vided  in  this  section.  Block  numbers  and  subtitles  have  been 
used  liberally  through  the  flow  charts  so  as  to  facilitate  an 
understanding  to  the  desired  computation.  Blocks  applicable 
to  each  subroutine  are  shown  in  Table  II. 
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INPUT  DATA  SUBROUTINE  -  HEADR 


o 


0 


O 


Figure  11.  Flow  Charts  of  the  Quantity-Distance  Computer 
Program 
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11.  TARGET  SELECTION  SUBROUTINE  -  CALLR 


Figure  11  (Continued) 
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III.  SPLIT-LEVEL  HOUSE  ANALYSIS  SUBROUTINE  -  HOUSE 


FIXED  PARAMETERS  ARE  SPECIFIED 
FEET  -  LENGTH  OF  ROOF  RAFTER  (£) 

AI  -  MOMENT  OF  INERTIA  OF  RAFTER  <Ia) 
C  -  DISTANCE  FROM  NEUTRAL  AXIS  TO 
EXTREME  FIBER  (C) 

G  -  GRAVITATIONAL  CONSTANT  {g) 


* ‘DETERMINES  NUMBER  OF  INPUT  PARAMETERS** 


•‘COMPUTES  CRITICAL  PERIODS  AND  CRITICAL  IMPULSES  FOR  HOUSE  ROOF** 


Figure  11  (Continued) 
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* 'COMPUTES  CRITICAL  PERIODS  AND  CRITICAL  IMPULSES  FOR  HOUSE  FRAME  SECTION** 


Figure  11  (Continued) 
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MODULUS  OF  RUPTURE 
0  -  HFW(III) -10000 


CRITICAL  PERIOD 


TC2 (1,11,111) 


2-1000 


CRITICAL  IMPULSE 

UI2d.II, III) 


* ‘COMPUTES  CRITICAL  PERIODS  FOR  HOUSE  BRICK  WALL** 


* ‘DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE  SIZE  AND  SET  OF  HOUSE  ROOF  PARAMETERS** 


Figure  11  (Continued) 
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♦VARIATION  OF  50  SCALED  DISTANCES,  *D* 


Figure  11  (Continued) 
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••MINIMUM  DISTANCE  TO  DAMAGE  HOUSE  ROOF  WITH  GIVEN  CHARGE  SIZE  IS  SELECTED* * 


* ‘CRITICAL  IMPULSE  AND  COMPUTED  IMPULSE  ON  HOUSE  FRAME  SECTION  ARE 
COMPARED  FOR  THE  SELECTED  DISTANCE  AND  GIVEN  CHARGE  SIZE** 


•BLAST  LOADING  ON  OUTSIDE  OF  FRAME  SECTION* 


Figure  11  (Continued) 
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•blast  loading  on  rear  surface  of  frame  section* 


* ‘CRITICAL  IMPULSE  ANO  COMPUTED  IMPULSE  ON  HOUSE  BRICK  WALL  ARE 


COMPAREO  FOR  SELECTED  OISTANCE  AND  GIVEN  CHARGE  SIZE** 


* ‘CRITICAL  IMPULSE  FOR  HOUSE  BRICK  WALL  IS  FUNCTION  OF  BLAST 
LOADING  ON  ROOF  OF  HOUSE** 


Figure  11  (Continued) 
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IV.  CHURCH  BUILDING  ANALYSIS  SUBROUTINE  -  CHURCH 


•DETERMINES  NUMBER  OF  INPUT  PARAMETERS* 


Figure  11  (Continued) 
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•COMPUTES  CRITICAL  PERIODS  AND  CRITICAL  IMPULSES  FOR  CHURCH  ROOF  BEAM* 


Figure  11  (Continued) 
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** DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE  SIZE 
AND  SET  OF  CHURCH  BEAM  PARAMETERS** 


♦VARIATION  OF  FIVE  CHARGE  SIZES* 


Figure  11  (Continued) 


♦MINIMUM  DISTANCE  TO  DAMAGE  CHURCH  ROOF  BEAM  WITH  GIVEN  CHARGE  IS  SELECTEO* 


♦♦CRITICAL  IMPULSE  ANO  COMPUTED  IMPULSE  ON  CHURCH  DECKING  ARE  COMPARED 
FOR  THE  SELECTEO  OISTANCE  ANO  GIVEN  CHARGE  SIZE** 


Figure  11  (Continued) 
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V.  SCHOOL  BUILDING  ANALYSIS  SUBROUTINE  -  SCHOL 


Figure  11  (Continued) 
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•DETERMINES 


NUMBER  OF  INPUT  PARAMETERS* 


•COMPUTES  CRITICAL  PERIOD  FOR  THE  SCHOOL  BRICK  WALL* 


Figure  11  (Continued) 
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♦COMPUTES  CRITICAL  PERIODS  AND  CRITICAL  IMPULSES  FOR  SCHOOL  ROOF1 


♦♦DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE  SIZE 
AND  SET  DF  SCHOOL  BRICK  WALL  PARAMETERS* * 


Figure  11  (Continued) 


•CRITICAL  IMPULSE  FOR  SCHOOL  BRICK  WALL  IS  FUNCTION  OF  THE 
BLAST  LOADING  ON  ROOF  OF  THE  SCHOOL* 


•BLAST  LOADING  ON  SCHOOL  BRICK  WALL* 


Figure  11  (Continued) 
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TS<CP<TEND 

CP<TS<DUR 

COMPUTE  AVERAGE  OVERPRESSURE 

_ y 

IF (CP -TEND) 

> 

tBfmf  rlH 

x> 

COMPUTE  AVERAGE  OVERPRESSURE 

OVER  TIME  INTERVAL  0 < t< CP : 

IN  TIME  INTERVAL  0<_t<.CP: 

AVER 

AVER 

CALCULATE  REFLECTED  IMPULSE 
IN  TIME  PERI 00  0< t<T£ND : 
CALI 


WRITE  CHARGE  SIZE,  RANGE 
CRITICAL  IMPULSE,  PI, 
PR,  AND  DUR 


•CRITICAL  IMPULSE  ANO  COMPUTED  IMPULSE  ON  SCHOOL  ROOF  ARE 
COMPARED  FOR  SELECTED  OISTANCE  AND  GIVEN  CHARGE  SIZE** 


WRITES  PROBABILITY  OF  GLASS  FRAGMENT 
INJURY  STATEMENT 


Figure  11  (Continued) 
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VI.  OFFICE  BUILDING  (HALF  WALL  KITH  WINDOWS)  ANALYSIS  SUBROUTINE  -  OHAFW 


© 


©" 


FIXED  PARAMETERS  ARE  SPECIFIED 
FEET  =  HEIGHT  OF  HALF  WALL  (h) 

DI  =  LENGTH  OF  HALF  WALL  (i) 

TW  -  THICKNESS  OF  WALL 

AIP  -  MOMENT  OF  INERTIA  OF  WALL  (Ia 

G  »  GRAVITATIONAL  CONSTANT  (g) 


•DETERMINES  NUMBER  OF  INPUT  PARAMETERS  * 


Figure  11  (Continued) 
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♦♦DETERMINES  MINIMUM  SEPARATION  DISTANCE  FDR  EACH  CHARGE  SIZE 
AND  SET  OF  OFFICE  HALF  WALL  PARAMETERS  * ♦ 


Figure  11  (Continued) 


* VARIATION  OF  FIVE  CHARGE  SIZES* 


•VARIATION  OF  50  SCALED  DISTANCES* 


•BLAST  LOADING  ON  REAR  SURFACE  OF  HA  F  WALL* 


Figure  11  (Continued) 


VII.  DFFICE  BUILDING  (FULL  WALL  WITH  NO  DPENING)  -  OFFICE 


Figure  11  (Continued) 
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* 'DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE  SIZE 


AND  SET  OF  OFFICE  FULL  WALL  PARAMETERS*  * 


* VAR I AT  I  V  OF  FIVE  CHARGE  SIZES* 


Figure  11  (Continued) 
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VIII.  PASSENGER  BUS  OVERTURNING  SUBROUTINE  -  BUS 


*OET ERMINES  NUMBER  OF  INPUT  PARAMETERS* 


Figure  11  (Continued) 
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** DETERMINES  CRITICAL  IMPULSES  FOR  BUS  OVERTURNING** 


♦DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE 


SIZE  AND  SET  OF  BUS  PARAMETERS* 


Figure  11  (Continued) 
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♦BLAST  LOADING  ON  REAR  SURFACE  OF  BUS* 


Figure  11  (Continued) 
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IX.  CAMPER  AND  PICKUP  TRUCK  OVERTURNING  SUBROUTINE  -  CAMPR 


•DETERMINES  NUMBER  OF  INPUT  PARAMETERS* 


Figure  11  (Continued) 
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** DETERMINES  MINIMUM  SEPARATION  DISTANCES  FOR  EACH  CHARGE 


SIZE  AND  SET  OF  CAMPER-PICKUP  PARAMETERS*  * 


DO  424  I  -  1,  K1 


CG  LOCATIONS 
CGT  -  CTG(I) 
CGC  -  CCG(I) 


DO  424  J  ■  1 ,  Kl 

WEIGHT  OF  UNIT 
WTR  -  TRW  ( J) 
WCP  -  CPW(J) 


WRITES  OF  KI2  SETS  OF 

CAMPER-PICKUP  PARAMETERS 


•VARIATION  OF  FIVE  CHARGE  SIZES* 


Figure  11  (Continued) 


76 


•BLAST  LOADING  ON  REAR  SURFACE  OF  CAMPER  IMCKUI  UNIT* 


Figure  11  (Continued) 
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X.  MOBILE  HOME  OVERTURNING  SUBROUTINE  -  TRUER 


Figure  11  (Continued) 
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* ‘DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE 
SIZE  AND  SET  OF  MOBILE  HOME  PARAMETERS** 


Figure  11  (Continued) 
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‘BLAST  LOADING  ON  REAR  SURFACE  OF  TRAILER* 


XI.  IGLOO  DOOR  ANALYSIS  SUBROUTINE  -  DOOR 


Figure  11  (Continued) 
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•DETERMINES  CRITICAL  IMPULSE  FOR  IGLOO  DOOR* 


♦‘DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE  SIZE 
AND  SET  OF  IGLOO  DOOR  PARAMETERS** 


♦VARIATION  OF  FIVE  CHARGE  SIZES* 


Figure  11  (Continued) 
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•VARIATION  OF  50  SCALED  DISTANCES* 


XII.  MAN  TRANSLATION  SUBROUTINE  -  MAN 


** DETERMINES  NUMBER  OF  MAN  TRANSLATION  PARAMETERS** 


Figure  11  (Continued) 
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“DETEKMTNES  minimum  separation  distance  for  each  charge  size 

AND  SET  OF  MAN  TRANSLATION  PARAMETERS** 


Figure  11  (Continued) 
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♦VARIATION  OF  SO  SCALEO  OISTANCES* 


XIII.  COMMERCIAL  JET  AIRCRAFT  ANALYSIS  SUBROUTINE  -  AIRP 


Figure  11  (Continued) 
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*  *  DETERMINES  FULL  NATURAL  PERIODS  OF  FRAME  MEMBER** 


* ‘DETERMINES  MINIMUM  SEPARATION  DISTANCE  FOR  EACH  CHARGE  SIZE 
AND  SET  OF  AIRCRAFT  PARAMETERS** 


•VARIATION  OF  FIVE  CHARGE  SIZES* 


♦VARIATION  1  50  SCALED  DISTANCES* 


Figure  11 


(Continued) 
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Figure  11  (Continued) 
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HF  ADR 


//  JOR 
//  DUP 
♦DFLFTE 
//  FOR 

*ONF  word  intfgfrs 

* AR I THMFT I C  TRACF 
♦TRANSFER  TRACF 
*  IOCS ( CARD. n 7?  PRINTER) 

COMMON  PO (SO )  , T AU 1  SO )  ,CPR (SO)  , HRE 1 4 )  , HRW ( 4 )  ,HRS( A )  ,HFF ( A )  , 

*HFW  <  4  )  »HFS( 4 ) »HBE<4) »HBT (4 ) »CRE  <4 ) »CRW(4 ) »CRS<4 ) »CDE  *4 ) »CDW(4 ) • 

♦  CDS(4) *SBWE( A) »SBT ( A ) *SRE( A) *SRW( A) *SRS( A) »OHF( A) *OHT  <  A ) *OFE ( A) * 

* C  TG 1 4 ) ♦ T  RW 1 A )  ,CPW ( A  )  .TRCGl A )  »  TRWT (A)  .DORS ( A ) *  WMAN  <  A ) *  PM AN ( A )  * 

*OFT (4 ) » I A( 11 ) »K  »NOT  »UI ( 4*4  *4) »  TC ( 4*4*4) *CS( S ) *RCG ( 4) »BW (4 ) *CCG ( A ) » 

*  A  I R  F  (  4 ) *AIRF( 10 ) *TC2(4*4*4)  *  T  C  S  (  4  *  4  *A) »  U  1 7  ( A*A*A) *  U  I  S  ( A*A*A) 

C  . 

C  PART  1.  INPUT  DATA  SUBROUTINE 

C  . 

C . CHARGF  SIZES . 

READ! 2  » 1  ) <CS( T ) *  1  =  1 *S ) 

C . INCIDENT  OVERPRESSURES.DURATIONS.AND  reflection  factors . 

RFAD1 7,1  ) (P0< I  )  ,  1  =  1 ,S0) 

READ!  7»  1  )  <  T  AU  <  I  >  »  I  =  1  .SO ) 

READ! 7,1)  ( CPR (  I  )  »  T  =  1  .SO) 

C . HOUSF  ROOF  PARAMETERS . 

R  FAD ( ?»?  )  (HRE (  I  )  .  1  =  1  .4  ) 

READ! ?,?  )  ( HRW (  I  ) .  1  =  1 .4  ) 

RFAD( 7.2)1  HRS (  I ) ,  1  =  1  , A ) 

C . HOUSE  FR  AM  F  SFCTIDN  PARAMETERS . 

READ! ?.7)  <HFF(  I ) *  1=1  *4 ) 

PFAD1 7 • 7 )  1 HFW 1  I) *  1  =  1 *4 ) 

RFAOl  r  »  ?  )  ( HF  S (  I  ),  1  =  1  ,4) 

C . HOUSE  PRICK  WALL  SFCTION  PARAMFTFPS . 

RFAD1 7. 7 )  (HPF  1  I  )  .  1  =  1  ,4 ) 

RFAD1?.?)  (  HR  T  (  I  )  »  I  *  1  ♦  4  ) 

C . CHURCH  ROOF  PARAMFTFRS . 

R  F  ADI  7.2 )  1  CRF  1  1  ) , I  =  1  ,  A  ) 

READ!  7 .7  )  1  C RW  1  I  )  ,  I  =  1  » A  ) 

READ!?. 7)  (CRS(I).  1  =  1.4) 

C . CHURCH  DECKING  PARAMETERS . 

R  F ADI  7. 7 )  1  CDF (  I ) ,  1  =  1  ,4 ) 

READ! 7.7)  1 CDWl  I ) ,  1  =  I  ,  A  ) 

RFAD1 7. 7 )  1  CDS!  I ) ,  I  =  1  ,  A  ) 

C . SCHOOL  BRICK  WALL  PARAMETERS . 

RFAD17.2)  1 SRWF 1  I  ) .  I =1  . A ) 

READ!  7.7)  1  SRT 1  I  )  ,  I  =  1  ,4 ) 

C . STHOOL  ROOF  PARAMETERS . 

READ(  ?, ?  )  ( SRE (  I  ) ,  I =1  ,4 ) 

READ!?,?)  1 SRW ( I ) » I = 1 >4 ) 

RFAD17.7) 1 SRS ( I ) ,1=1 ,4) 

C . OFF  I CF  BUILDING  HALF  WALL  PARAMETERS . 

RFAD1 7,7 )  1 OHF ( I ) , 1=1 ,4 ) 

R  F  ADI  7,2)  1 OHT (  I  ), 1  =  1  ,4) 

C . OF  F  I  CF  FULL  WALL  PARAMFTFRS . 

RFAD1 ?, 7 )  (OFF  1  I  )  ,  I  =  1  ,4 ) 

READ! ?,? )  (OFT!  I  )  ,  1  =  1  ,4) 

C . BUS  PARAMFTFRS . 

RFAD17.7)  1 BCGl I ) » I =1 ,4 ) 
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RFADt  2*2)  <BW< I) , 1=1,4) 

C . CAMPER  AND  PICKUP  TRUCK  PARAMETERS . 

R  F  AD(  ?  »  2  )  <CCG<  I  )  ,  1  =  1  ,4) 

R  FAD(  2*2  )  <CTG(I)  *  1=1  *4) 

RFADt  2.2 )  <  TRW ( I ) . 1  =  1 .4) 

RFADt  2*2  )  (CPW<  I )  , 1  =  1  ,4) 

C . MOBILE  HOME  PARAMETERS . 

R  FAD( 2  *  2 )  I TRCG (11*1=1*4) 

RFADt  2*2)  ( TRWT (  I)  , 1  =  1 ,4 ) 

C . IGLOO  DOOR  PARAMETERS . 

RFADt  2.2 )  (DORSt  11,1=1.4) 

C . MAN  PARAMFTFRS . 

RFADt  2,2)  (WMAN (  I  ) ,  1  =  1  ,4  ) 

R  FADt  2*2)  ( PMAN ( I  ) » I = 1  *  4 ) 

C . AIRPLANF  PARAMETERS . 

R  FADt  2  *  2 )  ( AIRFt I ) » 1=1 ,4) 

READt  2*1)  ( AIRF ( I ) , 1  =  1 » 10) 

C . NUMBER  OF  TARGETS  CONSIDERED  AND  THEIR  ORDER . 

R  EAD  (  2  *  1 )  NOT*  I  At  1  )  *  I  A<  2  )  *  I A  (  3  )  ♦  I  At  4  )  »  I  A(S  )  »  I  At  6  )  *  I A  (7  )  *  I  A  (8  )  » 

* IA(9) , T A 1 10) , I  At  1 1  ) 

K  =0 

CALL  LINK ( C ALLP ) 

1  FORMAT ( 5F1 2. 2 ) 

2  F  ORM A  T ( 4  F 1 2 • 2 ) 

3  F  ORM A  T (1212) 

FND 

//  DUP 

♦STORF  WS  UA  HF  ADR 

//  JOR 
//  DUD 

♦DFL  FT  F  CALLR 

//  FOR 

♦ONE  WORD  INTFGFRS 
♦ARITHMETIC  TRACF 
♦TRANSFFR  TRACF 
♦IOCS (CARD, 1122  PRINTFR) 

C  PART  2.  TARGFT  SUBROUTINE  SFLFCT I  ON 

COMMON  PO(SO)  , T  AU  t  SO )  , C  PR ( S  0 ) *HRF  t  4 )  *HRW(4) »HRS(4)  »HFF<4) » 

♦  HFW  (  4)  *HFS  (  4)  »  HRE  <  4  )  »  HRT  (  4  )  *  CRE  <  4  )  »CRWU  )  »CRS<4  )  *CDE  <4  )  *CDW(  4  )  * 
*CDS  t  4 ) *  SBWE  <  4 )  *SBT  t  4 )  *SRE  <  4)  * SRW  <4 )  ,SRS ( 4)  »OHE  <  4)  .OUT ( 4 )  *OFE  <  4  )  » 
♦CTG ( 4 ) »  TRW  <  4 )  »CPW  t  4 )  »TRCG  <  4 )  »  TRWT (4 ) , DORS ( 4  )  .WMAN <  4) »  PMAN ( 4 )  * 

♦  OFT ( 4) *  I A< 11  )  »K  »NOT ,  U  I (4*4*4)  ,Tc(4»4»4)  *C S ( S ) *BCG( 4) *BW<4 ) »CCG( 4) * 
♦A  I RF ( 4 )  *AIRF(lO)*TC2(4*4*4)  *TC3 ( 4  *4  *4  )  »U I  2  t  4  *4  *4)  *UI 3 ( 4*4  *4) 

K  =K  +  1 

IF(K-NOT)  1,1,2 
2  CALL  FXIT 

C . 

C . SELECTION  OF  TARGET  MODEL . 

C . 

1  J8= I  A ( K ) 

GO  TO  (  100 ,200 , 300 ,400 , SOD ,600 , 7P0 ,800 ,900  *  1000 , 1 100  ),  J8 
100  CALL  LINK  (HOUSF) 

200  CALL  LINK(CHRCH) 

300  CALL  LINK(SCHOL) 

400  CALL  LINK(OHAFW) 

SOO  CALL  LINK(OFULL) 

600  CALL  LINK (BUS) 
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100  CALL 
800  CALL 
900  CALL 

1  non  call 
1100  CALL 
END 

//  DUP 
♦STORE 
//  JOB 
//  DtJP 
*  DF  L  FT  F 
//  FOR 

♦ONF  WORD  IN  TFGFRS 
♦ARITHMETIC  TRACF 
♦TRANSFER  TRACF 
♦+OCS (CARP* 1137  PRINTFR) 

COMMON  PO(SO) ,  T  AU (SO)  *CPR  <  80 )  *HRF  <  4 )  *HRW( 4 )  *HRS ( 4 )  *  HFF ( 4 )  * 

*HFW (4)  tHFSU )  ♦  HRF ( 4 )  * HR T ( 4 )  *CRF ( 4 ) *CRW(4 )  *CRS< 4 )  *CDF  < 4  )  *CDW(4 )  * 
*CDS<4)*SrWf<4)  *SRT  <4)  »SRE<4) *SRW<4)  *SRS<4>  »OHFU)  *OhT<4)  *OFF < 4)  * 
*CTG<4>  *TRW<4>  *CPW  <  4  )  *  TRCG<  4  )  *  TRWT  < 4  )  *DORS  <  4  )  tWMAN<  4  )  *  PM  AN  <  4  )  * 

*OFT  <4  )  f  I  A<  11  )  tic  tNOT  tUl  <4  t4  *4)  *TC  <  tCS<  ■> )  ,  RCG  <  4  )'  *RW<  4  )  *CCG( 4  )  * 

*AIRE<4)tAIRF<l0)tTC2Ut4t4)tTC3<4t4t4)tUI?<4t4t4)tUI3<4t4t4) 


linkccampr  ) 
LINK ( TPLFR ) 
LINK(DOOR) 
LINK(MAN ) 
LINK! AIRP ) 


WS  UA  CALLR 


HOUSF 


PART  3.  HOUSE  ROOF  SURROUTINF 


C.  . 

c 

c . 

C . FIXFD  PARA^FTFR  VALUFS . 

A  =  1  •  /  3  • 

FFFT=?04. 

A  I  =  5  1  •  7 
C  =  3  •  7  5 
G  =  3  8  6  • 

C . DFTFRM I NES  NO.  OF  HOUSE  ROOF  PARAMETERS . 

DO  100  1=1*4 

I F ( HR  F (  I  )  )  101, 10]  ,100 

im  K  1  =  1—1 

GO  TO  14  6 
100  CONTINUF 
K  1  =4 

C . DETERMINES  NO.  OF  HOUSE  FRAME  WALL  PARAMETERS, 

146  DO  1?7  1=1*4 

I F  <  HFE  <  I  )  )  147*  147*  177 

147  K  7  = I- 1 

GO  TO  178 

177  CONTINUF 
K  7  =  4 

C . DFTFRMINFS  NO.  OF  HOUSF  B.  WALL  PARAMETERS..., 

178  DO  179  J=1 *4 

I F ( HRE ( J ) )  130,130, 179 

130  K3=J-1 

GO  TO  107 

179  CONTINUF 
K  3  =  4 


C . CRITICAL  PERIOD  AND  IMPULSE  FOR  HOUSE  ROOr, 

C . 

10?  DO  103  1=1  * K 1 

E=HRE ( I ) *1000000. 
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on  ini  1 1  =  1, ki 

W  T  =HRW ( I  I ) *8 • 5 
W=WI/102. 

DO  103  I  I  I =1 *K1 
S  I G  =HR S (  I  I  I  )*ioooo. 

PP=9. R7*SQRT ( (  E*A I *G ) / ( W*FEFT**4.  )  ) 

P=PP/5 .96 

T  C (  I ,  I  I ,  I  I  I  )  =  inon. / ( 4.*P ) 

WO= (SIG*R.*AI )/ (C*FFFT**?. ) 

U= ( ( WO**? , )*  (FFFT**5.))/(?40.*F*AI) 

wt=w*ffet 

H=SQRT(U*?.*WT/G) 

103  UI(  I*  I  I*  II  I  )  =  (H*1000. )  / 3 264 • 

WRITE! 3*3) 

WRITE! S,4) 

W  R  I  T  F  (  3 , 1  F  ) 

DO  1  n4  1=1 , K 1 
DO  104  J=1 *K 1 
DO  104  Jl=l ,K1 
E  =HRE (  I )*1000000. 

S I G =HR  S ( J]  )*100OO. 

104  WR I T E ( 3  *  3 )  E *HRW( J ) *SIG* TC(  I  * J* J]  )  *UI  (  I  * J* Jl  ) 

C . CRITICAL  PERIOD  AND  IMPULSE  FOR  HOUSE  FRAME  WALL 

DO  1  SI  I  =1  *K2 
F  =  hff ( i )*i0nooon, 

DO  1  SI  J  =  1  ,K? 

W  I  =  HFW( J ) *7 . S 
W  =W  I  / 90 , 

DO  1S1  I  I  1  =  ]  *  K  ? 

s  ig=hfs  (Hi  )*ioonn. 

PP=o. R 7* SORT ( ( F*s .  ??*G ) / ( W*qO. **4 . ) ) 

P=PP/5 .96 

TO? (  I  * J*  I  I  I  ) =] 000./ ( 4. »p ) 

WO=SIG*R .*6.??/ 1 4SRO. 

U  = (  (WO**? .)*(90.**5.)>/( ?40 ,*F*ft .?? ) 

WT=W*90. 

H  =SOR  T ( U*? .*WT/G ) 

1S1  UI? ( I , J,  I  I  I  ) =H* 1 OOO./l 440. 

WRITF! s, 9) 

WRITF! ^  *  4 ) 

WRI TF ( S  * 1R ) 

DO  is?  I =1 ,K? 

DO  1 S?  J  =  1  * K 2 
DO  IS?  Jl  =  l  * K ? 

E  =HFE(  I  ) *1000000. 

S I G=HFS ( J1 ) *10000. 

IS?  WRITF! *>  *  S  )  F *HFW( J) *SIG*TC? !  I  * J* J1  )  *UI? ( I  * J, J1  ) 

C . CRITICAL  PFRIOD  FOR  HOUSF  BRICK  WAIL . 

I  I  I  =1 

DO  1 S  S  1=1,  K3 
E  =  HRF (  I  )*1000000. 

DO  133  J  =  1 *  KS 
T  S  =  HRT ( J ) 

W  =  •  55 

P  P  =  9  * R 7* SORT (  ( E*4?.7*G ) / ( W*o0.**4 . )  ) 

P=PP/5 .96 

133  TCS ( I  * J*  I  I  I  ) =1000. / (4.*P ' 
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WRTTF(  ”3*10) 

WRI TEC  3* 1 1 ) 

WRI TFI It  19 ) 

DD  139  I =1 ,K3 
00  139  J=l,k3 
F=HRF( I i*ioooooo. 

139  WR I T  F ( 3  »  1  2  )  F.HRTt J)  ,TC3<  I , J, I  I  I  ) 

C . 

C . DFTFRM I NES  MINIMUM  STAND-OFF  STORAGE  D1STANCF  FOR  FACH 

C . CHARGE  SIZE  AND  SET  OF  PARAMETER  VALUES . 

C . 

C . VARIATION  OF  PARAMETER  VALUFS . 

DO  IDS  1=1  ,  K 1 
DO  105  J  =  1  * K 1 
DO  105  J  1  =  1  * K 1 
CI=UI( I* J* J1 ) 

CP  =  TC (  I. J. J1 ) 
f=hrei I )*ioooooo. 

S IG=HRS( J1 1*10000. 

WRITE(3>7)  E*HRW(J),SIG 

C . VARIATION  OF  CHARGF  SIZF . 

DO  1  Of,  L  =1  *  S 

C . VARIATION  OF  SCALFD  DISTANCF . 

DO  107  L  1  =  1  *50 

c . blast  paramftfrs  for  scalfd  distance 

P  I=P0(L1  1 

PR  =  ?.*PI*( 102.9+9.*P1  )/( 102.9+PI 1 
DUR  =  T  All  (  L  1  )*CS(L)**A 
I F ( Ll-37)  117*117.118 

117  S  =L 1 

R= ( S  +  3 .  )  *C  S  t  L ) **A 
GO  TO  119 

118  S  =  L 1 

R  = (  (S-37. )*5.^90.)*CS(L) **A 

119  Q0=2. 5*CPI**2. 1 /( 102.9+PI 1 

VFL  =  1 187.*SQRT ( 1 •  +  (  6  •  *P  I  )/(  102.9+PI )  ) 

C . TIMF  TO  STAG.  AND  MAX.  AVFR.  PRFSSURF  ON  ROOF 

TA=( 1 5 . 7 / ( 2 . *VF  L )  1*1000. 

TS=?.*TA 

PTS=PR*( 1 ,-T A/DUR )*FXP (-TA/DUR ) 5*00* ( 1 . - T A / DUR ) **? . *FX P ( -? . *T A 
* / DUR ) 

IE(CP-DUR)  120,120,121 

C . CASE  I.  CRITICAL  PERIOD  LFSS  THAN  DURATION . 

120  T IME=TS+.5*(CP-TS) 

T IMS=TIMF-TA 

P  T I =P I  * ( 1 ,-T I  MS /DUR ) *FXP ( -T I MS/DUR ) 3*00*  (  1  .-T I  MS /DUR 1**2. 
**FXP( —2. *T I MS/DUR ) 

CAL  I=.5*PTS*TS+.5* ( PTS-PTl  )*(TIMF-TS)+PTI*(TIMF-TS)  +  .5*PTI 
** ( CP~T I  ME ) 

GO  TO  122 

C . CASE  II.  CRITICAL  PERIOD  GREATER  THAN  DURATION . 

121  TIME=TS+.5*(DUR-TS) 

T IMS=TIMF-TA 

PT I =P I* ( 1 ,-T IMS/DUP )*FXP ( -T I  MS /DUR ) - . 3*00* ( 1 ,-T I  MS /DUR ) **2. 
**FXP( -2.*TIMS/DUR  ) 

CALI=.5*PTS*TS+.S*(PTS-PTI)*(TIMF-TS)+dTI*(TIMF-TS)+.S*PTI 
** (DUR-TIMF  ) 
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n  n  n  n  n  n  non 


12?  IF(CI-CALI)  107,126*126 

C . MINIMUM  DISTANCF  FOR  G1VFN  CHARGE  SFLECTFD . 

126  PRO  =P  TS 

WRITFf  3,8)  CS(L),R 
WRITFt 1, 17)  PI*  PR ♦ OUR 

....CRITICAL  IMP  AND  CALCULATED  IMP.  FOR  HOUSF  FPAMF  ARE 

• • • • COMPARFD  FOR  THIS  SFLFCTFD  CHARGF  AND  DISTANCF . 

WRITFf  3, 1 3 1 

....GLASS  FRAGMENT  INJURY  MODFL . 

WR I TE ( 7,53) 

I F ( P I -S . )  170,171,171 

171  WRITEf  7, SO)  PI 
GO  TO  17? 

170  FW=-.7B*PI+7.Q 
FV=18.7*PI+6?.5 
CV  =  ?30. *FX  P ( - 1 .  9*FW 1+103. 

IF(CV-FV)  173,174,174 
174  WRITEf  3,51  )  PI 
GO  TO  17? 

173  M  = (  (FV-CV) /50. ) 

S=M*.1+. 1 
WRITFf  3,5?)  PI  ,S 
17?  DO  138  11  =  1  ,K? 

DO  138  I2=1*K2 
DO  135  13=1,  K? 

....BLAST  LOADING  ON  OUTSIDE  OF  FRAMF  SFCTION . 

TS= ( 1 0. R /VFL ) *1 000. 

DD  =  TC?(  I  1  ,  I?,  13  ) 

PTS  =P I  * (  1 .-TS/DUR 1  *FXP ( -TS/DUR ) +00* ( 1 .-TS/DUR )**7.*FXP(-?. 
**TS/DUR) 

IF(DP-TS)  136,136,137 

• . . • CAS F  I.  CRITICAL  pfriod  lfss  than  time  to  stag . 

136  CALF=. 8* ( PR-PTS )*TS+TS*PTS 
GO  TO  140 

137  IF(DP-DUR)  138,138,139 

....CASF  II.  C  PFRIOD  LESS  THAN  DURATION . 

138  PCP=P  I* ( 1  .— DP/DUR 1 *FXP ( -DP/DUR ) +QO* (l.-DP/DUR)**7.*fXP(-?. 
**  DP /DUR ) 

CAL F=. 5* (PR-PTS )*TS+TS*PTS+.5* (PTS-PCP) * ( DP- TS ) +PCP* ( DP-T S ) 
GO  TO  140 

139  CALF=. 6* (PR-PTS ) *T S+ TS*P TS+ . 5*PT S* (DUR-TS) 

....BLAST  LOADING  ON  RACK  OF  FRAME  WALL . 

140  PTF  =P  I  * (  1 .-17. /DUR )*FXP( -17./DUR ) 

1 F ( DP-DUR )  148,148,149 

148  PCP=P 1*( 1 .-DP/DUR )*FXP ( -DP/DUR ) 

CALR=.5*PTF*17.+.S* ( PTF-PCP )*(DP-17.1+(D P-17.1* POP 
GO  TO  150 

149  CALB=.5*PTF*17.+.5*PTF*(DP-17. 1 

150  CAL I=CALF-CALB 

135  WPI TEf  3* 14)  UI? ( 1 1 , 12,13 ), CALI 

....CRITICAL  IMPIJLSF  AND  CALCULATED  IMP.  FOR  HOUSE  BRICK  WALL 

....ARE  COMPARED  FOP  SFLECTFD  CHARGE  AND  DISTANCF . 

WPITFf  3,15  1 
I  6=1 

DO  145  14=1. K3 
DO  145  I  5= 1 ♦ K3 

C . CRITICAL  IMPULSF  OF  BRICK  WALL  IS  FUNCTION  OF  BLAST 
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C . LOADING  ON  ROOF 

c=hoF(  I  4  )  *  1  OOOOOO  . 

T'-'  =  S  . 

TFT  =  HnT (IF) 

TS  =  (7.*7./VFL)*lOOO. 

DP  =  T  C  7  (  14*  IS  *  16  ) 

PV=1S.Q+.5*PR0*?R0.4 
01  =  S7.* ( TST*fi.+PV )  /Pin. 

Yl=  (S. *00**7. *Qi  ) / ( 7R4.*F*4? .7 ) 

07  = ( 4.* ( R.-Y1 )  * (7 .*PV+S0. ) ) /oo. 

r o?=( 4.* ( tw-7 .)*(?. *Pv+Fn.))/PO. 

WORK=(  s.*Yl  *01  )  +  (  .S*(  7.-V  1  )  *(0?-rQ7  )  )  +  C07*  (  7.-Y1  1 
V SO  =  WORK /  •  04  7 1 
V  =  SOR  T  (  V  SO  ) 

z i = ( ?s./G)*v+.oi?*v 
U  I  3  <  14*  I  S*  I  6  )  =  (ZI*1000./4S.  )*.*> 

IF(DIJR-TS)  141*141  *14? 

C . C ASF  I.  DURATION  OF  BLAST  LFSS  THAN  TI^F  TO  STA<S. 

141  CA(.  I  =  .  s*PR*DUR 
GO  TO  1 4  S 

C . CASF  II.TIMF  TO  STAG.  LFSS  THAN  DURATION 

147  PTSsPI*( 1 ,-TS/DUR)*EXP(-TS/DUR) +00* ( 1 ,-TS/DUR)**?.*FXP ( -7 . 

**TS /OUR ) 

C . CALGULATFS  point  OF  MINIMUM  effectivf  pressure . 

D I V= ( DUR-TS ) /SO. 

DO  14  7  I  7= 1  ,  SO 
S  =  I  7 

TFMD=TS+S*DTV 

PrNn  =  PI*(l.-TFND/DUR)*FXP(-TFND/DUR)-t-0O*(l.-TEND/DMR)**7.*FXP(-?. 

*  *  T  FND/ DUR ) 

I  F ( PEN  D-  •  7 )  144*144*147 

147  CONTINUE 

144  CAL  I  =  .S#TS*tPR-PTS) +TS»PTS+.S*PTS* ( TFND-TS) 

1  4  S  WRJTFt  7*14)  U  I  7  (  I  4  *  I  s  ,  I  f,  )  ,  f  A  L  I 
GO  TO  1  Oft 
107  CONTINUF 
loft  CONTINUF 
10S  CONTINUF 

CALL  LINK(CALLR) 

7  FORMAT ( 1 X , ' CR I T I  CAL  PFRIOD  AND  IMPULSE  FOR  HOUSF  ROOF',//) 

4  FORMAT ( 1 X, *F • *RX , 'WEIGHT • *RX , *M. OF  RUPT'*PX,'C  PFRIOD'*RX* 

*  •  r  I  mpi  il  SF  *  *  ) 

s  FORMAT ( FQ.O ,4F1 7. 7 ) 

7  FORMAT ( 1 X* • mini  mum  DIST.  FOR  F = , , FR . 0 » 1 X  ,  ' WF I GHT = '  ,  F 1 0  .  ? *  1  X  * 
*lX,'M.OF  RUPT  =  '  ,F10.7  ,/ ) 

R  FORMAT ( 1 X CHARGF  WE  I  GHT  =  ' , F R . 0 , 1  OX , ' D I  ST A NC F - ' , f 1 o . 2 , / ) 

9  FORMAT ( 1 X  CR IT ICAL  PFRIOD  AND  IMPULSE  FOR  H  FRAMF  SECTION**//) 

10  FORMAT ( 1 X* ' CR I T ICAL  PERIOD  FOR  HOUSF  RRICK  WALl',//) 

11  FORMAT ( 1 X *' F '  *1  OX *• TENSI LF  STRENGTH 1  OX *' C  PFRIOD',  ) 

17  FORMAT(F9.0,7F1 0.?) 

17  FORMAT ( 70X »' HOUSF  SRAMF  SFCTION',/) 

14  FORMAT  (  1  SX  *•  RFOIJI  RFD  I  MPUL  SF  =  •  ,  F  1  0  .  7  *  ?X  ,  '  C  A  L  C  .  I  mpi  IL  SF  =  •  ,  F  1  0  .  ?  *  ) 

15  FORMATt ?0X , 'HOUSF  RRICK  WALL',/) 

17  FORMATt SX, '  INC.  OVDR FS . =  ' * F 1 0 . 2  * • PF A K  RF FL . P RFS .  =  • , F 1 0 . ? * 

* • DUR • = ' , FI  0.7 , /  ) 

18  FORMAT ( 1 X *' LBS/SO  I N ' *SX ,' LBS/FT •  *SX *• LRS/SO  IN ' *SX * 'MSEC  '  * 

*RX , • PSI -MS ' , / ) 
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1R  FORMAT ( 1 Xi •  |  RS/SO  I N ' » 1  ox » • LRS/SQ  I N  •  » 1  OX » • MSFG • » / ) 

50  FORMAT t 1 X ,' FRAGMFNT  SIZES  AT  P I  =  *  if  1 0.7 » ' ARF  TOO  SMALL  To  CAUSF' 
*1X, • SFRIOMS  INJURY',/) 

51  FORMATt 1 X, • FOR  P I = • , F 1 0 . 7 » 7X , • PROB .  OF  SERIOUS  WOUMDS=0 . ' * / ) 

57  format t i x, ' for  di=* ,fio.7»7X» 'prob.of  serious  wounos=» *fio.2*/) 

S3  FORMATdX,  'PROBABILITY  OF  GLASS  FRAGMFNT  I  N  JURYZ  '  ,  ) 

FND 

//  m IP 
♦STORF 
//  JOB 
//  DIJP 
♦DFLFTF 
//  FOR 

*ONF  WORD  INTFGERS 
♦ARIThMETIG  TRAGF 
*TRANCFFR  TRAGF 
* IOCSI GARD, 11 3?  PRINT FR ) 

COMMON  PO  (  sn  )  ,T  All  (  sn  )  ,CPR  (  sn  )  ,HRF  (  A  )  ,HRW(  A  )  ,HRS  (  A  )  ,HFF  (  4  )  , 

*  HFW ( A )  , HFS  t  A )  , HBF  t  A )  ,H3T (A)  »CRF  C  4 ) * GRW ( A )  ,GRS(A)  ,CDF(A)  ,CDW(A)  , 

*  CDS ( A ) *SRWE  t  A)  ,  SRT (  4  1  '  SRF ( A )  ,  SRW ( A )  ,  SR  S ( A )  »OHF ( A )  *OHT ( A )  'OFF ( A )  » 
*CTG( 4 )  ,  TRW ( A )  »CPW ( A )  '  TRCG ( A )  , TRWT ( A )  »DORS ( A )  ,  WM AN t  A )  »  PM AN ( A )  » 

*  OFT ( A) » !A( 11 ) ♦< 'NOT'UI ( A  'A  »A)  »TC( A'A'A) *CS( S )  ,RCG<  A ) 'BW ( A  )  »CCG  t  A  )  » 
*AIPF(A)»AIRF<  l0)»TC7<4'4'4)  »TC3<4»4'4)  »  U I  7  <  4  ♦  4  »  A  )  »  U I  3  (  A  ♦  A  *  A  ) 


WS  l)A  HOI  ISF 


GHRCH 


PART  IV.  GHPGh  SUPROUTINF 


C.  . 

C 

C . 

C . FIXFD  PARAMFTFR  VALUES . 

A  =  1 ./o. 

FFFT  =  4S8 . 

A  I =?SAS. 

C=7 .RS 
G=3R6. 

C . DFTFRMINFS  mo.  OF  GHURGH  ROOF  PARAMTTFRF 

DO  7"0  1=1 » A 

I F ( CR F (  I  )  )  ?n  i  ,7^1  ,700 
701  K 1 =  I—  1 

GO  TO  273 
700  CONTINUE 
<1=4 

C . DFTFRMINFS  MO.  OF  GHURGH  ROOF  DFGXING  PARAMFTFRS. 


770  DO  ?i P  I =1 » A 

IFtGDFt  T  )  )  7 1 R  > 

7 1 R  <7=1  —  1 

GO  TO  707 
2  1  R  CONTINUF 
K  7  =  A 


71 R, 71  R 


C. 

C. 

C. 


.CRITIGAL  PFRIOD  AND  IMPUISF  FOR  GHURGH  ROOF, 


707  DO  703  1=1 1 

E=CRF ( I ) * ] Oooooo , 

DO  703  I  1  =  1  ,<1 
W  I  =CRW (  I  I )*37. 

W  =W  I  /  AS  R  . 

DO  703  111=1, K1 

s  i g  =  gr s (  ti  I  )*pnnn. 

PP  =  R. R 7* SORT (  t  E* A  I *G ) / (W*FFFT**4, 


)  ) 
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P=PP/4.96 

TC< I ♦ 1 1 ♦ I T I ) =ionn./ (4.*P) 

W0= (SIG*8.*AI >/ <C*FFET**?. ) 

U=  (  ( WO**? .  )  *  <FFFT**4.))/(240.*F*4I) 

WT=W*FFFT+140. 

H  =  SOR  T  <  U*?.*WT /G  ) 

?0  4  UK  I  »  T  I  *  I  I  I  )  =  (H*1  O-in,  )/fiu?UO. 

W  R  I T  p  (  4,6) 

WRITE! 4,4) 

WR  I  TF ( 4, 18 ) 

DO  PDA  1=1  ,K1 
DO  ?04  J=1 ,<1 
On  ?n u  J1=1*K1 
F  =r rf  < i ) * i onnono. 

S  IG=CRS( J1  ) *  1 ooon . 

?oa  wpitf(^,a)  f  »crw  (  j  ) ,  mg  ,tc  (  i  ,  j*  ji  )  »ui  (  t  ,  j  »  j  n 
c . critical  pfriod  and  impulse  for  church  decking . 

DO  ??0  I =1 »  K  ? 

E=CDF ( I )*lOOOOon. 

DO  ??0  J  =  1  ,  K  2 
W I =CDW( J ) *  1 4 . 

W=W  1/180. 
on  ??0  I  I  I = 1  ,  K  ? 

s  ig  =  cos  (iin*  i  pppo. 

PP  =  Q.  8  7*  SORT  (  ( F*46.64*G  > / <W*18  0. **4 . )) 

P=PP/4.R6 

TC?tI,J,ni)  =  lO00./(4.*P) 

W0=SIG*8 .*46.64/ ( 1.8*180.**?.) 
U=((WO**?.)*(18O.**4.))/(?40.*F*46.64) 

WT=w* 180. 

W-SORT  <  I  I*?  .  *WT/G  ) 

??0  U I  4  (  J  , J,  J I  I  ) =M* 1  POO. /?1  60. 

W  P  I  T  F  (  4 , 0  ) 

W R I Tf ( 4 , 4  ) 

WRI TF<  4 , 1 8 ) 

DO  ??1  I =1  ,  K  ? 

OO  ??1  J  =  1  ,(/? 

DO  ??1  J ]  =  1  , K  ? 

F=ro>r(  i  )*ioopppp. 

4  IG=COS ( j i ) *  1  POPP. 

??1  WR1TF(T»4)  F  ,COW( J)  , S IG , TO? (  I  , J ,  111  ,U  I  ? (  T  ,  J  ,  J  1  ) 

C . 

C . DFTFRMINES  MINIMUM  STAND-OFF  STORAGF  DTSTANCF  FOR  FACH 

C . CHAPGF  SI7E  AN|0  SFT  OF  PARAMETER  VALUES . 

C . 

C . VARIATION  of  PADAMFTFR  VALI'FS . 

no  ? o f  j  =i  ,i/i 
no  ? n 4  j  =  i  ,ri 
DO  ? o  4  j  i  =  i  , < i 
r  T  =  u  i  (  i,  j,  ji  ) 

C P  =  TC  (  I  ,  J, Jl  ) 

F  =  CRF ( I  )  *1 Oppppp. 

SIG=CRS( Jl  )*1P0P0. 

WRJTF(  4,7)  F  ,GRW(  J  )  ,S  jr, 

C . VARIATION  OF  CHARGF  S  I  7  F . 

Do  ?  o  6  1=1,4 

C . VARIATION  OF  SCALFD  DISTANCE . 
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DO  ?n 7  L 1  =  1  ,R9 

C . RLAFT  PARAMFTFRF  FOR  FCALED  DISTANCF 

P  T=PO(L 1  ) 

PR=7.*pI*(10?.o+a.*PI)/(1  0?  .  9  +  P  I  ) 

OllR  =T  At  I  (  LI  )*GF(  L)**A 
T  r ( l 1-87  )  ?08 , 7 OB ,  70o 

708  F=L1 

R  = (  F  +  8 .  I *GS ( L ) *  *  A 
GO  TO  710 

709  S  =  L1 

R  =  (  (S-87.)*8.+40.)*CF(L)**A 

710  00=7. F*( PI**?.  ) / (  1 97.0  +  P  I  ) 

vfl  =  i i n 7. * sort ( i .  +  < s  .*pi i/i  107.0+pi ) ) 

T A= ( 18.78/ (?.*VFL I 1*1 000. 

T  F  =  ?  .*T  A 

PTS=PR*  t  1  .-Ta/OUR  )*FXP  (-Ta/DUR  ) -. 7*00*(  1.-Ta/DUR)**?.*fXp(-2. 

*  *  T  A  /Ol  IR  ) 

IF(GP-DUR)  211.711.21? 

C . CAFF  I.  CRITTGAL  PFRIOO  LFFS  THAN]  DURATION . 

711  T  TmF  =  TRj-.8*(  GP-TF  ) 

T imf=T IMF-T  A 

PTI=PI*(  1  .-TTMf/OUP)*FXP(-TiNF/OUR)-.8*00*(  1  .  -  T  I M  S /OUR  )  *  *  ?  . 
**FXP(  -7  .  *T  I  MF/OtlR  I 

CAL  I  =  .F*PTS*TS+.S* ( PTF-PT I ) * ( T I  ME- T S ) +PT I  * ( T I MF-T S )  +  . S *PT I 
**(GP-T I^F ) 

GO  TO  718 

C . GAFF  II.  GRITIGAI  PFRIOD  GRFATFR  THAN  OIJRATION . 

717  T  IMF=TF+  .e*  (  ni  IR  — TF  ) 

T  I M  F=  T JME-T  A 

PT  I  =PT*  <  1  ,-T  IMS/DIJR  1  *FXP(  -TTMS/OUR)  -.  8*00*  (1  .  -  T  I  MS  /DUR  )  *  *  ?  . 
**FXP(  -?.  *T  I  MF/Ol  IP  ) 

GALI=.f*PTS*TS+.8*  (PTF-PT  l I  * ( T I MF-T S I +P T I  * ( T I MF-T  R  )  +  . R  *PT  I 

*  *  (  01  IP  — T  I  MF  I 

718  IF(GI-GALI)  707, 71  7, 717 

C . MINIMUM  O I  ST  ANOF  FOR  GTVFN  GHARGF  SFLFOTFO . 

’P  MRITFI  ■>  ,8  I  GF  (  L  I  «R 

WRITF(8,17)  PI.PP.OliP 

G . 

C . GRITICAL  IMP  AND  GALGULATFD  IMP.  FOR  OHUPOH  DECKING 

G . ARE  COMPARED  FOR  THIS  SFLECTFD  CHARGE  AND  DISTANCE . 

G . 

C . GLASS  FRAGMENT  INJURY  MOOFL . 

WRITFI 8,88) 

T  F ( P 1  — F . )  770,771  .771 

?  7 1  WRITE (7,80)  PI 
GO  TO  77? 

770  FW=-.78*PI+8.9 
FV=1 fl.7*P I +  67.  F 
CV=?80.*FXP ( -1 . 9*FW ) +108 . 

IF(CV-FV)  278, 774, 774 
774  WP ITFI8.Fl)  PI 
GO  TO  77? 

778  M  = (  ( FV-GV I / F  o .  , 

F  =  M*  . 1  +  .  1 
WR I TF(  8, 87  )  P  I  ,  F 
2  7?  WR  I  TF  (  8 , 1  0  ) 

DO  72?  I  1  =  1 ,K2 
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no  ???  i  ?=  i 
r>n  7  ??  I  8  =  1  ,  v  ? 

no-jr ?<  n  »pw^) 

TS=<1  #  ^  /  \/F  L  )  *  I  °  nn  • 

RT  S  =  PI  *  (  1  .-TS/DUR  )  *FXP  ( -TS/DUR  )  8*00* {  1  .  -T  S/DUP  )  *  *?  . 

**rXP<  -7#*TS/OlJR  ) 

TF(np-nuR)  ??a,??4,??s 

C . CASF  I.  CRITICAL  pfriod  lfss  than  duration . 

7  7  4  T  I  M  F=  T  S+  #  ^  *  (DP-TS) 

T  I  M  <;=  T  IMF 

DT  I  =D  I*  (  1  .-T  I  v.Q  /OUR  )  *FXP  (-T  I  MS  /DIJR  1  -  .  7*00*  (  1  .  —  T  I  ms /OUR  )  *  *  ”«  . 

**FXP { *T I MS/HHR ) 

C  A  L  I  =  .  8  *  ( PR-PT  S  )*  js+TS*PTS+  .  8*  < PTS-PT I ) *  <  T I MF-TS ) +DT I  * ( T I MF 

*  -TS)+.8*PTI *  ( DP-T  TMF) 

GO  TO  ??? 

C . CASR  II.  CRITICAL  D FR I OD  GRFATFR  THAN  DURATION . 

7?8  j  TMF=TS+.**  <  DUP-TS ) 

T  T  mf=T I mf 

DT  I  =P  T  *  {  1  T  IMC/OiiP  )  *CXP  (— T  I  Mr,  /  nitP  )  —  .  7*QO*  {  1  .  -  T  IMF  /  OUR  )  *«  ?  . 

*  *FXP  {  *T  I  ms/OI  IP  ) 

CALI=.^*<PR-PTS)*TS+T8*dTS+.8*  <  PTS-PT  J  ) *  <  T IMF-TS ) +PT I  *  <  T IMF 
*-TS  )  +  .  8*PTI*(DUR-T  IMF  ) 

???  WR  I  TF  (  7 , 11  )  Ul?m»I?»l^)»CALI 
GO  TO  ?06 

707  CONTINUF 
70^  C ON  T  I  Ml  IF 

708  rnMTINPF 

CAL  !  !  INMCAl  [_o  ) 

4  FORMAT  < 1 X  *  • F  •  *  KX  *  • WEIGHT •  ♦ RX  * • M.OF  RUPT  •  ♦  R  X  *  »  C  PF  R  I  OD  '  ♦  R  X  ♦ 

*  • C  I mphl  SF  *  *  ) 

8  FORMAT (FO.0.4F1  7.7) 

A  FORMAT (1  X  CR  IT  ICAL  PFRIOD  AND  IMPULSF  FOR  CHURCH  ROOF •  ♦ / / ) 

7  FOPMa  T ( 1  X ♦ 1  MINIMUM  D I  FT  #  FOR  F- •  * F R # 0  » 1 X ♦ • WF } GHT =  *  * F 1 0 . ?  *  1 X  « 

*  1  X  *  * M. OF  P  U  P T  =  *  *F1 /) 

R  FOPMA  T ( 1 X ♦ • CHAPGF  WF  T GHT  =  • ♦ F R . 0 ♦ 1  OX ♦ • D I S T AMCF= • *  F 1 0 . ? 9 / ) 

O  FORMAT  (  1  X  9  •  CR  I  T  IC  Al  nrqjor)  AND  IMPULSE  ^OR  C  nrc^IMG',//) 

I  0  FORMA  T ( 7 OX  < • CHURCH  DFfK IMG • 9 /  ) 

II  FORMAT  <1  SX  9  •  PFOIJIRFD  IMPULSF  =»»FlO*?*?Xi  •  C  ALC  T  MPU|  SF=»*FI0.2O 

17  FORMAT <8X9 •  INC.  OVpRFS.=#  »  F I  0 • 2  ♦  *  D  F  A  K  RF F L • PRFS .  =  •  * F 1 0 . ? ♦ 

*  tniJP.=  '  »FlO.?*  /  ) 

18  FORMAT (  l X ♦  • L  RS/ SO  I N  *  9*>X 9  ♦  l  RS / F  T  »  *6  X 9 • L  RS /So  I N • 9 A X 9 • M5F C  * 9 

*  AX  *  »D^  T-mf •  , / ) 

80  FORMAT ( 1 X FRAGvFMT  SlZFS  AT  P I = t 9 F 1 0 . 7 9 « A RF  TOO  SMALL  TO  CAUSE1* 

*  1 X  9  • SF  PIOUS  INJURY  *  9 / ) 

81  FORMAT ( 1 X 9' FOR  PI=t ,F10.?*?X* »PROR.  OF  SFRIOUS  W0UNDS=0. • * / ) 

8  2  FORMAT  (  ]  X  FOR  P  I  =  •  *F10.?*?X*  *  PROR  OF  SFRIOUS  WOUNDS  =  *  *F,10.?*/) 

87  FORMAT ( 1 X ♦♦ PRORARIL ITY  OF  GLASS  FRAGMFNT  INJURYZ**) 

FND 

//  D(fP 

*  S T OR p  MS  UA  CHRCH 

/  /  JOR 

//  Dud 

*DFLFTF  SCHOL 

/  /  FOR 

*ONF  UORD  INTEGFRS 

*  AR I T  HMFT  I C  TRACF 
*TRANSFrp  TRACF 
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o  n  o  n  n  n  noon 


*  IPrS ( CARD*11??  PRIMTPR) 

Common  PO(SO) »TaU<SO) »CPRC  SO) »HRF<  4) »HRW( 4) *HRS( 4 ) *HFE<  4) * 

*HFW ( 4 ) *  HFS ( 4 )  *  HRE ( 4 ) *  HRT ( 4 )  *CRE<4 ) »CRW<4 )»CRS(4)»CDE(4)  »CDW(4 )  » 
*CRS ( 4 ) *SrWE  <  4 )  *  SrT (4)  ♦  SRF  <  4 )  *  SRW ( 4 )  *  SR  S ( 4 )  *OHF  <  4 )  *OHT ( 4  )  *OFE<  4)  * 
*CTG( 4) *TRW ( 4 )  *CPW ( 4 )  »TRCG<4)  *  T  R  W  T (4 ) »DORS<  4 )  »WMAN ( 4 )  *  PMAN (4)  » 

*OFT (4)*IA<11)*K *  NOT  »  U  I (4>4*4)>Tc<4*4>4)*CS(S) »RCG<  4)  * DW ( 4 )  »CCG ( 4  )  » 
*AIRF<4)*AIRF<l0)*T0?<4»4»4)»TC3(4»4»4)»UI?<4*4*4)*UI3<4*4»4) 


PART  V.  SCHOL  St  IRROUT  T  N  F 


....FIXFO  DARAMf-JFR  VALUFS . 

A  =  1  ./^. 

FFFT=1?6. 

TW  =  8. 

A  I  =4?  •  7 
G  =  Tf>6. 

I  I  I  =1 

. . . .DFTFRM JNFS  NO.  OF  SCHOOL  R.  WALL  PAP  AMFTFRS 
DO  800  1=1,4 

I F ( srw F (  I )  )  sni , so i  ,snn 

SOI  K 1 =  1  —  1 

GO  TO  3?  9 
SOD  CONTINtlF 
<1=4 

. . . .OFTFRMTNFS  NO.  OF  SCHOOL  ROOF  P A R AH FT FR S . . . 
no  s?n  i=i,4 
IF ( SRF ( I ) )  SSI ,1S! , ssn 
SSI  <?  = I- 1 

GO  TO  SO? 

SSO  CONTINUF 
<7  =  4 

....CRITICAL  PFRIOD  FOR  SCHOOL  RRICK  WALL . 

on?  DO  sns  I  =1  ,<  1 

F  =  sRWF (  T  )*iOOOOnn. 

DO  SOS  11  =  1, <1 
T  S  =  S  R  T  (II) 

W  =  .  S  S 

PP=q. 87* SORT ( (  F*AI *G) / ( W*FEFT**4. ) ) 

P=PP/S .96 

°OS  TC<  I  ,  1 1 ,  1 1  I  )  =  1  or,o.  /  (  4.*p  ) 

WR I TF ( s  *  9) 

WRITF( s, 14) 

WR  I  TF ( S, S4 ) 

DO  S04  1=1  *  < 1 

DO  SO  4  J  =  1  , K 1 
I  I  1=1 

f=srwf< i )*innonon. 

s04  WRITF(s,ll)  F  *  SRT ( J ) ♦ TC  <  I »  J» 1 1  I ) 


CRITICAL  PERIOD  AND  IMPULSE  FOR  SCHOOL  ROOF 


DO  S3?  I =1 *<2 
E  =sre ( I )*ionooon. 
DO  33?  J=1*K? 

W I =  S R W ( J ) *?  8 . 
w=w I /SS6  • 

DO  33?  I  I  1  =  1  ,  K? 
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6IG=6R6(  n  i  )  *1  ooon  . 

PP  =  9.  8  7*  SORT  (  (F*817.A1*G)/(W*??6**A.  )  ) 

P=PP/S .96 

TC?< I » J » 1 1  I >=1ooo./(a.*pi 

W0= ( S IG*8.*M7.A1 ) / ( 6 .6*7?6.**?. ) 

11= ( ( VO**?. )*(776.**6.) 1  / ( ?AO.*E*61 7. 1 A1 
WTsW*11fi  . 

H=60RT(  I  !*?  .  *WT  /G  ) 

77?  UI? (  I , J, I  I  I ) =H*1000./806A. 

W R I TF( 7.17) 

WR I TE ( 7,16) 

VRITFI 7,65 ) 

DO  777  I =i  ,K  ? 

DO  77?  J=1 ,K? 

DO  777  i  I  1  =  1  ,K? 

F  =  6RF ( I )*loooooo. 

6IG=6R6(  I  I  I  )*10O00. 

777  WRITE! 7.17)  E*SRWIJ1»SIG*TC?(I»J»III)»UI?« I  *  J  *  I  I  I  ) 

c . 

C . DETERMINES  minimum  STAND-OFF  STORAGE  DISTANCE  FOR  PACH 

C . CHARGE  S  1 7  F  AND  SFT  OP  PARAMFTFR  VALUP6 . 

c . 

. . VARIATION  OF  PARAMPTFR  VALUPS . 

DO  706  1=1 
DO  70S  J=1 ,K 1 
CP=TC( I ♦ J»  I  I  I  1 
F  =  SBWF (  I  1*1000000 . 

WRI TF(  7  ,  1?  )  E  »  SRT ( J ) 

C . VARIATION  OF  CHARGE  SIZF . 

DO  706  I  =1 ,6 

C . VARfAT.ON  OF  SCALPP  DISTANCF . 

DO  707  L  1=1,60 

C . RL  A  6  T  PARAMPTPRS  FOR  SCALPD  DISTANCF 

P I =P0( L  1  ) 

Dl  IR  =T  AIM  L  1  )  *  C  c'  (  I  '  *  *  A 
I F (  L 1  - 7  7  )  708,70  6,709 

7  0  8  S  =  L  1 

P  = (  6  +  7 .  ) *C S ( L ) ** A 
GO  TO  710 

709  6  =L  1 

r=(  (S-77.)*6.+AO.)*CML)**A 

710  00  =  ? .9* (  PI **? .  ) /  (  1 0?  ,9  +  P I  ) 

VFl.  =  1  18 7.* SORT!  1 ,+(6.*PI  )/(  10?.9  +  PI  1  ) 

T  A= ( ?8. / ( ?.*VFL 11*1000. 

PTA=PI*(1.-TA/D'!R)*PXP(-TA/DI1R  ) -. A*00*(  1  .-T  A /DUR  )  **?  .  *F  X  P  (  -?  . 

**T  A  /Dl  IR  ) 

C . CPITICA!  IMPULSF  OP  SCHOOL  R  WALL  18  FUNCTION  OF  LOADING  ON  ROOF 

PV=7S.-L 1*1 . 

T  6T  =  8RT ( J ) 

01=(A.*TW)*(TST*TV+PV)/(?.*FEFT) 

Y 1  =  ( S.*FFFT**?.*01  1 / ( ISA  . *F *A  I  1 
07= ( 4 . * ( TW-Y1 1  * ( ?.*PV+7  0. )  ) /FFFT 
C0?= ( A. * ( TW-7 . ) *( 7. *PV+70 . ) ) /FFFT 

W  OR  K  =  (  . 6  * Y 1 *01  )  +  (.6*(7.-Yl  )*(0?-C0?)  ) +T0?* (  7 . -Y 1  ) 

V 60= WORK / . 06 A 
V  =  SOR  T  (  V  60  ) 

T I = ( 76 . /G  )*V+.0?8*V 
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n  n 


CI=TI*I0O0./63. 

pr  =  ?.*pi*( io?.9+4.*pi  )/( m?.9+Pl  ) 

t  s  =  ( ^ *  i  o.  s/v^l  i*nno. 

I  F  (  HIJR-T  9  )  37S,376,376 

c . c  a  s  f  i.  duration  of  blast  lfss  than  t imf  m  stag. 

3?S  CAI.  I  =  .  A*PR*niJR 
GO  TO  si  A 

C . CASE  II.TIMF  TO  STAG.  LFSS  THAN  DURATION 

376  DIV=(DUR-TS) /SO. 

C . CALCULATFS  POINT  OF  MINIMUM  EFFECTIVE  PRESSURE . 

DO  377  17=1, SO 
S  =  I  7 

T  fmd=TS  +  S*OI V 

P  fND  =  P I  * ( I  .-TFND/DUR ) *FXP( -TFND/DUR ) +  QO* (  1  . - T FND/DUR ) **7.*EXP ( -?. 
**TFND/DUR ) 

IF(PFND-.3  )  328,328*  377 

327  CONTINUE 

328  P  TS  =  P I  * (  1 ,-TS/DUR )*EXP (-TS/DUR )+QO* ( 1 . - TS /DUR ) ** 7 . *F XP ( -? . 

*  *TS / DUP ) 

I F ( TS-ro  )  311 ,3  n  *31? 

C . C.  PFRIOD  LFSS  THAN  TIMF  TO  STAG . 

312  A\/FR=(PTS  +  PR)  /7. 

IF(AVFR-.S)  310,313,313 

313  CAL  I  =  .S*< PR-PTS )*TS  +  TS*PTS+.S*PTS* ( TFND-TS ) 

GO  TO  31a 

C . C.  PFRIOD  GREATFR  THAN  TIMF  TO  STAG . 

311  I F ( CP-TFMD )  330,320,371 

370  P CP =P I* ( 1 . -CP/ DtlR  )*FXP  < -CP /DUR ) +00* (l.-CP/DUR)**o.*FXP(-7. 

*#co  /  m  >p  i 
AVFR=(PR+PCP)/7. 

IF(AVFR-.S)  319,316,316 

316  CAL  I =. s* (PR-PTS ) *TS+TS*PTS+PCP* ( CP-TS 1 + . S* ( P TS-PCP ) * ( CP-TS ) 
*+.S*PCP* (TFND-CP ) 

GO  TO  3 1  A 
3?l  A  VF  R  =  PR / 7 . 

IF(AVFR-.S)  319,373,373 

373  CAL  I*.  **  (PR-PTS )*TS  +  TS*PTS+.s*PTS*( TFNO-TS  ) 

3 1  A  IF(CI-CALT)  307,31  9,31  0 

C . MINIMUM  DISTANCE  FOR  GIVFN  CHARGE  SFLECTFD . 

319  WRITF(3,1S)  CS!L),R,CI 
WRITE! 3,47)  PI, PR, DUR 


CRITICAL  IMPULSF  and  C  ALCUL a  Tfd  IMPULSF  FOR  SCHOOL  POOF 
ApF  COMPARFD  FOR  THIS  SELECTED  CHARGF  AND  DISTANCE . 


C . GLASS  FRAGMFNT  INJURY  MOOFL 

WRITE! 3, S3) 

I F ( P I  — S . )  370,371,371 

371  WRITE!  3, SO)  PI 
GO  TO  377 

370  FW=-.7**PI+3.9 
FV=18.7*PI+6?.S 
CV=730.*FXP ( -1 . o*fw ) +1 03 . 
IF(CV-FV)  373,374,374 
374  WRITE! 3, SI )  PI 
GO  TO  377 

373  M  = (  (FV-CV) /50. ) 
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-S=M*.  1  +.  1 
WRITFf  i,  57)  PI  ,  S 
177  VPITF(I.IO) 

in  I  1  =  1 ,k? 

on  114  17=1, K7 
DO  114  1 1  =  1  »K2 
DP=TC? <  I  1  , I?  ♦  II  ) 

T  A  = ( ?fi . /VFL  >  #  1  OOP . 

r . T  I m F  and  max.  AVER.  ORFSSURF  on  school  roof . 

T  Z  =  (  78. / ( 7. *VFL )  )* IDOO. 

PTZ=PI*( 1 .-TZ/ntiP)*FXP<-TZ/DUR)-.i*QO*( 1 TZ /DUP 

*  *T7  /Di  |R  ) 

IF(DHR-TA)  115,115,116 

115  CAL  I  =  .  5*  (PTZ*DIIR/TA  )  *DUP 
00  TO  114 

116  IF(DP-TA)  117,H7,118 

117  CAL  I  =  .  5  * (PTZ*DP/TA )  »DP 
OO  TO  114 

HP  I  F  (  DP  — DI  IR  )  IIP*  IIP,  140 
IIP  T I MF=DP— T7 

P  T  I  =P  I  *  (  1  ,-T  IMF/DUR  ) *FXP ( -T IME/DUR ) -.  1*00*  (  1  ,-T  I "  F  / D U R  )  **?  . 

**FXP( -?.*T IMF /OUR ) 

CALI=.R*PTZ*TA+.5*(PT7-PTI )*(DP-TA)+DTI*(DP-TA) 

00  TO  114 

140  CALI  =  .8*PTZ*TA4-.5*PT7*(DIJR-TA) 

114  VP  I  TF  (  i  ,  1  8  )  1)1  ?  (  I  1  ,  I  ?  ,  I  1  )  ,CAL  I 

OO  TO  loft 
10  7  COM  TIN!  IF 
106  CONTINHF 
105  CONTINHF 

CALL  L  INK ( CALLR  ) 

Q  FORMAT (  1 X ,  •  CR  I  T  I  C AL  PFRIOD  FOR  SCHOOL  PRICK  WALL’,/) 

11  FORMAT ( 1F1 5. ? ) 

1?  FORMAT  (  1  X  ,  t  M I  N  I  mij  M  DISTANCF  for  F  =  »,F10.?,*T  STRFMOTH=  •  ,  E  1  0.?  ,/  ) 

11  FORMA  T  (  1  X  ,  •  CR  I  T  ICAI  DFPIOn  AN~  IMO||L$F  FOR  SCHOOL  ROOF  »  ♦  /  ) 

14  FOPma T (  1 X  ,  • F •  , 1  OX  ,  • TONS  I Lf  STRFNOTH  •  , 1  OX  ,  • C  PERIOD*,) 

15  FORMAT ( 1 X, • CHARGE  WF I GHT  =  * , F 1 0 . ? , • D I  ST ANCF= *  »F 1 0 . ? ,  *UN  I  T  IMPULSE^*, 

*  F  1  o  »?  »  /  ) 

16  FORMAT (  1 X  ,* E  '  *8X  ,* WF  IGHT *  ,8X  ,* M.OF  RUPT*,8X,*C  PFRIOD', 8X, 

*  ' C  IMPUL  SF  '  ,  ) 

17  FORM  AT ( FP.O ,4F 1 ?.  7 1 

IQ  FORM  A T (  1  OX ,' SCHOOL  ROOF  srcTION',/) 

18  FORM a T (  1 SX , ' RFQUI RFD  I MPUL IF = '  , F I  0 . 7  ,  7X ,  •  C A l C  .  I MP"L SF =  '  ,  F ] 0 . ?  ,  ) 

47  FORMAT (5X* • INC.  OVPRFS.  =  • ,F 1 0. ? ,  * PFAK  RFFL . PRES . =  •  , F 1 0 . ?  , 

*  • DUR.  =  ' , FI  0. ?  ,  /  ) 

50  FORMAT  (  1  X  ,'  FRAGMENT  SIZES  AT  d  I  =  •  ,  F I  0 . 7  ,  •  A  RE  TOO  SMALL  TO  CAUSF  '  , 

*  1 X  ,  ' SFR I OUS  INJURY',/) 

51  FORMA T ( ) X, ' FOR  DI  =  '  , F10. 7  ,?X , • PROR.OF  SFRIOUS  WOUNDS =0 .',/ ) 

5?  FORMAT () X, ' FOR  PI  =  •  ,F10. 7 , ?X , • PROR.  OF  SFRIOUS  WOUNDS =', E 1 0 . 2 »/ ) 

51  FORMAT ( 1  X  ,' PRORAR IL I TY  OF  GLASS  FRAGMFNT  INJURY?',) 

54  FORMAT ( 1 X ,' LOS/SO  IN'  , 1  OX  ,  ' LRS/SO  IN'  » 1  OX , ' MSFC '  , / ) 

58  FORM A T (  1 X  ,  ' L RS/SO  I N '  , 6 X ,  ' l  RS /F T  •  , 6 X , ' l RS /SO  IN  •  »6X  *  •  MSFC  '  , 

*  6X ,  •  PS  I -MS '  , / ) 

FND 

//  DUP 

♦STORF  WS  UA  SCHOL 

//  JOR 
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OHAFW 


//  DUP 
♦DELETF 
//  FOR 

♦ONE  WORD  INTEGERS 
*AR I THMFT IC  TRACF 
♦TRANSFER  TRACF 
*  IOCS ( CARD, IIS?  PRINTFR) 

COMMON  PO(s0),TaU(s0),CPR<S0),HRF<4)  »HRW( 4) »HRS( 4 )  ,HFE<  A ) , 

* HFW ( 4 ) , HFS  t  A ) » HRF ( 4 ) ,HRT ( A ) »CRF ( A ) * CRW ( A ) » CRS ( A ) » CDF ( 4 )  * CDW ( 4 ) » 

♦  CDS! 4) »SrWF( 4 ) »SRT ( 4 ) »SRF ( 4) »SRW (4 ) *SRS( 4) »OHF<  4 ) ,nHT  (  A ) »OFF  (  A) , 
*CTG(A) »TRW(A)  *CPW(A)  *  T  R  C  G ( A )  *  TrWT (4  )  ,DORS ( A )  *  W M  AN ( A )  »PMAN ( A )  » 
♦OFT(4)*IA<ll)»K»NOT,Ul(4»4»4)»TC(4»4»4)»CS(S)  » BCG ( A) »RW( A)  ,CCG<  4)  , 
*AIRF<4)»AIRF<10)»TC?<4»4»4)»TC3(4»4»4)»UI?<4*4»4),UI3<4*4»4) 

c . 

C  PART  VI.  OFFCF  SURROUTINF 

c . 

C . FIXFD  PARAMFTFR  VALIJFS . 

A=l./0. 

FFFT  =  R6 • 

D  I  =  2 1  6  . 

TW  =  fl. 

G=0B6. 

AIP=A2.7 
I  I  I  =1 

C . DFTFRMINFS  NO.  OF  PARAMFTFRS  FOR  OFFICF  HALF  WALL . 

DO  APO  1  =  1  ,4 
IF(OHF(I))  401,401*400 
AO  1  K 1 = I  —  1 

GO  TO  AO? 

AOO  CONTINUF 
K  1  =  4 

c . 

C . CRITICAL  PFRIOD  FOR  OFFICF  HALF  WALL . 

C . 

AO?  DO  ADO  1=1  ,K  1 

E  =OHF (  I >*1000000. 

DO  AOO  I  1  =  1 ,K1 
TF=OHT (II) 

S I G  =  1 OOOOOO. 

W  =  .  ? 

PP  =  o. «7*S0RT ( (  F*AIP*G) / (W*DI**A. )  ) 

P  =  PP/S  .<96 

AOO  TC(  I  *  1 I»  1 1 1 ) =  < 1000./ (A.*P) )*.67 
WRITE( 0*13) 

WRITE! 3, 14) 

WRITE! 0,19) 

DO  AO A  1=1 »K1 
DO  AO A  J=1 , K 1 
I  I  I  =1 

F  =OHF (  I  )*lOOOOOO. 

AOA  WRITE! 0,11  )  E , OHT ( J ) , TC (  I  ,  J , I  I  I  ) 

c . 

C . DETERMINES  MINIMUM  STAND-OFF  STORAGE  DISTANCE  FOR  EACH 

C . CHARGE  SIZE  AND  SET  OE  PARAMETER  VALUES . 

C . 

C . VARIATION  OF  PARAMFTFR  VALUFS . 

DO  AOS  1=1,10 
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u  u 


DO  J  =  1  ,K1 

C  P  =  TO  (  I  ,  J, I  I  I  ) 

E  =OHE (  I )*10000nn. 

WRITFIB.l?)  E  *OHT ( J ) 

C . VARIATION  OF  CHARGF  SIZF . 

DD  AOfi  L=1 .S 

C . VARIATION  OF  SCALFD  DISTANCE . 

DO  AO 7  L  1  =  1  »  S  0 

C . BLAST  PARAMFTFRS  FOR  SCALFD  DISTANCE 

P I=PO( Ll  ) 

Dl)R  =T  AIJ  (  L 1  )*CS(L)**A 
I F ( L  1-07 )  A0  8  ,  A  OB  *  A  OR 
AOB  S  =  L 1 

R  =  ( S  +  B .  ) *rS( L ) **A 
GO  TO  AIO 
A  0  q  S  =L  1 

R= (  (  S  — 07 . ) *S.+AO.  )*CS(L ) **A 
AIO  00  =  ?. S*( PI**?.  ) / ( 10?.9  +  P  I  ) 

VEL  =  1  I B7. *SQRT ( 1 .  +  ( 6 . *P I) / ( 1 0? .9+P I )) 

PR  =  ?.*PI*( 10?.9  +  A.*PI  )/( 10?. P  +  P  I  ) 

.CRITICAL  IMP.  ON  RRICK  WALL  IS  FUNCTION  OF  BLAST  LOADING 

.ON  TOP  OF  WALL . 

PV=?6.-L 1*. I S 
TST  =OHT ( J) 

Ol  =  (A.*TW)*( TST  *TW+P V ) / (B.*FFFT ) 

Y 1=  (S.*FEET**B.*01  ) / (B8A.*E*AIP ) 

0?=  (  A. * ( TW-YI  ) *(?. *PV+q.B )  ) /FFET 
C0?= (A.* ( TW-S. ) *( ?.*PV+q.B ) ) /FFET 

WORY=  (  .  S*Yl  *01  )  +  (  .  S*  (  7.-Y]  )  *  (  0?-CQ->  )  )  +C0?*  (  7  .-Yl  ) 

V SO = WORK / . OsB 
V  =  SOR  T  t  V  SO ) 

T  I  =  (?S.6/G)*V-t-.0?A*V 
CI=TI*1000./AB. 

T  S  = ( 0 • *  A . / VFL )* 1 000. 

P  TS  =P I  * ( l.-TS/DUR)*EXP(-TS/DUR)+  00* ( I  .-Ts/DUR)**?.*FXP(-?. 

**TS/D(JR  ) 

C . CALCULATES  POINT  OF  MINIMUM  FFFECTIVF  PRESSURE . 

D  I  V=  (  DllR-TS  )  /SO  . 

DO  AO  1  17=1* sn 

S  =  I  7 

TFND=TS+S*DIV 

PFND=P I*  (  1 .-TEND/DUR)*FXP(-TEND/DUR)+Q0*(1 .- TFND/DIIR )**?.* EXP (-?. 
**T  FND/DUR ) 

I F ( PEND- • ? )  ASB.ABB.ABl 
A  0 1  CONTINUE 

ABB  IF(TS-CP)  All, All. Al? 

C . C.  PFRIOD  LESS  THAN  TImf  to  STAG . 

Al?  A VFP= (PTS+PR ) /? . 

I F ( AVER- • S )  A?B,A1B,A1B 

A  1 B  CALF  =  . S*TS* ( PR-PT S ) + TS*P T S+ . S*P TS* ( TEND-TS) 

GO  TO  A  1 9 

C . C.  PFRIOD  GREATER  THAN  TIMF  TO  STAG . 

All  IF(CP-TFND)  A16.A16.A1S 

A  1 6  PCP  =  Pl*(l.-rP/DURl*FXP(-CP/DUR)+  OO* (  1  . -CP /OUR ) **? . * F x P ( -? . 

**CP/DUR ) 

AVFR=r(PR  +  PCP  )  r? . 

IF(AVFR-.S)  A?R  , ABO, ABO 
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430  CALF=TS*PTS+.5*  ( PR-P TS ) * TS+ . 5*  (PTS-PCP)*  (CP-TS  )  +PCP*  (  CP-TS  ) 
*+.5*PCP*(TFND-CP> 

GO  TO  419 
415  AVE  R  =  PR /  2  • 

IFIAVFR-.5)  428*418,418 

418  0 ALF=.8*( PR-PTS)*TS+TS*PTS+.5*PTS* <  TFND-TS ) 

c . blast  loading  on  rfap  of  half  wall . 

419  TA=( 16.67/VEL >*1000. 

P  T  A  =P  I  *  ( l.-16./DUR)*FXP(-16./DUR)-.3*QO*(l-16./DUR)**2. 
**FXP(-32./DUR> 

IF(TFND-TA)  480,480,481 

480  CALR=.5*(PTA*TFND/TA)*TFND 
GO  TO  482 

481  T5=TA+.5*< TFND-TA ) 

PT5=PI*  (  l.-T5/DUR)*FXP(-T5/DUR)-.3*00*<  I  ,-T  5  /DUR  )  **2  • 

**FXP( -2.*TS/DUR ) 

CALR  =  . 6*PTA*TA+.5*  < PTA-PT5 )* ( T  5 -T  A ) +P  T  5* ( T 5- T A )  +  . 5 *PT 5 * ( T FND-T 5 ) 

482  CAL  I  =  CALF-CALB 
IF(CI-CALI)  407,428,428 

C . MINIMUM  DISTANCF  FOR  GIVFN  CHARGE  SFLECTFD . 

428  WRITF(8,18)  GS(L)*R*CI 
WRITFI 8,17)  PI, PR, DUR 

C . GLASS  FRAGMENT  INJURY  MODFL . 

WRITE! 8,53) 

IFIPI-5.)  470,471,471 

471  WR I TF ( 3,50)  PI 
GO  TO  472 

470  FW=-.75*PI+3.9 
FV=18.7*PI+62.5 
CV=230.*FXP(-1.9*FW)+103. 

IF(CV-FV)  473,474,474 
474  WRITF! 3,51)  PI 
GO  TO  472 

473  M  =  (  (FV-CV) /50. ) 

S=M*. 1+. 1 

WR  I  T  F ( 3 , 52 )  PI,S 

472  GO  TO  406 
407  CONTINUF 
406  CONTINUE 
405  CONTINUE 

CALL  LINK! CALLR ) 

11  FORMATI3F15.2) 

12  FORMAT  (  IX,  •  MINIMUM  DISTANCF  FOR  E='*F10.2»'T  STRFNC,Th=  •  ,  F 1  0, 2  ♦  /  ) 

13  FORMAT (1 X, ' CRIT ICAL  PFRIOD  FOR  OFFICF  HALr  WALL',/) 

14  FORMAT ( IX, • F' *  I  OX , • TFNSILF  STRFNGTh* *10X,'C  PERIOD'*//) 

15  FORMATMX,  'CHARGE  WEI  GHT  =  '  ,  F9. 0  ,  •  D  I  ST  ANCF=  '  ,  F  1 0 . 2  *  •  UN  I  T  IMPULSE='* 
*F10«2,/ ) 

17  FORMAT (5X»  •  INC.  OVPRES,  =  ' *F10.2*  'PEAK  REFL . PRES .  =  ' , F 1 0 . 2  * 

* ' DUR.  =  • *  FlO.2,/ ) 

19  FORMAT ! IX, ' LRS/SO  I N ' , 1  OX , ' LBS / SO  I N '  , 1  OX , ' MSEC ' , / ) 

50  FORMATIlX, 'FRAGMENT  SIZFS  AT  P I = • »F 10 . 2 » • ARE  TOO  SMALL  To  CAUSE', 

*  1 X  » 'SERI OUS  INJURY',/) 

51  FORMAT ! 1 X, • FOR  P  I  =  ' , F 1 0 . 2 , 2X , ' PROR .  OF  SERIOUS  WOUNDS  =  0 . '  , / ) 

52  FORMAT ! IX,  ' FOR  P  I  =  ' * F 1 0 , 2  * 2X , ' PROR . OF  SFRIOUS  WOUNDS= •  , F 1 0 , 2  * / ) 

53  FORMAT (IX, 'PROBABILITY  OF  GLASS  FRAGMFNT  INJURY',  ) 

END 

//  DUP 
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*STORF  WS  IJA  OHAFW 

//  JOB 
//  nup 

•DELATE  OFIJLL 

//  FOR 

•ONE  WORD  INTEGERS 
•ARITHMETIC  TRACF 
•TRANSFER  TRACF 
•+OCS < CARO, 1  1  A?  PRINTFR) 

COMMON  PO( SO) , T AU(SO) ,CPR<SO) »HRF< A ) »HRW<  A ) tHRS<  A ) *HFF<  A ) , 

*  HFW ( A )  *HFS ( A ) *HPF<  A )  *HRT ( A )  *CRF  < A  > *CRW( a ) * C R S  t  A ) »COF< A )  ♦  C D W ( A )  » 

*  CDS ( A ) ♦ S  RW  F I  A )  *SBT ( A  )  *  SRF I  A ) »SRW( a )  *  SR  S ( A )  *OHE<  A ) *OHT ( a )  *OFF  <  A  )  * 

*  CTG  t  A  )  »  TRW  (  A  )  *CPW  t  A  )  »  TRCC.I  A  )  »  TRWT ( A )  *  DORS  t  A  )  *WMAN  C  A)  *  PM  AN  (  A  )  » 

»OFT(A)»IA<ll)*K*NOT»UT(A»A»A)»TC(A*A*A)*CS(S)  ,rcG ( A ) *BW( A ) *CCG  < A  )  * 
*AIRF<A)*AIRF<10)*TC2<A*A*A)  *  TC A  <  A  *A  *  A  )  »UI? I  A  *A  *A )  »UI  A( A* A  » A) 

c . 

C  PART  7  OFFICF  HILL  WALL  SUBROUTINE 

. . 

C . FIXFD  PARAMFTFP  VALUES . 

A  =  1  ./A. 

FFFT  =  Rf> . 

A  I  =  A?  •  S 
G=AS6. 

I  I  1  =  1 

DO  Son  1=1, A 

C . DFTFRMINFS  NO.  OF  OFFICE  FULL  WALL  PARAMETERS . 

I f ( OFF ( T  )  )  so i ,sni  , s no 
sn i  <i=i—l 

GO  TO  SO? 
son  CONTINUE 
K  1  =  A 

C . 

C . CRITICAL  PFRIOO  FOP  OFFICF  FULL  WALL . 

C . 

so?  no  sos  i =i ,< i 

f=opf ( i ) * i nnnnor . 

DO  sns  I  1  =  1 ,<1 
T  S  =  OF  T (  I  I  ) 

W=.  ?1 

PP=q. fl 7* SORT t ( F*A I *G > / ( W»FEFT**a  . ) ) 

P=PP/S.06 

so-5  TC(  I  ,  1 1  ,  1 1  I)=inno./u.*P) 

WRITF( s, 17) 

WP I TF ( S , 1  A ) 

WR  I TE ( A  » 1 R ) 

DO  Soa  1=1  , < 1 
no  SOA  J=1 , K 1 
1 1 1  =1 

f=off( i )*ioonnnn. 

SOA  WP I TE ( s , 1 1 )  F *OFT t J ) , TC t I , J  ,ITI) 

c . 

C . DFTFRMINFS  MINIMUM  STAND-OFF  STORAGF  DISTANCF  FOR  FACH 

C . CHARGF  SIZE  AND  SFT  OF  PARAMFTFR  VALUFS . 

C . 

C . VARIATION  OF  PARAMETER  VALUES . 

DO  SOS  1  =  1, <1 
DO  SOS  J=1,<1 
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C  P=  T  0 (  I«J*III) 
e=ofF( I )*ionoonn. 

WRITE! 3*1?)  E,OPT(J) 

C . VARIATION  OF  CHARGF  SIZE . 

DO  S06  L  =  1 , 5 

C . VARIATION  OF  SCALED  DISTANCE . 

DO  SO  7  L  1  =  1 , 5  0 

C . BLAST  PARAMETERS  FOR  SCALED  DISTANCE 

PI=PO( LI ) 

DUR  =  TAIJ(L1)*CS(L)**A 
IF ( Ll-37 )  SOB, 508, SOP 
508  S*L1 

R= ( S+3. ) *CS ( L ) **A 
GO  TO  510 
SOP  S=L 1 

R  = (  (S-87. )*5 .+40. )*CS(L )**A 

510  QO=?.5*(PI**2.)/(102.9+PI) 

VFL=1187.*S0RT ( 1 .+(6. *PI ) / < 102.9+PI ) ) 

PR  =  2  »  *P  I  *  (  10?.9+4.*PI  )/< 102.9+PI  ) 

W0RK=708. 

VSO=WORK/.0163 
V=SORT ( VSO) 

T 1= (9.4/6 >*V+.008*V 
CI=TI*l0nn./48. 

TS= ( 8. *8  ,/VFL )*1 0  00. 

IF(DUR-TS)  511, 51?, 51? 

C . CASE  I.  DURATION  OE  BLAST  LESS  THAN  T I MF  TO  STAG. 

511  CALI=.5*PR*DUR 
GO  TO  S 1 4 

C . CASE  II. TIME  TO  STAG.  LFSS  THAN  DURATION 

SI?  PTS  =P I  * ( 1 ,-TS/DUR ) *FXP ( -TS/DUR) +00* ( 1 ,-TS/DUR )**•>.  *rXP  ( -2  . 
**TS/DUR ) 

D I V  = ( DUR-T  S ) /S  0 . 

DO  S 1 3  17=1, SO 
S  =  I  7 

TEND=TS+S*DIV 

PEND=PI*( l.-TFND/DUR) *EXP( -TEND/DUR) +00*( 1 .-TFND/DUR ) **2 .*EXP ( -2 . 
*  *TEND/DUR ) 

I E ( PEND- • 2 )  S19,S19,S13 

513  CONTINUE 

519  IE(TS-CP)  515,515,516 

C . CASE  I.  CRITICAL  PERIOD  LESS  THAN  TIME  TO  STAG . 

516  AVER=(PTS+PR)/?. 

I E ( AVER-5. )  517,518,518 

518  CALI  =  .S*(PR-PTS  ) * TS  +  TS*P TS+ . 5*PTS* ( TEND-TS ) 

GO  TO  514 

C . CASE  II.  C.  PFRIOD  GREAT  FR  THAN  T IMF  TO  STAG . 

515  I F ( CP-TFND )  S?0,S20,S?I 

520  PCP=PI*(1.-CP/DUR)*EXP(-CP/DUR)+00*(1.-CP/DUR)**?.*FXP(-?. 

**CP  /Dt)R  ) 

AVER=(PR+PCP)/?. 

I E ( AVER-5 • )  517,  522,522 

522  CALI=.5*(PR-PTS )*TS+TS*P TS+PCP* ( CP-TS ) + . 5* ( PTS-PCP ) * ( CP-TS ) 
*+.5*PCP* ( TEND-CP ) 

GO  TO  S 1 4 

521  AVFR=PR/2. 

I F ( AVER-S . )  S17,5?3,S23 
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523  CALI  =  .5*(PR-PTM*TS+TS*PTS+.5*PTS*(TFND-TS) 

c . minimum  distancf  for  givfn  charge  selected . 

51A  IF(CI-CALI)  507,517,517 

C . MINIMUM  DISTANCE  FOR  GIVFN  CHARGE  SELECTED . 

517  WRITE(R*15)  CS(L>*R*CI 
WR  I  TF  (  2 ,  37  )  P  I  »PR  *DUR 
GO  TO  5 D 6 
507  CONTINUF 
506  CONTINUF 
505  CONTINUE 

CALL  LINK(CALLR) 

11  FORMAT (F10.0*2F15.2) 

1?  FORMAT ( IX* 'MINIMUM  DISTANCF  FOR  E= • Fp . 0  * 5X , ' T  S T REMG TH= ' , F 1 0 . ? » / ) 

14  FORMAT ( I  X » • F ' ♦ 1  OX » ' TENS  I LF  STRENGTH 10X C  PFRIOD'*//) 

15  FORMAT! 1 X. 'CHARGF  WF I GHT= • ♦ FO . 0  ,  ' D I  ST ANCF= » , F 1 0 . 2  * 

♦  'UNIT  IMPULSE*'  ♦FlO.2*/) 

17  FORMAT! 1X*'CRIT  ICAL  PERIODS  FOR  OFFICE  FULL  WALL'*/) 

IP  FORMAT ( IX* 'LBS/ SO  IN' *  1  OX  * ’ LRS/SO  I N ’  *  1  OX  * ' MSEC '  *  1 

37  F ORM A T ( 5  X i '  INC.  OVPR FS .  =  ' * F 1 0 . 2 * • PE AK  REFL . PRES .  =  ' . F 1 0 . 2  * 

* • DUR • = ’ *FlO. 2, // 1 
FND 

//  DUO 

* STORE  WS  UA  OFULL 

//  JOR 
//  DUP 

♦DELETE  BUS 

//  FOR 

♦ONE  WORD  INTEGERS 
♦ARITHMETIC  TRACE 
♦TRANSFFP  TRACE 
♦lO^StCARO*! ,3?  PRINTER) 

COMMON  PO (50) , T  AU( 50 1 *CPR  (  50) *HRF ( A ) ,HRW( 4 ) , HRS ( 4  )  * HFF ( 4 1  * 

*HFW(4)  ,HFS<4)  ♦  HRF  (  A  )  *HBT  (4) *CRF<4)  *CRW(4)  *CRS(4)  ,CDE<A)  ,CDW(4)  * 

♦  CDS(4)*SRWF(4)  *SrT ( 4  1  *SRF ( A )  ,SRW( 4 )  »SrS(4)  »OHF (4)*OHT(4)»0FE<A)» 
*CTG( 4) *  TRW! A )  *CPW ! a  )  »TRCG( A )  »TRWT ! 4 )  *DORS ( 4  )  *WMAN ( 4 )  ,  PMAN ( 4 )  * 

♦  OFT (4 ) ,  I  A (  11  )  ,  K  ,N0  T  tUl  (4*4*4)  ,TC(A*A,4)  * CS  (  5 )  »BCG l A )  *  BW ( A )  *CCG  <  A  1  * 
♦AIRF<A)»AIRF<lO)»TC2<A»A*A)*TC3(A»A»A)*UI2<A*4*A)*UI3<4»A*A) 

c . 

C  PUS  OVERTURN  SUoPOUTTNF 

c . 

C . FIXED  PARAMETER  VALUFS . 

A  =  1  .  /  3  . 

RASF=42. 

I  I  I  -1 

C . DFTFRMINFS  NO.  OF  BUS  PARAMFTFRS . 

DO  600  1=1*4 
I F ( RCG! I ) 1601,601 ,600 

601  K 1 = I— 1 

GO  TO  602 
600  CONTINUE 
K  1  =  A 

C . 

C . CRITICAL  IMPULSE  FOR  RUS 

C . 

602  DO  603  1=1  , XI 
CG=  RCG (  I  ) 

DO  603  I  1  =  1, K1 
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W  =  BW ( I  I ) 

DIST=SQRT(BASE**2.+CG**2. ) 

R =D I  ST* .939 
S  =  D  I  ST *.  342 
R  I S  F  =  R  — C  G 
SLID=BASF-S 

WORK  =  BW( I  I ) *R ISF-SL ID*.364*BW(  I  I ) 

AI  =  (RW(  I  1 1/386.  )*t(  R6.5*#2.+90.  **?,  )  / 1 2  -  +60.**2  .  ) 

A  I A  =• 5*A I 
0MGA=W0RK/AIA 
OMGA=SORT (OMGA ) 

PUL  SE= ( A  I *OMGA ) /72. 

60S  UT(T*TT*TTI ' =PUL-SF*1000. /SO4OO . 

WR I TF( 3  » 18 ) 

WR I TE (3,19) 

W R  I T E (  3,25) 

DO  604  1=1, K1 
DO  604  J=1 ,K1 
W=BW ( I ) 

CG=RCG( J ) 

CI=UI ( I , J,  1 1  I  ) 

604  WRITE( 3,20)  W,rG*GI 

C . 

C . DETERMINES  MINIMUM  STAND-OFF  STORAGE  DISTANCE  FOR  EACH 

C . CHARGE  SIZE  AND  SET  OF  PARAMETER  VALUES . 

C . 

C . VARIATION  OF  P ARAMETFR  VALUES . 

DO  60S  I =1  ,K1 
W  =BW (  I  ) 

DO  60S  J=I  ,K1 
CG=RCG ( J ) 

WRITE! 3*24)  W  »  CG 

C . VARIATION  OF  CHARGE  SIZE . 

DO  606  L  =  1 » 5 

C . VARIATION  OF  SCALED  DISTANCE . 

DO  607  L  1  =  1 ,50 
C I  =  UI ( I , J, I  I  I ) 

c . blast  paramftfrs  for  scaled  distance 

P  I  =  PO( LI ) 

PR  =  2  »  *P I  * ( 102,9+4, *P I )/( 102,9+PI ) 

DUR=TAU(L1 )*CS(L>**A 
IFtLl-37)  608,608,609 

608  S  =L 1 

R  =  ( S+3 . ) *CS  t  L ) **A 
GO  TO  610 

609  S=L1 

R  =  t  (S-3  7. )*S.+40. ) *CS ( L ) **A 

610  Q0=2.5*(PI**2. )/( 102.9+PI ) 

VEL=11R7.*SQRT( 1 • + ( 6. *P I ) / t 102,9+PI ) ) 

TS= (3,*10, 5/VEL 1*1000. 

I  Ft DUR-TS 1611,61 2,612 

C . CASE  I.  DURATION  OF  BLAST  LESS  THAN  T I MF  TO  STAG. 

611  CAL  I=.5*PR*DUR 
GO  TO  622 

C . CASE  II.TIMF  TO  STAG.  LESS  THAN  DURATION 

612  PTS=PI*(  l.-TS/DUR)*EXP(-TS/DUR)+QO*< 1 ,-TS/DUR ) **2 ,*FXP (-2 . 
**TS/DUR ) 
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C . CALCULATES  POINT  OF  MINIMUM  EFFFCTIVF  PRESSIJRF . 

0  I  V  =  (  01 IR-TS  )  /SO. 

DO  613  I  7=1  ,  «S0 
S  =  I  7 

T  FND=T  S  +  S*D I V 

PEND=P I* ( 1 . -TFND/ DUR ) *FXP (-TEND /DUR )+Q0* ( 1 TEND/DUR ) ** 2 . *EXP ( -2 
**TFNP/DUR> 

I F ( PEN  D-  •  2 )  614,614*613 

613  CONTINUE 

614  CALF=.S*(PR-PTS) * T S+TS*P TS  + . 6* (PTS-.l )*< TFND-TS )+ ( TFND-TS) *PFND 

C . BLAST  LOADING  ON  RFAR  OF  RtJS 

TA=((7.2+42.)/VFL)*l000. 

T I MC  =  T  A- (7.7 /VFL ) *1000. 

PT  A  =P I  * (  l.-T IMF/DUR ) *EXP < -T I  ME /DUR ) 3*00* (  1  ,-T IMF/DUR )  **  2. 

**FXP ( -2 .*T I MF/DUR ) 

IF(TEND-TA)  616,619,619 
616  PT I =PTA*TFND/T IMF-7. 2*PT  A  /  4  2  • 

CALR=. S*PTI*(TFND- (7. 2*1000. /VFL)) 

GO  TO  621 

61°  CALR=.p*PTA*TIMF+.6*PTA*(TFND-TIMF) 

621  CAL I=CALF-CALB 

C . MINIMUM  DISTANCF  FOR  GIVEN  CHARGE  SFLECTED . 

622  IF(CI-CALI)  607,620,620 

C . MINIMUM  DISTANCF  FOR  GIVEN  CHARGE  SELECTED . 

620  WRITE! 3,21 1  C S ( L )  , R 

WPITF(3,17)  PI, PR, DUR 

C . GLASS  FRAGMENT  INJURY  modFL . 

WRI TF< 3, S3 ) 

I F(  PI-6.  )  670,671  ,671 

671  WPITF(3,S0)  PI 
GO  TO  677 

670  FW=— . 7S*P I +3 .9 
FV= 1 8. 7*P I +62. S 
CV=?30.*FXP(-l.o*FW 1+103. 

IF(CV-rV)  673,674,674 
674  WRI TF ( 3 , 61 )  PI 
GO  TO  67? 

673  M=< (FV-CV ) /SO. ) 

S=M*. 1  +  .  1 
WRITE!  3,S2)  S  ,P  I 
67?  GO  TO  606 
607  CON  T I NUF 
606  CONTINUF 
60S  CONTINU‘D 

CALL  LINK(CALLR) 

17  FORMATfSX,'  INC.  OVPRFS .  =  • , F 10 . 2 , *  PE AK  RE F L . P RFS .  =  '  ,  F 1 0 . 2  . 

* ' DU  R.=  ',F10.?,/) 

18  FORMAT! 1 X  ,' CRI T ICAL  IMPULSE  FOR  BUS  OVERTURNING',//) 

19  FORMAT ( 1 X  ,' BUS  WF I GHT •  , 1  OX  ,  ' C  GR A V  I T Y '  ,  1  OX  ,  •  C  IMPULSE',) 

?6  F  ORMA  T ( 1 X  ,  ' LRS '  ,10X,'  INCHES' ,10X,  'P  SI -MSFC ' , // ) 

20  FORMAT ( 3F1 0. ? ) 

21  FORMAT ( 1 X, ' CHARGE  WF I GHT = ' , F 1 0 . ?  ,  1  OX , ' D I  ST ANC F= '  , c  1  0 . 2  ,  /  ) 

24  FORMAT! IX, 'MINIMUM  DISTANCF  FOR  WF I GHT  ='  ,r10.?.10X»'CG='  ,  F  10  »  2  , 
SO  FORMAT ( l X , ' I NSUFF IC I FNT  DATA  CONCERNING  PLATE  GLASS  FRAGMENTS', 
*1X,  'AT  PI  =  ' ♦F10. ?,/) 

31  format ( 1 x ,< usi ng  window  g.  data, prop,  of  serious  wounds  is  o.'» 

*1X,  'AT  PI=* ,Fin.?,/) 
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52  FORMATt IX* 'USING  W.  G.  DAT  A  *  PROR  .  OF  SFR.  WOUNDS= ' , F 10 . 2 » 

♦  IX,  'AT  PI  =  '  ♦FlO.?,/  ) 

53  F0RMAT(1X*'PR0+A+ILITY  OF  GLASS  FRAGMENT  INJURY',  ) 

END 

//  DUP 

♦STORE  WS  UA  RUS 

//  JOR 
//  DUP 

♦DFLFTE  CAMPR 

//  FOR 

♦ONE  WORD  INTEGERS 
♦ARITHMETIC  TRACE 
♦TRANSFFR  TRACF 
*  IOCS ( CARD* 1 132  PRINTFR) 

COMMON  PO (50)  »  T  AU ( 30) ,CPR(50)  ,HRF ( 4)  ,HRW( 4 ) ,HRS( 4 ) *  HFF ( 4 )  » 

♦  HFW (4)  ,HFS<  4) * HRE ( 4 ) *HRT ( 4)  ,CRF ( 4 ) *CRW(4 ) ,CRS(4 )  ,CDE<4 ) *CDW(4  )  * 

♦  CDS ( 4 ) »  SRWE ( 4 ) *  SRT ( 4 ) , SRE ( 4 ) *SrW( 4 )  ,SrS ( 4)  »OHF (4)  *OHT ( 4 )  *nFE  <  4 ) » 

♦  CTG( 4) *  TRWt  4)  *CPW(4  )  *TRCG( 4)  , TRWT (4 )  , DORS (4 )  , WMAN (4)  ,PMAN (4)  * 

♦  OFT (4) *  I  A ( 11 ) *K  ,NOT ,UI ( 4*4*4)  ,TC(4,4,4)  »CS(5  )  »BCGt  4)  ,BW(4 ) »CCG( 4) , 

♦  AIRE<4)*AIRF(  10),TC2<4*4*4)  *TC3<4*4*4)  »UI ? (4  *4*4)  ,UI3<4*4*4) 


CAMPER  OVFRTURN  SURROUTINF 


C.  , 

c 

c . 

C . FIXFD  PARAMFTER  VALUFS . 

A  =  1  .  /  3  • 

BASF=33. 

I  I  I  =] 

C . DETERMINES  NO.  OF  PARAMETERS  FOR  CAMPER. 

DO  700  1=1,4 
I F (  CCG ( I ) )  701,701 ,700 
701  K 1=1-1 

GO  TO  702 
700  CONTINUF 
K1  =4 


C . CRITICAL  IMPULSE  EOR  CAMPER  AND  TRUCK  UNIT . 

C . 

702  DO  703  I = 1  * k 1 
CGT=CTG (  I  ) 

C  GC  =  CCG (  I ) 

DO  703  I  1  =  1, K1 
WTR  =  TRW(  I  I ) 

WCP  =  CPW(  I  I ) 

D I  ST 1 =SQRT ( CGT^»? ,+BASE**? . ) 
DIST2=SQRT(CGC»»2.+RASF»*2. ) 

R  ISE1=DIST1-CGT 
R  ISF2=DIST2— CGC 
WORK  =WTR*R  I  SF  1  +WCP*R  I  SE  2 

A  I  =  t  (WTR+WCP ) /3R6. )* (  t 7? . *♦? . +54 . **? .  )/  12. +  2612.) 

AIA=.5*AI 

OMGA=WORK/AI A 

OMGA=SQRT ( OMGA ) 

PULSE=  A I *OMGA/60 • 

70  3  UI  t  I  ,  I  I  , I  I  I  )=PULSE» 1000./ 17640. 

WRI TE( 3 ,20 ) 

WRITFt  3*22) 

WRI TE( 3,25 ) 
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DO  704  I =1 ,K1 
DO  704  J=1 ,K1 

7O4  WRITE<8»11>  CPW(  I)  ,TRW<  I  )  ,CCG<  J)  *CT6<  J)  .UK  I  »J,  1 1  I  ) 

c . 

C . DETERMINES  MINIMUM  STAND-OFF  STORAGF  DISTANCE  FOR  FACH 

C . CHARGF  SIZE  AND  SFT  OF  PARAMETER  VALUES . 

C . 

C . VARIATION  OF  PARAMETFR  VALUES . 

DO  70S  I =1 »K1 
CGT=CTG< I ) 

CGC  =  CCG (  I ) 

DO  705  J=1 ,K1 
W TR  =  TR W  (  J) 
wrP=CPW( J) 

WRITE! 8  ,  ?8  )  WCP ,WTR ,CGC ,CGT 

C . VARIATION  OF  CHARGF  SIZF . 

DO  706  L  =  1  ,  5 

C . VARIATION  OF  SCALFD  DISTANCF . 

DO  707  L 1 = 1 , 50 
C I =UI ( I ,  J, I  I  I ) 

C . BLAST  PARAMFTERS  FOR  SCALED  DISTANCF 

P I  =  PO( LI  ) 

PR  =  ?.»PI*  < 10?.9  +  4.*P I )/ < 1 0?.9+PI) 

DUR=TAU( LI )*CS(L)**A 
I F ( L 1 -37 )  708,708,709 

708  S  =  L  1 

R  =  ( S+3  . ) *CS( L ) **A 
GO  TO  710 

709  S  =  L  1 

R  = (  ( S-’7 . ) *S.+40.  )  *  C  S  (  L ) *  *  A 

710  00=?.S*( PI**?. )/ ( 10?.q+Pi ) 

VFL  =  1 1 87. *SORT ( I . - ( 6 . *P I  ) / ( 1 0?  . 9  +  P I  )  ) 

TS= ( 8. *7 .5 /VFL ) *1000. 

I F ( DUR-TS )71  1,71?, 71? 

C . CASE  I.  DURATION  OF  BLAST  LESS  THAN  TIME  TO  STAG. 

711  CALI=.5*PR*DUR 
GO  TO  7?? 

C . C  A  S  F  II.TI^F  TO  STAG.  LFSS  THAN  DURATION 

71?  PTS  =  PI*(1.-TS/DUR)*FXp  (  -  TS/nUR  )  +  00*  <  1  .-TS/DUR  )  **  •> .  *FXP  (  -5  . 

**TS/DUR ) 

D I V  = ( DUR-T  S ) /SO. 

DO  7] 8  17=1,50 
S  =  I  7 

TFND=TS+S*DIV 

P  FNO  =  P I  * ( 1  . -T FN n/ DUR ) *FX P ( -T  END /DUR ) +00*  <1  .- T END/ DUR ) * *? . *F XP ( - ? . * 
*T  FND/OIIR  ) 

I F ( PFN  0- . ? )  714,714,717 

718  CONTINUF 

714  C ALF= . 5* (pR_PTs )*TS  +  TS*PTS+. 5* ( PTS- . 1  )* ( TFND-TS ) +PFND* ( T  ND-TS  ) 

C . BLAST  LOADING  ON  RFAR  SIDE  OF  CAMPFR  AND  TRUCX  UNIT . 

Ta=  t  I9.+4?. ) /VFL ) *1000. 

T  IME  =  TA-( 9. /VFL )*1000. 

PTA=PI  *(  l  . -T  IMF /DUR  )  *FXP  (  -T  IMF /DUR  )  -.  8*00*  (  -\  .-T  I  ME /DUR  1**7. 

**FXP(  -?.*T  IMF/DIJR  ) 

I F  t  TEND-T  A )  716,719,719 

716  PTI =PTA*TFND/TlME-9.*PTA/4?. 

CAL R=.S*PTI*(TFND-t9. *1000. /VFL)) 
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GO  TO  7?  1 

719  CALR=.A*PTA*TIMF+.5*PTA*(TFND-TIMF) 

771  CAL  I=CALF-CALR 

77?  IF(CI-CALI)  707,770,770 
C . MINIMUM  DISTANCE  FOR  GIVFN  CHARGF  SFLFCTFD . 

720  WRITFt?*?! )  CS( L) *R 
WRITF(3»17)  PI, PR, DIJR 

C . GLASS  FRAGMFNT  INJURY  MODFL . 

WRITEt  3, S3) 

I  F  (  P  I  — 5  •  )  770,771,77] 

771  WRITEt  3, SO)  PI 
GO  TO  772 

770  FW=-.75*PI+3.9 

FV=18.7*PI+67.5 
CV=?30.*FXP(-l.9*FW)+103. 

IF(CV-FV)  773*774*774 
774  WRITEt  3, SI)  PI 
GO  TO  77? 

77?  M=( ( FV-CV ) /SO. ) 

S=M*. 1+. I 
WRITEt1*, S2)  S,PI 
77?  GO  TO  706 
707  CONTINUE 
706  CONTINUE 
705  CONTINIJF 

CALL  LINK(CALLR) 

11  FORMAT ( S  F 1 0 . 2 ) 

17  FORMATtSX, •  INC .  OVPRF S.  =  • , F 1 0. 2 , • PF A K  RFF L . PRFS ,  =  •  , F 1 0 . 7  , 

* • DUR.=  •  , FI  0. ?  ,  /  ) 

20  FORMAT ( 1 X  ,  «CRIT ICAL  IMPULSE  FOR  CAMPFR  O VFR TURN  I NG • ♦ / / ) 

21  FORMATt IX, 'CHARGE  WE  I GHT  =  • , F9 . 0  , 1  OX , • DI ST AN CF= *  , F 1 0 . 2  * / ) 

22  FORMAT  (  IX ,  *C  Wf,T  •  ,  RX*»T  WGT',RX,  •  f  CG‘*RX*'T  CG*»8X*'C  IMPULSE', 
*  ) 

23  FORMATtlX, 'MIN  DIST  FOR  C  WGT = • , F8. ? ,  •  T  WGT = • , FR . 2 » ' C  CG='»F8.2* 

*  I  X , ' T  CG= ' , F8 . 7  ,  / ) 

7  S  FORMAT ( I  X  ,  • LBS'  ,RX ,  • LRS '  ,RX  ,  •  INCRFS  •  ,PX  ,  '  I NCHF  S ' , 8X ,  ' PS  I -MS'  ,/ ) 

50  FORMAT t I  X  ,  •  I NSUFF IC I FNT  DATA  CONCFRN I NG  PLATF  GLASS  FRAGMENTS', 
*1X,  'AT  PI=' ,F1o,?,/) 

51  FORMATt IX, 'USING  WINDOW  G.  DAT  A  ,  PROR,  OF  SFRIOUS  WOUNDS  IS  0.', 

♦  IX  » • AT  PI  =  ',F10.?,/) 

52  FORMATt IX, 'USING  W.  G.  DATA, PROR.  OF  SFR.  WOUN DS= • , F 1 0 . 2  , 

*1X, 'AT  PI=' ,Fin.2»/> 

53  FORMAT ( 1 X ,' PRORARIL ITY  OF  GLASS  FRAGMFNT  INJURY',  ) 

FND 

//  DUP 

♦STORF  WS  UA  CAMPR 

//  JOR 
//  DUP 

♦DELETF  TRLER 

//  FOR 

*ONF  WORD  IMTFGFRS 
* AP I T HMFT I C  TRACF 
♦TRANSFFR  TRACF 
♦IOCStCAPD, 113?  PRINTFR) 

COMMON  P0( SO) ,TAU( SO) *CPR( SO) ,HRE (4) *HRW(4) *HRS( 4) * HFE ( 4 ) ♦ 

*  HFW ( 4 ) ♦ HFS  t  a ) ,HRE  < 4 ) 'HBT  <4 ) 'CRE  < 4 ) »CRW<4 )»CRS(4),CDE<4) 'CDWt  4 ) * 

♦  CDS  t  4) *SRWE<4)  ,SRT  t  4 ) 'SRE ( 4)  , SRW ( 4  )  ,  SRS ( 4 )  »OHE ( 4 )  >OHT (4 )  >OFE ( 4) * 
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»CTG<  4 )  *  TRW (4 )  »GPW ( 4  )  »  TRCG ( 4 )  ♦ TRWT (4 )  *DORS ( 4 )  *WMAN ( 4 )  »  PM  AN  <41  * 
*0FT(4)*IA<11)*K*N0T*Ul<4»4»4)*TC<4*4»4)*CS(S) »BCG < 4 )  * BW ( 4 ) *CCG <  4 ) 
*  A  IRF  *  4  )  »  A  I RF  <  l0i»Tr?(4»4»4)  »  Tc8  <  4  *4  *4 )  *U  I  2  <  4  *  4  *4 )  »U 13(4*4*4) 


HOUSE  TRAILFR  SUBROUTINE 


....FIXFD  PARAMFTFR  VALUES . 

A  =  1  •  /  3  • 

R  A  S  F  =  4  4  • 

I  I  I  =1 

. . . . DFTFRMINFS  NO.  OF  MORILF  HOM  F  PARAMFTFRS 
DO  800  1=1*4 
IF ( TRCG ( I ) )  80] ,801 ,800 

801  K 1 =  I— 1 

GO  TO  802 
800  CONTINUE 
K  1  =4 


CRITICAL  I MPULSF  FOR  MORILF  HOMF 


802  DO  808  I =1 ,K1 
CG  =  TRCG (  I ) 

DO  808  I  1  =  1  ,K1 
W=TRWT (II) 

D I S  T  =  SORT ( RASE**?.+CG**2. ) 

R ISF=DI ST-CG 
WORK  =  R I SF*W 

AI=(W/8fl6.)*(  <1?0. * *2.  +  10R. **?.)/  12. +7 l?s  .  ) 
A  I A  =  . S*A  I 
OMGA=WORK/ A  I A 
OMGA  =  SOR  T ( OMGA ) 

PULSE= A I*OMGA/84. 

808  U I (  I ,  I  I ,  II  I  ) =PULSF* 1 000 • /fi l 200 • 

WRITE! 8,24) 

WPITFI 8,2S  ) 

WRITFt  8,28) 

DO  804  1=1  * K 1 
DO  804  J  =  1  » K  1 
CG=  TRCG (  I ) 

W=TRWT< J ) 

C  I  =111  (  I  *  J  »  I  I  I  ) 

8  04  WRITFt 3, 11)  CG  »W  *C I 


DFTFRMINFS  MINIMUM  STAND-OFF  STORAGE  DlSTANCF  FOR  FACH 
CHARGE  SIZE  AND  SET  OF  PARAMFTFR  VALUES . 


VARIATION  OF  PARAMFTFR  VALUES 
DO  80S  1=1  * K 1 
C G  =  TRCG (  I  ) 

DO  80S  J=1 *K1 
W=TRWT ( J ) 

WRITFI 8 , 29 ) W  »CG 

C . VARIATION  OF  CHARGE  SIZF . 

DO  806  L  =  1  *  S 

C . VARIATION  OF  SCALED  DISTANCE. 

DO  807  L  1  =  1  .SO 
C I =UI (  I »  J* I  I  I  ) 
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C . RLART  PARAMETERS  FOR  SCALFD  DISTANCF 

P  I  =  PO( LI  ) 

PR=?.*PI*( 10?.9+4.*PI  )/ ( 102.9+PI  ) 

DUR  =  T  AU ( LI )*CS( L  )  **  A 
IF(Ll-37)  808*808,809 
808  S  =L  1 

R  = ( R+3  .  ) *C  R ( L ) *  *  A 
GO  TO  810 
8  09  S  =  L  1 

R  =  (  (  R  -  3  7  • )*6.+40. ) *  C  R ( L ) ** A 

810  00=2.6*! PI**2. > / ( 10?.9+P I > 

VEL=] 187. *SORT ( 1 . - ( 6 . *P I ) / ( 1 0? . Q+P I ) ) 

TR=<  3.*10./VEL ) *1000. 

IF ( DUR-T  R )  811,812,812 

811  CALI=.6*PR*DUR 

C . G ARE  I.  DURATION  OF  RLART  LFRR  THAN  T I MF  TO  RT  AG. 

GO  TO  822 

C . CARF  I  I • T I MF  TO  STAG.  LFRR  THAN  DURATION 

81?  P  TS  =  P I  * ( 1 ,-TS/DUR )*EXP (-TS/DUR) +00* ( 1 . - TS / DUR ) ** ? . *F XP ( -? . 

* *T R /DUR  ) 

C . CALCULATER  POINT  OF  MINIMUM  effective  pressurf . 

0  I V  =  ( DUR-T S ) /50 . 

DO  813  17=1*80 
R  =  I  7 

TFND=TR  +  R*DI  V 

PFND=P I  * ( 1  .-TFND/DUR ) *FXP ( -TEND /DUR ) +00*1 1  . -TEND/ DUR ) **? . *FXP ( -? . 
**TEND/DUR ) 

I F ( PEND-  .2 )  814,814,813 

813  CONTINUE 

814  C ALF  =  . 6*( PR-PTR ) *TR  +  TR*PTR+ . 5* ( P TR-PFN D ) * ( TEND-TR )  +  ( TFND-TS) *PFN D 

C . RLART  LOADING  ON  R^AR  OF  MORILF  HOMF 

TA=( ( 1 0, +40. )/VOL)*inno. 

TIMF=TA-(lO./VFL)*100n. 

PT A  =P I  * ( 1 .-TIMF/DUR)*FXP(-TIMp/DUR)-. 3*00* ( 1  .-T I MF/DUR ) **? . 

**FXP ( -?. *TI MF/DUR ) 

I F ( TEND-T  A )  816,819,819 

P 1 6  PTI  =PTA*TFND/T IMF-10. *PT A/40. 

C  AL R= • R*P T I  * ( TEND- ( 1 0 . *  1 000 . / VE L )  > 

GO  TO  821 

819  CALB=.f*PTA*TIMF+.R*PTA*(TFND-TIMF) 

8 ? 1  CAL  I  =  CALF-CALR 

8?3  IF(CI-CALI)  807,8  ?0 , 8  ?0 

C . MINIMUM  DISTANCE  FOR  GIVFN  CHARGE  RFLFCTFD . 

8?0  WRITE! 3,?1 )  CS(L) ,R 

WRITF(3,17)  PI, PR, DUR 

C . GLASS  FRAGMENT  INJURY  MODEL . 

WRITE! 3, 63 ) 

IF! PI-6. )  870,871 ,871 

871  WR I T  Ft  3,80)  PI 
GO  TO  87? 

870  FW=-.78*PI+3.9 
FV=18.7*PI+62.6 
CV=?30.*EXP(-1.9*FW)+103. 

IF(CV-FV)  873,874,874 
874  WRI TE( 3*  61 )  PI 
GO  TO  87? 

873  M  = (  (FV-CV ) / 6 0 . ) 
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S=M*.1+. 1 

WP I T  F ( 3*5?)  S.PI 
87  2  GO  TO  806 
807  CONTINUE 
806  CONTINUE 
806  CONTINUE 

CALL  L  INK (CALLR ) 

11  F  OR  M  A  T (  8  F  1  0 , 2  ) 

17  FORMAT (  5  X , '  INC .  OVPRFS.= ' ,F 10. 2 , 'PFAK  RFFL . PRES . =  • . F 10 . 2 ♦ 

♦»DUR.=' »FlO.?*  /  ) 

21  FORMAT ( IX. ' CHARGF  WEIGHT  = • *  F9*0  , 1  OX , • D I  ST  A NCE=  '  > F  1  0 . 2  .  /  ) 

24  FORMAT ( 1 X. • CR  IT  ICAL  IMPULSE  FOR  HOUSF  TRAILFR'.//) 

26  FORMAT < 1 X, ' TRA ILER  CG • ♦ 1  OX « t TRA ILER  WGT'.lOX.'C  IMPULSF'.  ) 

26  F  ORM  A 1 ( ]Xi • INCHFS' i 1  OX i • LRS • » 10X. • P  S I -MS ' * / ) 

20  FORMAT ( 1 X. 'MINIMUM  DISTANCE  FOR  WGT = ' ♦ F 1 0 . 2 » ' CG= ' . F 1 0 . 2  * / ) 

80  FORMAT ( 1 X. ' INSUFF ICIFNT  data  CONCERNING  PLATE  GLASS  FRAGMENTS'. 

*1 X,  '  AT  P  1= ' , F10.?  ,  /  ) 

81  FORMAT! IX. 'USING  WINDOW  G.  DAT  A » PROB •  OF  SERIOUS  WOUNDS  IS  0.'. 
*1X»  'AT  PI  =  '  .F10.2./  ) 

52  FORMAT! IX,  'USING  W.  G.  DAT  A , PROR  .  OF  SER.  WOUNDS* '  , F 1 0 . 2  , 

♦  IX, 'AT  PI='» Fin.?,/) 

68  FORMAT ( 1 X, ' PRORARIL I TY  OF  GLASS  FRAGMENT  INJURY',  ) 

FND 

//  DUP 

♦STORE  VS  UA  TRLER 

//  JOR 
//  DUP 

♦DELETE  DOOR 

/  /  FOR 

♦ONF  WORD  INTFGFRS 
♦ARITHMPTIC  IRACF 
♦TRANSFFR  TRACE 
♦IOCS! CARD. 11 82  PRINTFR) 

COMMON  PO ( 80 )  ,  T  AU ( 80 )  ,CPR(60)  .HR  E ( A )  ,HRW <  A ) , HRS ( A ) »  HFE ( A  )  , 

»HFW ( 4 ) *HFS<  4 ) .HRF  <  4 ) *HAT (4 ) »CRE  <4 ) *CRW( 4 ) *CRS( 4 ) <CDE  <  4 ) «CDW( 4 ) « 

♦  CDS ( A ) ,SBWE ( A )  , SRT ( A )  ,  SR  F ( A )  , SRW ( A )  , SRS ( A )  »OHF ( A )  ,OHT ( 4 )  .OFF  ( 4  )  , 

♦  CTC.  (4  )  »  TRW!  4  )  ,  CPW  (  a  )  ,TrcG<A)  ,TrwT  (a)  ,DORS(A  )  »WMAN  (  4  )  ,Pman(a)  , 

♦  OFT  ( A )  ,  I A(  1 1  )  ,K  .NOT  ,U  I  ( A  ,  A  ,  A )  , TC ( A , A  ,  A  )  , CS (  6  )  ,RGG(  A  )  ,RW ( A  )  ,CCG ( A  )  , 

♦  AIRF(a)  ,AIRF(  10)  ,TC2(A,A,A)  , TCS( A  .4  .4 )  »U I  2 ( A  » A  ,  A )  »  U I  8 ( A  ,  A  ,  4 ) 

C . 

C  IGLOO  DOOR  SUBROUTINE 

C . 

C . FIXFD  PARAMFTFR  VALUFS . 

A  =  1 ,/8. 

FFFT=1?0. 

AI=1?.R? 

C . DFTFRMINFS  NO.  OF  IGLOO  DOOR  PARAMFTPRS . 

DO  son  I =1 , A 
IF! DORS (  I)  )  901  .901  ,900 

SOI  K 1  =  I  - 1 

GO  TO  902 
900  CONTINUE 
K  1  =4 

C . 

C . CRITICAL  IMPULSF  FOR  IGLOO  DOOR . 

C . 

11  =  1 
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1 1 1  =  1 

902  DO  9 no  J=1 ,<l 

S I G  =DORS (J)*100O0» 

RFND=SIG*AI/2.71 
WORK=RFND*.3075 
S  =  1  .1 

H=SORT ( 2 • *1 • 1 *WORK ) 

900  UI  ( J, I  I ,  I  I  I ) =H*1 000./?8fl0. 

WRITE! 3*29) 

29  FORMAT ( 1 X,  ' CRIT ICAL  IMPULSF  FOR  IGLOO  DOOR'*/) 

WRITE! R,30) 

WRI TF( 3  *40  ) 

DO  90/,  j  =  ]  ,«1 
STG=DORS( J)*10OOO. 

904  WRITE!  3*  11  )  S I G  ,IJ  I  (  J  *  I  I  *  I  I  I  ) 

C . 

C . DETERMINES  MINIMUM  STAND-OFF  STORAGE  DISTANCE  FOR  FACH 

C . CHARGE  SIZE  AND  SET  OF  PARAMETER  VALUES . 

C . 

C . VARIATION  OF  PARAMFTFR  VALUFS . 

DO  905  I =1  * K 1 
C  I=!JI  (  I,  1 1  *  1 1 1  ) 

S IG=DORS ( I )*10000. 

C . VARIATION  OF  CHARGF  SIZF . 

WRITE! 3*31 )  SIG 
DO  906  L  =  1 »  5 

C . VARIATION  OF  SCALED  DISTANCE . 

DO  907  L 1 =1*50 

C . BLAST  PARAMFTFRS  FOR  SCALED  DISTANCF 

P I =P0( LI  ) 

PR=?.*PI*(102.9+4.*PI)/(10?.9+PI) 

DUR=TAU(L1 )*CS(L)**A 
IF ( Ll-37 )  917,917*918 

917  S=L1 

R= ( S+3 • )*CS(L)**A 
GO  TO  919 

918  S  =  L 1 

R  =  (  (S-37.  )*5.+40.  )  *  C  S  (  L  )  **A 

919  00=2.5* ( PI**2.  )/( 102. 9+P  I  ) 

VEL=1 187. *SORT ( 1 .+ ( 6. *P I ) / ! 1 02 .9+P I ) ) 

TS=  !3.*10./VEL)*1000. 

IF! DUR-TS 1920,920 *921 

C . CASE  I.  DURATION  OF  BLAST  LFSS  THAN  TIME  TO  STAG. 

920  CAL  I  =  . 5*PR*DUR 
GO  TO  922 

C . CASF  IT.TIMF  TO  STAG.  LFSS  THAN  DURATION 

921  PT  S=P I  * ( 1 .-TS/DUR)*EXP(-IS/DUR)+QO*( 1.-TS/DUR)*#2.#EXP(-2. 

**T  S/DUR ) 

DIV=(DUR-TS)/50. 

DO  1140  17=1*50 
S=I  7 

TFND=TS+S*DIV 

PFND  =  PI*  (  1  .  —  TFND/DUR  )*FXP(-TFND/DUR)  +00*  (  1  .  -  TFND/DIJR  )  **7  .  *  EXP  ( -2  . 
*  *TFND/DUR ) 

I F ( PEND- • 2 )  1141,1141,1140 

1140  CONTINUE 

1141  CAL  I  =  . 5* (PR-PTS )*TS+TS*PTS+. 5*PTS* ( TEND-TS) 
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922  I  r  (  r  I -cal  I  >  907,924,974 

C . MINIMUM  01  ST ANCF  FOR  GIVFN  CHARGE  SFLECTFD . 

924  WRITF(3,?R)  CS<L)»R 

WRITF<3,17)  PI»PR»DUR 
GO  TO  906 
907  CONTINUE 
906  CONTINUE 
905  CONTINUE 

CALL  LINK(CALLR) 

11  FORMAT ( 2F1 5.? ) 

17  FORMAT(5X,' INC.  OVPRES.=  •  ,E10.2, 'PFA<  RFFL . PRFS . =  ' , F 1 0 . 2 , 

* ' OUR • = • >F10,7 i /  /  I 

?B  Format ( 1 X. ' CHARGE  WFIGHT=' ,F9.0  ,  • D I STANCE= '  »F10.2 ,/ ) 

31  FORMAT ( IX MINIMUM  DISTANCE  FOR  STR  ESS  = ' »  F9  •  0  ,/) 

30  FORMAT ( IX MODULUS  OF  RUPTURE ', 1  OX ,' C  IMPULSF',  ) 

40  FORMAT ( 1 X LRS/SO  IN  '  » 1 5X , • PSl -MS • ,/ ) 

FND 

//  DUP 

♦STORE  WS  UA  DOOR 

//  JOB 
//  DUP 

♦DELETE  MAN 

//  FOR 

♦ONE  WORD  I  NT  EGFRS 
♦ARITHMETIC  TRACF 
♦TRANSFER  TRACF 
♦IOCS(CARD*I 132  PRINTFR) 

COMMON  PO  <50 ) »TAU( 50) »CPR ( 50) *HRF( 4 ) *HRW( 4 ) *HRS(4 ) tHFF  <  4 1 ♦ 

♦  HFW ( 4 ) »  HFS ( 4 ) »  HBF ( 4 )  ,  HBT  (4)  »C  RE ( 4  )  *  CR W ( 4 )  » CR S ( 4  )  »  C  DE ( 4  )  »  CDW ( 4 )  * 
*CDS<4 ) ♦ SrWE ( 4) ’SrT ( 4  )  *SrE  <4)  *Srw (4)  »SrS( 4)  *0HF( 4)  »0HT (4 ) ’OFF  <  4)  * 
♦CTG ( 4 ) *  TRW ( 4 )  »CPW ( 4 )  »  Tr  CG  (  4 )  »  TRWT (4 )  »DORS ( 4 ) » WM AN <  4 ) *PMAN (4)  * 

♦  OFT (  4  )  »  I  A  (  11)  »<»N0T»UI<4»4»4)  »  Tc  <  4  »4  »4  >  »CS<5)»BCG<4)’BW(4)»CCG<4)» 

♦  A  I RE ( 4 ) *  A I RF( 10 ) ♦ TC? (4*4*4) »TC3<4*4*4) *UI ?  <  4  *4  »4 )  *UI 3( 4»4*4) 

C . 

C  MAN  TRANSLATION  MOOFL 

C . 

C . FIXED  PARAMETER  VALUES . 

A  =  1 

I  11=1. 

C . DETERMINES  NO.  OF  MAN  TRANSLATION  PARAME  TFRS . 

DO  1000  1=1,4 

I F ( WM  AN (  I )  )  100  1,100  1,1000 

1001  <1=1-1 

GO  TO  1002 
1000  CONTINUE 
<1  =  4 

C . 

C . CRITICAL  IMPULSF  FOR  MAN  TRANSLATION . 

C . 

1002  DO  1003  1  =  1  , <1 
W=WMAN ( I ) 

DO  1003  J=1,<1 
P  A=  PMAN ( J) 

1  00  o  UK  I  ,  J,  I  I  I  )  =  (72000.*W/3B6.  )/(  PA*144.  ) 

WRITE(3,30) 

WRITE! 3, 31 ) 

W  R I T  E (  3,40) 
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DO  1004  I  - 1  »  K 1 
DO  1004  J  =  1  »  K  1 
W= WMAN ( I ) 

PA=PMAN ( J ) 

1004  WR I TE( 3*  3? )  W »P A *U I  (  I  »  J  »  1 1  I  ) 

C . 

C . DETERMINES  MINIMUM  STAND-OFF  STORAGE  DISTANCE  FOR  EACH 

C . CHARGE  SIZE  AND  SET  OF  PARAMETER  VALUES . 

C . 

C . VARIATION  OF  PARAMETER  VALUES . 

DO  1005  I  =  1  » K  1 
DO  1005  J  =  1  » K  1 

ci=ui(i*j*im 

W  =W  MAN  (  I  ) 

P A=PMAN ( J ) 

WR I TE(  3,  S3 )  W,PA 

C . VARIATION  OF  CHARGF  SIZF . 

DO  1006  L  =  1  »  5 
S0=0. 

C . VARIATION  OF  SCALED  DISTANCE . 

DO  1007  L 1  =  1 ,50 

C . BLAST  PARAMETERS  FOR  SCALED  DISTANCE 

P I=PO( LI ) 

Dl)R  =  T All  (  LI  )*CS(L)**A 
I F( Ll-37 )  1008,1008,1009 

1008  S=L1 

R  = ( S+3. )*CS(L)**A 
GO  TO  1010 

1009  S=L 1 

R  =  (  (S-37)*5.+40.)*CS!L)**A 

1010  00=?.5*( PI**?. ) /( 10?.9+PI ) 

P  R  =  ? • *P I  * ( 10?, 9  +  4 • *P I )/( 10? , 9  +  P I  ) 

. EAR  DRUM  RUPTURF  IS  I NCURRFD  IF  REFLFCTFD  PRFSSURF  IS 

. GRFATFR  THAN  5.  PSI . 

I F ( PR-5, )  1011, 101?, 101? 

101?  R ANGE=R 

GO  TO  1007 

1011  S=QO*DUR*. 5 

IF(S-CI)  1013,1013,1014 

1014  SO=R 

GO  TO  1007 

1013  I F ( SO )  1015,1015,1016 

C . MINIMUM  DISTANCE  FOR  GIVEN  CHARGE  SELECTED . 

1015  WR I T E ( 3,34)  CS  (  L ) *  RANGE 
GO  TO  1006 

C . MAXIMUM  TRANSLATION  DISTANCE . 

1 C 1 6  DIST  =  .7*3. 5*DUR / 1 000 , 

C . MAXIMUM  ACCELERATION . 

ACC=00*PA/ (W/386. ) 

WRITE! 3,36)  CS(L)»RANGF 
WRITE! 3,35)  S0,DIST,ACC 
GO  TO  1006 
1007  CONTINUE 
1006  CONTINUE 

1005  CONTINUE 

CALL  L  INK ( CALLR ) 

30  FORMAT! IX, 'CRITICAL  IMPULSF  FOR  MAN'/) 
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31  FORMAT! IX. 'MAN  WF I GHT  =  • , 1  OX , • MA N  P  AREA= ' *  1  OX , • C  IMPULSE='*  ) 

AC  FORMAT (] X. • LBS* , 15X *• SO  IN'  *15X, 'PSI-MS*  */) 

3?  FORMAT (  3  F 1 0 . 2 ) 

33  FORMAT! IX* 'MINIMUM  DISTANCE  FOR  MAN  WGT = • * F 1 0 . 2 ♦ 1  OX  * ' P  ARF  A= 1  * 

*F 10 . 2  * / ) 

3 A  FORMAT ( IX* • FOR  C  SIZF='*F9.0  *'ONLY  FAR  DAMAGE  AT  R = ' * F  1 0 . 2 » '  I S  •  , 

*  '  I NCIJRRFD '  *  /  ) 

35  FORMAT!  5X * ' TRANSLA T I  ON  AT  R = • . F 1 0 . 2  *  IX  * • T  D I  ST = i , F 1 0 . 2 * ' ACC  =  '  * 
♦F10.2,/) 

36  FORMAT! IX* 'FOR  C  SIZE='*F9.0  * • EAR  DAMAGE  AT'.FIO.?*) 

END 

//  DUP 

♦STORF  WS  UA  MAN 

//  JOR 
//  DUP 

♦DFL  F  T  E  AIRP 

//  FOR 

*ONE  WORD  INTEGERS 
♦ARITHMETIC  TRACE 
♦TRANSFFR  TRACE 
♦IOCS ( CARD. 11 32  PRINTFR) 

COMMON  PO(50)*TaU<S0)*CPR(50),HRE<A)*HRW<A)*HRS(A)*HFF(A)* 

♦  HFW ( A ) *HFS<  A ) *  HRF ( A )  *HRT<A)’CRF<A) *CRW ( a ) *CRS( A)  *CDE<  A )  *CDW( A )  * 
*roS ( A) *  SRWE ( A ) *  SrT ( A )  *  SRF ( A)  *  SRW ( A )  *SRS ( A)  *OHF ( A)  *OHT ( A  )  *OFF  <  A  )  * 
*CTG(A)*TRW(A)*CPW(A)*TRCG(A)*TRWT(A) *DORS ( A ) *WMAN( A) »PMAN (A) * 
*0FT(A)*IA<11)*K ’NOT  *Ui<a*A’A)’Tc(A’A’A)’CS<5) ’BCG<  A ) »BW( a ) *CCG<  A ) * 
♦AIRF<A)*AIRF<10)’TC2<A*A*A)*TC3<A*A*A)*UI2<A*A*A)*UI3<A*A*A) 

C  AIRCRAFT  SUPROUTINF 

A  =  I  .  /  3  . 

A2=.3DAA 
A  1  =  .3079 
G  =  R  86 • 

DO  1100  1=1, A 
I F (  A IR  E (  I  )  )  110  1*110  1*1100 

1101  K 1 = I  —  1 

GO  TO  1102 
1100  CON  TIN  UF 
<1  =A 

1102  WRITF(3,12) 

WRITE! 3, 13 ) 

WRITE! 3, AO) 

DO  1103  1=1, <1 
E=A IRE ( I )»1000000. 

P  =  36  .  +  (  I -1 ) *2. 1 3 

1103  WRITF! 3, 1A)  F*P 
DO  110a  1=1, <1 

F  =A  TRF (  I )»1000000. 

DO  1105  L  =  1  *  5 

DO  1106  L  1  =  1  *50 

I  F ( L  1-37 )  1111,1111,1112 

mi  s=li 

R  = ( S  +  3 • )*CS!L)**A 
GO  TO  1113 

1112  S=L 1 

R= (  ( S-3  7.  )*5.  +  A0.  )  *CS ( L  )  *♦ A 

1113  PI=PO( LI  ) 

DUR  =T  AU ! LI ) *CS ! L ) **A 
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P=06.+( 1-1 )*?. 13 
F=DUR/P 

I F ( F-.5)  1114,11)4,1115 

111^  FAC=1. 

GO  TO  1109 

1115  I F ( F-5 • )  1107,1108,1108 

1108  F  AC  = 1,81 
GO  TO  1109 

1107  S  =  F /  •  5 
M  =  S 

DIF=5-M 

FAC=AIRF(M)+DIF*(AIRF(M+1 )-AIRF(M) ) 

1109  PRFS=7.20/FAC 

I F ( PRFS-CPR ( L 1 1 )  1106,1110,1110 

1110  WRITFt  0,18)  CS ( L )  ,R 

18  FORMAT ( IX, '  CHARGF  WGT= • , F 10 . 2 , 10X , • DI STANCE= • , F 10 . 7 » > 

WR  I  TE ( 0,17)  PI  ,CPR ( LI ) ,DUR 
GO  TO  1105 
1106  C  ON  TIN  UF 
1105  GONT  INI  IF 
1104  GONTIMUF 

GALL  LINK(GALLR) 

12  FORMAT(lX,'FULL  NATURAL  PFRIOD  OF  FRAMF  MFMRFR*,/) 

10  FORMAT(lX,'E* ,16X,'FULL  PER  I ODS ( MS )  '  , ) 

40  FORMAT ( 1 X, '  LRS/SQ  IN' » 1 5X  * • MSEC '  , / ) 

14  FORMAT (F10.0,F8#0) 

17  F ORM  AT ( 6 X , '  INC.  OVPR FR .  =  •  , F 1 0 . 2  ,  •  PF A K  RFFL . PRFS .  =  '  , F 1 0 . 2  , 
*  •  r>!JR.  =  '  ,  FI  0. 7  ,  /  ) 

F  NO 

//  DUP 

*  S  T  0  R  F  WS  UA  AIRP 
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Computer  programmed  methodology  has  been  formulated  which  may  be 
used  to  estimate  minimum  separation  distances  for  the  storage  of  five 
different  quantities  of  stored  explosives  with  respect  to  minimizing 
damage  to  ten  specific  types  of  civilian  targets.  These  procedures  have 
been  based  on  the  results  of  more  recent  analytical  and  experimental 
blast  damage  assessment  programs  in  which  structures  or  structural  ele¬ 
ments  subject  to  blast  could  be  related  to  the  civilian  targets  of 
interest  through  similar  structural  elements  and  construction  techniques. 

Analytical  techniques  have  been  programmed  which  will  develop  a 
minimum  separation  distance,  for  a  selected  target  and  quantity  of 
explosives,  which  minimizes  the  probable  risk  that  the  blast  damage 
resulting  from  an  accidental  detonation  will  exceed  a  predetermined 
acceptable  level.  A  discussion  of  the  program  and  operating  instruc¬ 
tions  are  presented  in  this  manual. 
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